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LASER  ANNEALING  OF  ION  IMPLANTED  SEMICONDUCTORS. 
MECHANISM  OF  AMORPHOUS  CRYSTALLINE  TRANSITION. 


1.  INTRODUCTION 


One  of  the  most  fascinating  debates  in  the  field  of  Semiconduc¬ 
tors  Physics  in  recent  years,  centers  on  the  fundamental  interpre¬ 
tation  of  laser  annealing.  The  question  of  the  energy  transfer 
from  an  intense  beam  to  a  disordered  material,  such  as  amorphous 
silicon,  resulting  in  the  crystallization  of  the  amorphous  substan¬ 
ce  has  been  approached  from  two  different  points  of  view,  both 
referring  to  a  set  of  fairly  clear  experimental  results.  On  one 
hand,  a  claim  has  been  made  that  the  laser  beam  simply  heats  the 
sample  up  to  melting  the  amorphous  material  which  on  cooling  crys¬ 
tallizes  from  melt^  ^  .  On  the  other  hand,  one  has  considered  that 
amorphous  to  crystalline  phase  transition  can  occur  at  low  tempera¬ 
tures  without  passing  through  the  molten  stated.  Many  arguments 
have  been  developped  in  support  of  these  interpretations.  A  large 
number  of  publications  have  appeared  in  the  litterature. 


The  results  we  are  reporting  here  constitute  a  contribution 
toward  the  elucidation  of  the  mechanism  of  amorphous  to  crystalline 
transition.  Our  effort  has  been  focussed  on  the  two  extreme  time 


scales  :  infinite  time  observation  of  the  structural  modifications 


which  have  been  retained  after  laser  annealing  and  very  short  time  : 
the  evolution  of  the  system  in  the  femtosecond  range.  We  have  attemp¬ 
ted  to  correlate  the  steady  state  observations  with  the  sub-picose¬ 
cond  modification  of  the  system.  If  the  general  conclusion  is  that 
random  nucleation  in  the  high  density  plasma  fluid  state  governs  the 
final  crystallization  of  the  system,  we  are  using  essential  links 
to  demonstrate  such  an  opinion.  The  description  of  the  system  in 
the  highly  excited  non  linear  regime  maintained  during  the  irradiation 
of  100  fs  is  not  available.  The  description  of  the  relaxation  pro¬ 
cesses  in  the  fs  range  is  not  clearly  established.  The  relation 
between  the  excited  short  living  fluid  and  melting  is  not  known. 

The  theory  of  melting  and  recrystallization  in  this  regime  is  not 
establ i shed. 
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We  are  presenting  here  a  gread  deal  of  interesting  and  important 
results.  They  are  contributing  to  ask  the  right  questions  more  than 
they  have  solved  the  essential  problems.  As  a  result  of  this  effort 
nevertheless,  we  are  convinced  today  that  the  radiation  energy  in¬ 
troduced  into  the  condensed  matter  is  retained  in  a  highly  excited 
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phase  with  more  than  10  antibonding  states  as  long  as  the  radia¬ 
tion  field  continues  to  arrive  on  the  system.  The  evolution  of  this 
particular  fluid  is  the  interesting  question  with  which  one  should 
be  concerned  now. 

In  addition  to  the  fundamental  interest,  laser  annealing  has 
been  considered  to  have  a  strong  potential  in  the  technology  of 
semiconductor  doping  by  ion  implantation.  For  all  these  reasons, 
the  attention  of  a  large  audience  has  been  focussed  on  this  problem 

We  believe  that  the  results,  we  are  presenting  here  show  the 
possibilities  and  set  the  limits  of  the  use  of  laser  annealing  in 
the  semiconductor  technology. 


2.  RAMAN  SPECTROSCOPY  OF  GALLIUM  ARSENIDE  AMORPHOUS  TO  CRYSTALLINE 
TRANSITION  INDUCED  BY  LASER  ANNEALING . 

Raman  spectroscopy  of  laser  annealed  semiconductors  rendered 
amorphous  by  ion  implantation  gives  the  possibility  to  follow  the 
rec rys tal 1 i zat ion  of  the  material.  This  technique  which  has  succes- 
fully  been  applied  to  silicon1  *  '  can  as  well  be  used  to  characte¬ 
rize  ionic  semiconductors  as  Gallium  Arsenide. 

2.1.  Experimental  techniques. 

Semi  insulating  GaAs  have  been  implanted  and  laser  annealed  at 

the  Eaboratoiro  Central  de  Recherches  de  la  Thomson  C.S.F.  (Corbe- 

ville).  The  dose  chosen  was  2.10^As/cm^  +  2.10^  Ga/cm^  under 

300kcV  in  order  to  keep  a  correct  stoechiometry  in  the  samples. 

Laser  annealing  was  achieved  by  using  a  YAG  laser  equipped  with  a 

frequency  doubler.  This  laser  deliveres  energy  pulses  in  the  energy 

2 

range  0.3  to  0.6J/cm  at  a  wavelength  of  530nm.  Raman  spectra  were 
recorded  by  using  the  488nm  line  of  an  Argon  laser.  The  diameter 
of  the  spot  of  the  Argon  laser  was  of  the  order  of  80  m.  This  size, 
compared  to  the  1mm  diameter  spot  of  the  YAG  laser  permits  a  rather 
good  Raman  characterization  of  all  the  surface  of  the  annealed 


3 


region  by  sweeping  the  Raman  laser  spot  inside  the  annealed 
zone . 


2.2.  Experimental  results. 


Figures  1  and  2  present  the  Raman  spectra  of  two  GaAs  samples 
in  a  backscattering  geometry  on  faces  (100)  and  (110).  On  these 
two  different  faces,  the  selection  rules  are  complementary,  the  TO 
phonon  being  forbidden  on  the  (100)  face,  the  LO  phonon  forbidden 
on  the  (110)  face.  For  perfectly  crystalline  samples,  the  selection 
rules  are  verified,  the  residual  forbidden  peaks  in  each  case  are 
due  to  residual  misalignment.  The  frequencies  of  the  lines  are 
respectively  =  268  cm-1,  win  =  291  cm' \  their  widthes,  measured 


LO 


is  equal  to  3.6  cm1  for 

A  A  /~) 

both  peaks.  After  an  implantation  of  2.10  ion  Gallium  per  cnr  and 


with  an  instrumental  resolution  of  2  cm 

1  P 
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2.10  ion  Arsenic  per  cm  ,  the  samples  are  rendered  amorphous  and 
the  spectra  observed  are  typical  for  this  kind  of  amorphous  material 
^  i.e.  a  density  of  states  spectrum. 


After  laser  annealing,  the  spectra  are  drastically  modified. 

For  a  low  energy  density,  the  forbidden  peaks  are  respectively 
activated  for  each  face  ;  this  is  especially  spectacular  for  the  (100) 
face  for  which  the  forbidden  TO  is  larger  than  the  allowed  LO.  More¬ 
over,  the  frequencies  of  these  peaks  are  shifted  towards  the  low 
frequency  region  and  their  widthes  are  broadened  to  more  than  6  cm  1 . 
As  the  annealing  energy  is  increased,  the  forbidden  peaks  intensities 
tend  to  decrease,  their  frequency  shift  and  their  widthes  reduce. 

However,  the  normal  characteristic  are  not  completely  recovered  at 
2 

0.6J/cm  especially  for  the  (100)  face. 


The  spectacular  behaviour  of  the  (100)  face  for  low  energy  annea¬ 
ling  has  lead  us  to  study  the  Raman  spectra  at  the  periphery  of  the 
annealed  zone  by  sweeping  the  Raman  laser  spot  across  the  limit  of 
the  annealed  region  where  the  irradiation  energy  continuously  varies. 

As  can  be  seen  in  Figure  3,  the  spectra  are  continuously  modified, 
the  amorphous  band  at  250  cm" ^shifts  and  get  structured  changing 
into  the  TO  phonon.  The  allowed  LO  phonon  appears  only  as  a  shoulder 
on  this  band.  The  polarization  analysis  has  been  made  on  the  (100) 
face  sample,  with  the  incident  electric  field  parallel  to  the  <01 1> 
axe,  the  scattered  field  been  either  parallel  or  perpendicular  to 


the  incident  field.  For  this  last  polarization  both  peaks  TO  and 
LO  are  forbidden  in  a  perfectly  crystalline  sample  whereas  the 
LO  phonon  normally  appears  when  both  fields  are  parallel.  In 

Figure  4,  it  can  be  seen  that  the  selection  rules  are  well  observed 

• 

when  the  analysis  is  made  at  the  center  of  the  annealed  zone  but 
are  not  at  the  intermediate  region,  where  spectra  are  rather  iden¬ 
tical  for  both  polarization. 


The  more  important  feature  presented  here  seems  to  be  the  gradual 
transition  from  the  amorphous  spectrum  to  a  crystalline  one.  This 
fact  involves  a  continuous  variation  of  the  degree  of  order  versus 
the  annealing  density  energy.  The  late  apparition  of  the  longitudi¬ 
nal  optical  phonon  when  increasing  the  annealing  energy  confort  this 
hypothesis  as  the  existence  of  this  phonon  is  due  to  long  range 
Coulombian  forces.  A  way  of  understanding  such  a  continuous  varia¬ 
tion  is  to  consider  the  presence  of  microcrystallites  in  the  annealed 
region.  Such  crystallites  have  yet  been  observed  in  laser  annealed 
Silicon  and  Gallium  Arsenide ^ ^ .  There  are  very  few  published 
results  on  Raman  spectroscopy  of  small  crystals.  The  only  extensi¬ 


vely  studied  case  is  the  graphite 


(11) 


where  only  very  small  size 


crystallites  are  available.  In  this  material,  new  peaks  as  well  as 
a  shift  of  the  high  frequency  E2g  first  order  peak  are  reported  and 


attributed  to  the  wave  vector  selection  rules  relaxation.  The  shift 
observed  is  toward  the  high  frequency  region.  This  is  not  contradic¬ 
tory  with  our  observation  of  a  low  frequency  shift  in  GaAs  if  the 
special  shape  of  the  dispersion  curves  of  graphite  is  considered. 


In  this  material,  the  highest  branch  of  the  dispersion  curves  reaches 


-1. 


a  region  of  higher  frequency  1 600  cm  ),  than  the  zone  center  E 


2g 


mode  (1581  cm  ) ^  Furthermore,  this  E?  mode  is  transformed 

'■S  _  1 

into  a  band  situated  at  higher  frequency  (1600  cm  )  in  the  amorphous 


state  of  graphite  (glassy  Carbon)  which  is  an  opposite  behaviour  as 
compared  to  GaAs  or  Si. 


Only  a  complete  calculation  of  the  lattice  dynamic  of  a  micro¬ 
crystallite  would  fully  justify  the  hypothesis  of  the  size  effect. 

A  similar  computation  has  already  be  done  at  the  laboratory  on 
Silicon^.  The  extension  of  this  model  to  GaAs  is  not  straight¬ 
forward  because  of  the  presence  of  the  long  range  Coulombian  forces 
which  ensure  the  stability  of  this  ionic  crystal.  This  extension  is 
the  object  of  the  following  chapter. 


WAVE  NUMBER  [cm-i] 


Raman  spectra  of  a  laser  annealed  sample  of  GaAs  (100) 
from  various  points  of  the  sample.  The  incident  field 
is  parallel  to  a<100  ax^j  the  scattered  field  is  perpen 
dicular  or  parallel  to  the  incident  field. 


3.  THE  VIBRATIONS  OF  A  GALLIUM  ARSENIDE  MICROCRISTALLITE . 


3.1.  Lattice  dynamics  of  thin  ionic  slab  of  GaAs . 


The  complete  calculation  is  presented  in  three  publications 
added  in  annexe  (A1).  We  just  summarize  here  the  method  and  the 
results. 


A  full  lattice  dynamics  of  a  microcrystallite  is  outside  our 
possibilities  as  it  should  need  the  exact  knowledge  of  the  geome¬ 
tric  shape  of  the  grain  and  as  it  would  need  the  diagonalization 
of  matrix  of  order  3N  where  N  is  the  number  of  atoms  inside  the 
crystallite.  Even  for  very  small  crystallites  N  exceeds  some 
thousands  so  that  computation  duration  should  be  excessive.  That  is 
the  reason  for  which  the  crystallite  has  been  modelized  by  a  thin 
layer  of  material  inifinite  along  two  dimensions  and  finite  along 
the  third.  The  finite  dimension  should  be  sufficient  to  introduce  a 
shift  on  the  frequencies  of  vibration  as  compared  to  a  crystal  infi¬ 
nite  along  3  dimensions. 


To  resolve  this  problem,  it  is  necessary  to  solve  the  equations 
of  motion  of  the  atoms  of  the  slab,  that  is  to  say  to  diagonalize 
the  corresponding  dynamical  matrix.  To  simplify  the  problem,  the 
choice  of  a  cristallographic  orientation  of  the  slab  is  of  premium 
importance  as  symetry  consideration  can  reduce  the  computation.  In 
our  case,  the  slab  is  parallel  to  the  plane  (111).  With  classical 
notations,  the  equation  of  motion  are  : 
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The  slab  being  infinite  along  two  dimensions,  the  cyclic  condition 
of  Born  Van  Karman  can  be  applied  along  these  two  directions.  The 
solution  of  equation  (1)  can  be  expressed  in  the  form  : 


exp{-iu)t+2iriy  .x(lk)  } 
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where  the  vector  v(l3,  K)  depends  only  on  lj  the  parameter  which 
labels  the  cells  along  the  limited  dimension. 

y  is  a  two  dimensional  vector  of  the  two  dimensions  B?‘llouin 
zone  associated  to  the  structure.  By  reporting  (2)  into  (1),  we 
obtain  : 


'.(‘j")  'W3k:1jk,h’Vl3k,> 


with 
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Da 3  is  the  dynamical  matrix,  0^  being  the  matrix  of  force 
constants . 


Since  cyclic  boundary  condition  does  not  apply  to  the  limited 
dimension,  no  further  reduction  is  possible  except  those  associa¬ 
ted  to  symmetry  considerations.  Under  these  consideration  the 
3x3  matrixes  D  ( 1 ^ »  K,  L'^,  K')  which  describe  the  interactions 
between  the  plane  lattices  ake  the  form  : 


D'Clj,  K,  L'3,  K')  = 


A-B  0 
0  A+2B 


which  shows  that  the  modes  of  vibration  can  be  separated  in  two 
groups,  one  doubly  degenerated  which  describes  the  vibrations 
parallel  to  the  slab,  the  other  non  degenerated  which  describes 
the  vibrations  perpendicular  to  the  slab. 

To  go  further,  we  need  now  to  choose  a  model  of  force  constants. 
The  long  range  interactions  and  the  short  range  interactions  must 
be  processed  differently  as  it  is  well  known  that  Coulombian  forces 
introduce  discontinuity. 


,-.\s 

_»  «  .  •  «  »  ' 


« -  -  •'.*  >; 


>  W*  - 

•>  v  •-»  • 

.V, 

>  ,  .  . 


-/  %•  •/  • 
.•  v  v  vv  i 

•S-'.vSvv’ 


/  / 


7 


For  short  range  interaction,  the  model  choosen  is  the  Valence 
Force  Field  force  constants  with  10  parameters  describing  the 
linear  and  angular  bondings  between  first  and  second  neighbours. 
The  potential  energy  for  this  model  take  the  form  : 
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E  (Ar^)  +r.  p 
Ga-As  3  ' 


Ga-Ga 


Arik  + 
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4rij  4rjk- 


where  A,  p,  y,  a,  v,  kQ,  k^,  krr  are  the  ten  parameters  and 
rQ  and  r1  the  first  and  second  nearest  neighbour  distances. 

3.2.  Results  and  discussion. 

The  frequencies  of  vibrational  modes  of  a  thin  slabs  have  been 
calculated  for  a  thickness  varying  from  2  to  25  unity  cells. 

The  Figure  5  shows  the  variation  of  frequency  of  six  modes 
versus  the  thickness  of  slabs.  These  modes  are  grouped  in  two 
optical  modes  (TO  and  LO)  and  four  surface  modes.  The  two  opti¬ 
cal  modes  are  the  more  interesting  as  they  are  the  modes  that 
correspond  to  the  Raman  active  TO  and  LO  modes  in  the  limit  of 
large  number  of  unit  cells  (infinite  crystal).  As  can  be  seen, 
the  frequencies  of  these  modes  tend  very  rapidly  to  their  limits. 
After  10  unit  cells  the  frequencies  have  reached  a  stable  value 
corresponding  to  the  frequencies  of  the  infinite  crystal.  For  a 
small  number  of  cells  than  ten,  the  frequencies  of  the  TO  and  LO 
modes  are  inferior  to  the  bulk  frequencies,  result  which  corre¬ 
lates  the  Raman  spectra  observed.  The  surface  modes  are  not  easily 
seen  by  Raman  spectroscopy  so  that  their  interest  is  only  theoric 
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Frequency  variation  of  vibrational  inodes  versus  the 
thickness  of  the  slab. 
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as  they  cannot  be  compared  to  experimental  results.  Another  inte¬ 
resting  theoretical  deduction  is  that  for  a  given  slab  of  definite 
thickness,  all  the  modes  except  the  surface  ones,  fall  on  the 
dispersion  curves  of  the  infinite  structure  of  GaAs  (Fig. 6).  The 
four  missing  point  at  the  edge  of  the  Brillouin  zone  correspond 
to  the  four  surface  modes.  These  surface  modes  can  thus  be  consi¬ 
dered  as  a  superposition  of  the  four  missing  modes.  This  fact 
implies  that  the  continuum  of  frequencies  obtained  by  several 
authors  when  calculating  solution  for  different  wavevectors  on  the 
plane  of  the  slab,  is  the  projection  of  the  dispersion  surface  on 
the  two  dimensional  Brillouin  zone. 


4.  EXPERIMENTAL  DETERMINATION  OF  THE  TEMPERATURE  OF  SILICON  AT 
THE  SPOT  OF  A  CONTINUOUS  LASER. 


The  possibility  to  characterize  the  laser  annealing  during  the 

impulsion  by  using  the  Raman  spectroscopy  was  evidenced  by  Compaan 
f  1  31 

et  al/  Their  measures  showed  that  the  temperature  reached 

immediately  after  the  end  of  the  pulse  was  of  the  order  of  40fl°C. 
This  result,  which  is  still  controversed ,  is  the  last  experimental 
evidence  for  a  non  thermal  annealing  during  manosecond  scale.  To 
test  the  validity  of  this  kind  of  measurements,  we  have  measured 
the  temperature  of  a  Silicon  sample  during  a  continuous  laser  irra¬ 
diation  by  using  Raman  spectroscopy.  Three  quantities  have  been 
studied  :  the  Stokes  to  anti-Stokes  ratio,  the  shift  of  the  normal 
mode  and  the  broadening  of  the  normal  mode . 


4.1.  The  Stokes  to  anti-Stokes  ratio. 


Raman  scattering  consists  of  an  inelastic  scattering  of  the 
light  by  quantified  excitations.  In  the  case  of  a  semiconductor, 
the  elementary  excitations  responsible  of  the  scattering  are  the 
phonons.  Two  processes  coexist  in  the  scattering. 


Tn  the  Stokes  process,  incident  photons  of  frequency  u>.  create 
phonons  of  frequency  fi  and  are  scattered  with  a  frequency  w 


In  the  anti-Stokes  process 
ar.ihilate  the  phonons  and  are 


the  incident  photons  interact  and 
scattered  with  a  frequency  u)^s=u>^+S3. 
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Raman  spectra  are  thus  formed  by  a  set  of  peaks  symetrically 
disposed  in  comparison  with  the  incident  frequency  .  The  inten¬ 
sity  ratio  of  these  two  processes  is  governed  by  the  statistic 
of  Bose-Einstein  : 

I  .  , 

=  £_J.  where  n  =  1  /  (exp  [Kft/kT  |  - 1 ) 
lAs  n 


In  fact,  this  term  only  describes  the  thermodynamical  part  of  the 
intensity  ratio  and  in  order  to  be  accurate,  it  is  necessary  to 
take  into  account  some  corrective  terms. 


First,  the  Raman  effect  is  a  scattering  phenomenon  and  then, 
the  scattered  intensity  varies  as  the  fourth  power  of  the  frequency 
of  the  incident  light.  More  precisely,  it  varies  as  the  product 

3 

u>i  .  <ijs  where  and  u>s  are  the  incident  and  scattered  frequencies 
(uj  *  uj.  -  n.,  fl  being  the  phonon  frequency)  .  This  term  introduces 

si  3 

the  corrective  factor  (“s/wA  )  which  is  temperature  independent. 


Another  dependence  upon  the  incident  frequency  is  due  to  a 

resonance  effect.  Raman  scattering  intensity  increases  as  the 

incident  light  comes  close  to  the  direct  gap  of  the  material. 

This  resonance  term  is  difficult  to  compute  so  we  have  chosen  to 

use  experimental  results  obtained  at  the  laboratory ^ .  This  correc- 
S(Wi  ws) 

tive  factor  - ■gj- ■  * ^ — y  is  temperature  dependent  by  the  intermediary 
of  the  energy  ga^’v^fiation.  However,  this  temperature  dependence 
can  be  neglected  when  the  energy  gap  is  far  from  the  energy  of 
incident  wavelength  which  is  the  case  in  our  conditions. 


A  third  correction  is  due  to  the  difference  between  the  values 
of  the  absorption  coefficient  at  the  incident  and  scattered  energy. 
This  gives  a  correction  of  the  form  : 


ai  *  ^s 
ai  +  aAs 

Besides  all  those  theoretical  corrections,  these  is  an  experi¬ 
mental  correction  induced  by  the  variation  of  the  overall  response 
of  the  spectrometer.  This  factor  K  can  be  measured  with  a  standard 
spectral  lamp.  With  all  these  corrective  factors,  the  ratio  of 


“i  *  “As  “s  .3  s(“i.V 
“i  *  “s  “As  '  S(“i'“As' 


In  the  case  of  Silicon,  the  entire  correction  has  been 
computed  and  for  an  energy  of  of  20.492  cm”1  is  equal  to  0.83, 


0.83  =0.83  exp  (Jlfi/kT) 


A  measure  of  the  Raman  intensities  of  Stokes  and  anti-Stokes 
peaks  versus  temperature  has  been  realized  on  a  Silicon  heated 
in  an  oven.  The  points  on  Figure  7  represent  the  ratio  of  I.  /Ig. 
On  the  same  figure  is  represented  by  a  continuous  line  the  theore¬ 
tical  function  I^s/Is  =  1»2  exp  (-74.8. 3/T)  .  As  can  be  seen  on  the 
figure,  the  fit  is  only  good  at  low  temperature.  At  high  tempera¬ 
ture,  this  curve  is  significantly  lower  than  the  experimental 
points.  The  reasons  of  this  discrepancy  is  the  incertitude  on  some 
of  the  coefficients  which  form  the  correction  stated  above.  In 
order  to  increase  the  precision  of  the  measure,  we  have  choosen 
to  realize  a  least  square  best  fit  of  the  experimental  results 
with  a  theoretical  curve  of  the  form  : 


W1*  =  A  exp  (-748. 3/T]  where  A  is  the  parameter  to  be  adjus¬ 
ted.  This  computation  gives  a  value  for  A  of  1.284  so  that  we  can 
practically  measure  the  temperature  by  using  the  formula  .* 


lXs/ls  =  1.284  exp  (-748. 3/T) 

This  curve  is  represented  on  Figure  7  by  the  interrupted  line. 
The  temperature  obtained  by  this  method  is  exact  with  a  precision 
of  about  20  % . 

4.2.  Anharmonic  properties  of  Silicon. 


The  ratio  between  the  Stokes  and  anti-Stokes  intensities  is 
not  the  only  quantity  which  can  be  used  to  determine  the  temperature. 
Anharmonic  properties  which  involve  a  shift  and  a  widening  of  the 
Raman  peaks  an  be  also  used.  These  effects  are  due  to  interaction 
between  phonons.  The  simpler  model(15»  1 6)  supposes  that  the  optical 
phonon  anihilates  giving  rise  to  two  acoustical  phonons  with  half 
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the  frequency  of  the  optical  mode.  In  this  case,  the  width  of 
the  optical  phonon  responsible  of  the  Raman  peak  can  be  evaluated 
to  : 


r(T)  =  r(0)  |i+  \  (i) 

e  - 1 

The  damping  constant  and  the  frequency  shift  have  been  inves¬ 
tigated  systematically  as  a  function  of  temperature.  Figure  8 
gives  the  temperature  variation  of  the  damping  constant  T(T) 
between  5  and  HOOK.  The  dashed  curve  represents  r(T)  calculated 
from  the  relation. 

This  equation  (1)  is  an  approximate  expression  for  the  tempera¬ 
ture  dependence  of  the  damping  constant  based  on  three-phonon  pro¬ 
cesses  (cubic  anharmoni city  in  second  order)  and  the  simple  Klemens 
model*-16).  It  seriously  underestimates  the  damping  constant  at  high 
temperatures.  We  attribute  this  discrepancy  at  least  in  part  to  the 
neglect  of  four-phonon  processes. 

It  is  of  interest  to  investigate  whether  this  discrepancy  can 
be  eliminated  by  generalizing  Equation  1  to  include  the  contribution 
of  four-phonon  processes.  Following  the  approach  of  Klemens*-11),  we 
write  the  kinetic  equation  for  the  net  rate  of  decay  of  an  incident 
phonon  into  three  thermal  phonons  in  the  form  : 

(6n0)  =  -B  |  (6n0  +  n0)  (n^l)  (n2+1)  (n3+1)  ^ 

-  (6nQ  +  n3+1 )n1n2n3 |  , 

where  5nQ  is  the  deviation  of  the  incident  phonon  occupation  number 
from  its  thermal  equilibrium  value  np  and  B  is  a  constant.  Using 
the  equilibrium  condition  : 


n0(n1+l)(n2+1)(n3+1)  -  (n0+1)n1n2n3 


0, 


(3) 
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We  can  rewrite  Eq.(3.5)  as  : 

^  (6nQ)  =  B(n1n2+n1n3+n2n3+n1+n2 

CD 

+n3+1) 6nQ. 


Energy  conservation  can  be  satisfied  in  the  simple  Klemens 
fashion  by  setting^  =u,2  =  ^3  »  ***0/3.  Consequently,  n^=n2=n3- 
The  generalization  of  Eq.(3.4)  to  four-phonon  processes  then 
takes  the  form  : 


r(T)=A  (  1+  — -2—  ) 
'  ex  -  1  ' 


(ey-1) 


(5) 


Where  y  =  Wv»q/3KbT  and  A  and  B  are  constants.  In  the  high-tempe¬ 
rature  limit,  the  factors  multiplying  A  and  B  in  Eq. (5)  vary  as 
7 

T  and  T  ,  respectively. 


Equation  (5)  has  been  used  to  fit  the  experimental  data  presen¬ 
ted  in  Figure  8  by  suitably  choosing  the  constants  A  and  B.  The 
best  values  of  A  and  B  are  found  to  be  1.295  and  0.105  cm  \  res¬ 
pectively,  and  the  resulting  plot  of  T(T)  vs  T  is  given  by  the  solid 
curve  in  Figure  8.  We  see  that  the  agreement  between  the  calculated 
curve  and  the  experimental  points  is  now  quite  good. 

The  experimental  results  for  the  line  position  f2(T)  as  a  func¬ 
tion  of  T  are  shown  in  Figure  9.  Also  shown  is  the  fit  to  the  data 
(solid  curve)  specified  by  the  expressions  : 


ft(T)  =  w0  ♦  ACT) 


(6) 


and 


A(T)  =  C  (  1  + 


ex  -  1 


iy  -  1  (ey- 1 ) 


)  ,  (?) 


Where  C  and  D  are  constants  with  the  values  528,  -  2.96,  and 
-  0.174  cm  ,  respectively.  Equation  (7)  is  the  analog  of  Equation 
(5)  and  specifies  the  contributions  of  three-phonon  and  four-phonon 


processes  to  the  frequency  shift.  The  agreement  between  the  expe¬ 
rimental  points  and  the  solid  curve  is  seen  to  be  good. 

I f  we  try  to  fit  the  experimental  data  with  three-phonon 
processes  only  by  omitting  the  'term  in  Eq.(3.10)  with  the  factor 

D,  we  obtain  the  dashed  curve  in  Figure  9  with  n  =  529  cm  ^  and 

- 1  u 
C  =  4.24  cm  .  Although  this  curve  fits  the  data  well  at  tempera¬ 
tures  up  to  600K,  it  is  clearly  inadequate  at  higher  temperatures. 
This  demonstrates  the  necessity  of  including  terms  corresponding 
to  four-phonon  processes  in  the  expression  for  A(T) . 

In  principle,  the  four-phonon  contributions  in  Equations  (5) 

and  (7)  should  include  terms  arising  from  difference  processes. 

We  have  omitted  such  terms  on  the  grounds  that  their  inclusion 

2 

would  simply  introduce  additional  terms  varying  as  T  and  as  T 

in  the  high-temperature  limit  and  would  not  add  any  new  qualitative 

features. 

The  temperature  determined  by  measuring  the  width  and  the  frequen 
cy  shift  of  the  Raman  peak  and  using  the  formulas  with  the  fitted 
parameter  can  give  a  good  correlation  with  the  temperature  measured 
by  using  the  Stokes  to  anti-Stokes  ratio. 

As  an  example  of  using  these  methods,  we  give  now  a  determina¬ 
tion  of  the  temperature  reached  at  the  spot  of  an  Argon  laser 
focussed  onto  thin  layers  of  Silicon  deposited  on  Silicon  substrates 
and  on  Silica  substrates. 

4.3.  Temperature  of  a  thin  layer  of  Silicon  under  continuous 
laser  irradiation. 

Thin  layers  of  Silicon  have  been  deposited  by  C.V.D.  technique 
(Chemical  Vapour  Deposition)  at  a  temperature  of  700°C.  The  layer 
of  Silicon  is  then  a  polycristalline  material.  Two  kind  of  substra¬ 
te  have  been  used  :  Silica  and  monocristal 1 ine  Silicon.  The  diffe¬ 
rence  of  thermal  conductivity  of  these  two  materials  gives  rise  to 
a  different  behaviour  under  laser  irradiation.  The  laser  used  to 
irradiate  the  samples  is  a  continuous  Argon  laser.  The  power  density 
of  the  irradiation  has  been  varied  between  0  and  lOkW/cm  . 
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Figure  9 


Variation  of  the  shift  of  the  zone  center  phonon  of 
Silicon  versus  temperature. 


4.3.1.  Fused  Silica  substrate. 


The  Raman  spectra  for  this  layer  versus  the  irradiation  power 
are  presented  on  Figure  10.  The  spectra  shows  a  dissymetry  towards 
the  low  frequencies  region  whioh  is  due  to  the  presence  of  strains 
inside  the  thin  layer.  This  fact  prevents  the  use  of  the  frequency 
shift  and  the  width  of  the  Raman  peak.  These  strains  are  due  to 
the  difference  between  the  dilatation  coefficients  of  Silica  and 
Silicon  and  this  problem  could  be  neglected  in  the  case  of  an  homo¬ 
geneous  sample. 

The  temperature  of  the  layer  has  thus  just  be  determined  using 
the  curve  of  the  Figure  7. 

4.3.2.  Silicon  substrate. 


The  Raman  spectra  of  a  deposited  layer  of  Silica  on  Silicon 
are  presented  on  Figure  11.  In  this  case,  the  thermal  conductivity 
is  so  high  that  the  temperature  elevation  is  very  small.  The  Figure 
12  shows  the  temperature  for  both  substrates  versus  the  power 
irradiation.  As  can  be  seen,  the  temperature  exceeds  1000°C  for 
the  Silica  substrate  whereas  the  temperature  does  not  reach  100°C 
under  the  same  conditions  for  the  Silicon  substrate. 

5.  RAMAN  SPECTROSCOPY  OF  VERY  HEAVILY  DOPED  SILICON. 

Ionic  implantation  followed  by  laser  annealing  can  produce  new 
samples  doped  with  a  concentration  of  impurities  which  can  exceed 
the  solubility  limit.  This  new  material  could  present  a  technolo¬ 
gical  interest  in  the  field  of  solar  cells  or  I.I.L.  logic  devices. 

The  Raman  spectroscopy  can  characterize  the  doping  effect  by 
controlling  two  different  characteristic  of  the  spectra.  The  first 
one  is  the  occurence  of  new  peaks  in  the  spectra  ;  these  peaks  are 
related  to  the  position  of  the  implanted  impurity.  The  second  one 
is  the  interaction  between  free  carriers  and  normal  or  impurity 
modes.  This  effect  is  responsible  of  a  deformation  of  the  line  shape 
of  these  modes. 
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Raman  spectra  of  a  polycristalline  layer  of  Sil 
deposited  on  Silicon  versus  irradiation  energy. 
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5.1.  The  vibrational  inodes  introduced  by  the  defects. 

When  defects  are  introduced  in  a  perfect  crystalline  lattice, 
the  translation  symetry  is  suppressed  and  the  vibrational  proper¬ 
ties  of  the  crystal  are  modified.  The  calculation  of  the  new  fre¬ 
quencies  is  rather  difficult.  The  simplest  method  to  achieve  it 

f  1  71 

is  the  use  of  the  Green  functions  first  introduced  by  Lifshitz'-  ’ 
and  then  used  by  Dawber  and  Elliott ^ . in  a  model  describing  the 
vibrations  of  a  unique  substitutional  impurity  in  the  Silicon. 

Depending  upon  the  mass  of  the  impurity  two  kind  of  vibrations 
are  distinguishable.  When  this  mass  is  higher  than  that  of  the 
matrix,  the  vibrations  are  of  the  ’’band  mode"  sort.  These  vibra¬ 
tions  can  propagate  along  the  lattice.  When  the  mass  of  the  impu¬ 
rity  is  lower  than  that  of  the  matrix,  the  vibrations  are  of  the 
"localized  mode"  sort  which  are  located  in  the  neighbourhood  of 
the  impurity.  The  Raman  spectroscopy  is  usually  more  sensitive  to 
localized  modes  than  to  band  modes.  In  the  Silicon,  the  only  light 
impurity  with  a  high  solubility  is  the  Boron  so  this  impurity  has 
been  widely  investigated  by  optical  spectroscopy.  The  localized 
mode  induced  by  this  atoms  in  substitutional  positions  have  been 

observed  some  times  ago  in  diffused  or  implanted  materials  by 

VlSI  f ? o  2 1 1 

infrared  spectrometry^  J  or  Raman  spectroscopy'-  ’  . 


5.2.  Free  carriers-phonons  interaction  in  heavily  doped  Silicon , 


The  simultaneous  presence  of  two  kind  of  excitations  a  discrete 
one  and  a  continuous  one  gives  rise  in  the  Raman  spectrum  of  Si 1 i — 
con  to  specific  effect v  .  The  discrete  excitation  corresponds 

to  the  phonon  and  the  continuous  one  to  the  free  carriers.  Both 
excitations  are  Raman  active,  the  phonon  giving  a  sharp  peak  contri¬ 
bution  to  the  spectrum,  the  free  carriers  giving  a  very  broad  conti¬ 
nuum.  The  coupling  between  these  two  excitations  modifies  the  Raman 
spectrum  in  such  a  way  that  what  is  observed  is  not  only  a  simple 
superposition  of  the  sharp  peak  upon  a  broad  continuum  but  an  inter¬ 
ference  which  modify  the  shape  of  the  phonon  line  (Fano  effect). 

This  shape  which  is  Lorentzian  in  Lorentzian  in  absence  of  free 

— ^ - ;)  becomes  assymetric  in  presence  of  elec- 


carriers  (I  (w)  =  — , 

r*  ~ 

lowing 

=  ia(q+e) 


tronic  following  ^  ^0-*  the  equation  : 

2 

where  q  is  proportional  to  the  ratio  of  the 


tensors  for  purely  lattice  and  purely  electronic  Raman  scattering, 
e  =  where  Q  is  the  "dressed"  phonon  frequency  and  F  the  phonon 
damping  in  presence  of  free  carriers  but  without  interference,  .’he 
sign  of  the  factor  q  governs  the  sense  of  the  assymetry  observed  on 
the  spectrum.  This  sign  is  itself  depending  upon  the  type  of  the 
free  carriers  (electrons  or  holes).  The  doping  with  P  type  impuri¬ 
ties  (Boron,  ...)  induces  an  assymet  ry  towards  the  higli  frequency 
whereas  doping  with  n  type  impurities  gives  rise  toassymetry  towards 
the  low  frequency  region. 

5.3.  Experimental  results . 


5.3.1.  Laser  annealing  of  _E  looped  _SH  icon  . 

The  sample  used  for  this  study  is  a  Silicon  (111)  implanted  with 
16  2 

Boron  (10  lons/cm  ,  70  kV) .  Boron  being  a  light  ion,  the  material 
is  not  continuously  amorphized.  After  implantation,  the  Raman  spec¬ 
trum  of  the  sample  is  still  a  cristalline  one.  When  compared  to  the 
Raman  spectrum  of  the  sample  before  implantation,  only  a  decreasing 
in  intensity  and  a  small  assymet i y  towards  low  frequency  are  noticea¬ 
ble  (Fig.  13).  This  assymet  ry  is  due  to  strain:-,  induced  by  the  implan¬ 
tation  process  and  to  small  iTnstrf<  ~  f  amorphous  material.  A  very 
small  peak  appears  at  018  cm  1  which  is  the  frequency  of  the  locali¬ 
zed  mode  due  to  isotope  in  substitutional  position.  After  i  }  ra¬ 

diation  with  a  YAG  laser  (O.hJ/cm*-  at  1.06pm  +  0,2J/cm^  at  0,53pm), 
the  center  of  zone  mode  presents  a  strong  assymet  ry  towards  the  high 
frequencies  and  a  rather  high  intensity  mode  appears  at  618  cm 
This  peak  too,  presents  an  assymetry  towards  the  high  frequencies. 

The  appearance  of  the  618  cm  '  peak  is  a  proof  of  the  high  concentra¬ 
tion  of  substitutional  Boron.  The  assymet ry  of  both  peaks  (normal  modi 
and  localized  mode)  is  the  indication  of  a  high  free  carrier  concen¬ 
tration.  The  comparison  of  laser  annealing  and  thermal  annealing  is 
done  on  Figure  14.  The  thermal  annealing  was  achieved  in  an  oven  at 
a  temperature  of  1000°C  during  2  hours.  As  can  be  seen,  the  Raman 
spectra  are  more  or  less  similar,  the  intensity  of  the  Boron  local 
mode  being  a  little  more  intense  in  the  case  of  the  laser  annealing. 
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5.3.2.  Lase r _ anne a ling_o f_comgens at ed_ Silicon . 

Phosphorus  ions  in  substitutional  site  do  not  give  rise  to  local 
mode  because  of  their  mass  which  is  of  the  order  of  the  Silicon  mass. 
These  ions  should  give  a  band  mode  but  these  band  modes  are  usually 
difficult  to  observe  by  Raman  spectroscopy.  However,  when  a  Silicon 
sample  is  doped  with  Boron  plus  Phosphorus,  the  Phosphorus  atoms  can 
modify  the  local  modes  due  to  Boron  atoms.  As  a  matter  of  fact,  the 
local  mode  of  isolated  Boron  in  substitutional  positions  is  triply 
degenerated.  If  a  Phosphorus  atom  is  situated  in  the  neighbourhood 
of  the  Boron  atom,  the  symetry  is  modified  and  the  triply  degenerated 
mode  at  618  cm  ^  is  split  into  two  modes  at  628  cm  ^  and  600  cm  ^ 

The  Figure  15  presents  the  Raman  spectra  of  Silicon  implanted  with 
2.10^  Boron  per  cm^  under  50  kV  and  2.10^  Phosphorus  per  cm^ .  The 
implantation  energies  have  been  choosen  in  order  to  approximately 
balance  the  penetration  depthes  of  both  impurities.  After  implanta¬ 
tion  a  spectrum  characteristic  of  amorphous  material  is  observed. 
After  laser  annealing  (0,6J/cm^  at  1,06  m  +  0,2J/cm^  at  0,53  m) ,  the 
Raman  mode  of  the  zone  center  is  dissymetric  towards  the  lower  fre¬ 
quencies  region,  proof  of  a  surcompensation  of  free  carriers  as  com¬ 
pared  to  the  sample  implanted  with  Boron  only.  The  concentration  in 
electrons  is  thus  superior  to  that  of  holes  on  the  region  sampled 
by  ‘.he  Raman  laser  beam.  The  local  mode  does  not  see  its  frequency 
modified.  The  intensity  of  this  mode  is  higher  than  for  a  sample 
implanted  with  Boron  only  because  of  the  reduction  of  the  interfe¬ 
rence  due  to  the  partial  compensation  of  free  carriers.  The  fact  that 
the  frequency  of  this  mode  has  not  moved, seems  to  prove  that  there 
are  very  few  Boron  and  Phosphorus  atoms  on  nearest  neighbours  sites. 
On  the  contrary,  in  the  case  of  the  thermal  annealing  of  the  same 
sample  (1000°C  -  2  hours),  the  local  mode  has  moved  to  626  cm  1 
which  is  a  frequency  characteristic  of  a  pair  mode  B-P.  The  other 
peak  of  a  pair  mode  situated  around  600  cm  ^  is  not  observed.  The 
normal  mode  is  slightly  assymetric  towards  the  lower  frequencies. 

The  sample  is  thus  rather  well  compensated  with  just  a  weak  over 
concentration  of  electrons  on  the  path  of  the  Raman  laser  beam. 
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Raman  spectra  of  Boron  and  Phosphorus  implanted  Silicon 
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Laser  annealing  and  Thermal  annealing  have  thus  different  effects 
on  Boron  +  Phosphorus  implanted  Silicon.  In  the  case  of  thermal 
annealing,  atoms  of  Boron  and  Phosphorus  have  enough  time  to  pair 
on  neighbours  lattice  sites  whereas  in  the  case  of  Laser  annealing 
this  is  not  possible. 

5.3.3.  Laser _anne a ling_ of _n-doged_ Silicon . 

The  samples  used  were  cut  along  a  (100)  face  and  implanted  with 
Arsenic  dose  varying  from  5.10^  cm~^  to  5.10^  cm  After  implan¬ 
tation,  the  samples  were  subjected  to  the  irradiation  beam  of  a  Q 
switeched  YAG  laser  equiped  with  a  frequency  doubler.  The  duration 

2 

of  the  pulses  was  typically  100ns  and  energy  density  up  to  2.5J/cm  . 
The  processing  of  the  samples  was  done  at  the  "Groupe  PHASE"  at  the 
C.R.N.  (Strasbourg). 

After  implantations,  the  amorphous  character  of  each  sample  has 
been  verifyied  by  detecting  the  classical  "density  of  states"  Raman 
spectrum  of  Silicon.  This  typical  broad  band  situated  at  480  cm-1 
is  due  to  the  scattering  of  the  light  induced  by  all  the  phonons  in 
the  Brillouin  zone  in  contrary  to  the  sharp  peak  at  520.5  cm-1  obser¬ 
ved  on  crystalline  Silicon  for  which  only  the  mode  of  the  center  of 
the  zone  is  active.  For  the  samples  implanted  in  the  range  5.1015  to 
3.10  As/cm  ,  the  large  amorphous  band  is  well  observed  and  no  thin 
crystalline  peak  can  be  observed  in  addition  to  this  band.  This  signi¬ 
fies  that  these  samples  are  amorphous  on  a  thickness  larger  than  the 
penetration  depth  of  the  laser.  For  the  last  sample,  a  thin  crysta¬ 
lline  peak  at520.5  cm-1  is  detected,  superposed  to  the  broad  amorphous 
band.  This  peak  could  be  due  to  a  part  of  the  sample  which  has  been 
partially  recrystallized  during  the  implantation  process  itself.  The 
high  dose  of  Arsenic  implanted  (5.10  /cm  )  imposed  a  very  high  beam 
current  in  order  to  keep  reasonable  the  time  of  processing.  This 
high  flux  induces  then  an  elevation  of  the  temperature  of  the  sample 
sufficient  to  recrystallize  a  thin  thickness  of  the  just  amorphized 
layer  (Ion  beam  annealing). 

After  laser  annealing,  the  disappearence  of  the  broad  amorphous 
band  confirms  the  recrystallization  of  each  sample.  Two  features 
of  interest  arc  visible  on  the  spectra  presented  on  Figure  16.  A 
peak  at  520.5  cm  '  is  due  to  the  light  scattered  by  the  substrate 
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through  the  perturbed  layer.  Besides  this  line  and  shifted  towards 
the  lower  frequency  region,  a  broader  and  slightly  asymetric  peaks 
can  be  observed.  The  frequency  of  this  peak  varyies  from  515.5  cm-1 
to  511.5  cm"1  as  the  implantation  dose  increases  from  5.1015  to 
5 . 1 01 ^As/cm2 . 

Two  processes  can  explain  the  presence  of  this  shifted  peaks  in 
the  spectra.  First,  a  size  effect  could  be  induced  by  the  presence 
of  microcrystallites  in  the  annealed  area^1^.  This  effect  can  be 
ruled  out,  considering  the  good  monocrystalline  quality  of  the  sam¬ 
ples,  determined  by  Rutherford  backs-scattering  measurements.  This 
other  effect  is  the  coupling  of  the  continuum  of  electronic  transi¬ 
tions  with  the  discrete  phonon  state  (Fano  effect)  discussed  in 
paragraph  5.2.  This  effect  is  usually  weak  in  N  type  semiconductors 
but  the  special  processing  of  these  samples  (implantations  and  laser 
annealing)  can  lead  to  a  free  carriers  concentration  far  above  the 
saturation  and  thus  to  a  very  large  Fano  effect. 

A  ? 

On  Figure  17  are  presented  the  spectra  of  an  1.5  10  As/cr 
implanted  sample  laser  annealed  before  and  after  a  thermal  annealing 
(600°C  -  1  jour).  The  shift  peak  has  disappeared  after  the  thermal 
processing.  This  thermal  processing  has  released  the  out  of  equili¬ 
brium  heavily  doped  state  of  the  sample  thus  leading  to  the  disapa- 
rition  of  shifted  peak  of  the  Raman  spectrum. 

From  the  value  of  the  frequency  shift,  it  is  possible  to  calcu¬ 
late  the  free  carriers  concentration.  The  table  summarizes  the  re¬ 
sults  obtained  for  4  samples  and  compared  to  the  values  of  concen¬ 
tration  measured  by  resistivity  at  Strasbourg. 


TABLE 


- 

0.51016  cm-2 

1  1016  cm" 2 

3  lO^cm-2 

5  10^6  cm-2 

Raman  shift 
(cm-1) 

Raman  con¬ 
centration 

(cm- 3) 

Resistivity 

concentrator 

(cm-3) 

5 

2 1 

1.6  10  1 

2  1 

0.37  10*' 

l 

5.5 

2  1021 

0.52  1021 

8 

5.5  1021 

1  .2  1021 

9 

7.5  1021 

2.25  1021 
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As  it  can  be  seen,  the  concentration  measured  by  the  mean  of 
Raman  scattering  gives  higher  values  than  the  concentration 
obtained  by  electrical  measurements.  This  discrepancy  can  perhaps 
be  explained  by  the  hypothesis  mode  to  compute  the  concentration 
from  the  Raman  shift.  The  model  used  to  achieve  this  computation 
suppose  the  parabolicity  of  the  band  structure*-^.  For  very 
heavily  doped  Silicon,  the  Fermi  level  reaches  regions  where  this 
parabolic  approximation  is  no  longer  valid.  The  flattening  of  the 
bands  increases  the  effective  masse  of  electrons  which  in  turn 
increase* the  joint  density  of  states.  A  complete  calculation  using 
the  real  shape  of  the  band  structure  would  lowers  the  free  carriers 
concentration  and  thus  gives  a  closer  agrement  with  the  values 
obtained  by  electrical  measurements. 

6.  SUB-PICOSECOND  SPECTROSCOPY. 

The  experiments  presented  in  the  previous  chaptei„  .'ere  made 
after  laser  annealing.  This  kind  of  experiments  can  characterize 
the  samples  and  are  thus  of  interest  for  devices  development. 
However,  it  is  not  possible  to  get  information  on  the  process  of 
laser  annealing  itself  from  such  studies.  In  order  to  understand 
the  laser  annealing  process,  time  resolved  spectroscopy  is  required. 
Since  the  real  problem  is  the  understanding  of  the  relaxation 
processes  during  the  first  hundredths  of  femtoseconds,  the  labora¬ 
tory  has  now  built  a  sub-picosecond  laser  source. 

6.1.  The  sub-picosecond  optical  source . 

The  oscillator  part  of  this  source  is  of  the  "colliding  pulse 
modelocked"  (CPM)  type  described  by  C.V.  SHANK  .  A  Rodhamine 
6G  dye  laser  is  continuously  pumped  by  a  C.W.  Argon  laser  in  a  ring 
geometry  cavity.  In  the  path  of  the  cavity,  a  thin  jet  of  saturable 
absorbant  (DODCI)  allows  the  synchronization  and  the  shortening  of 
the  pulses  (Fig. 18).  This  geometry  automatically  involves  that  the 
two  pulses  travelling  in  both  directions  in  the  cavity,  collide  on 
the  absorber  where  the  sum  of  their  intensity  is  sufficient  to  reach 
the  threshold  of  transmission.  Pulse  widthes  of  90  fs  (10”^s)  are 
easily  achieved  in  those  conditions  with  this  Colliding  Pulse  Mode  - 
locking (C. P. M. )  oscillator. 
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This  oscillator  is  now  working  and  is  extensively  tested  in 
order  to  characterize  the  mean  power,  peak  power  and  duration 
of  the  pulses. 

The  second  part  of  the  experimental  set  up  is  the  amplifier. 

The  construction  of  this  part  is  now  in  progress.  This  is  a  four 
stages  amplifier,  each  stage  being  separated  from  the  others  by 
absorber  jets  to  provide  isolation.  The  3  first  gain  cells  are 
pumped  transversally  by  a  YAG  laser.  The  last  stage  is  pumped 
longitudinally  by  the  same  laser.  At  the  output,  the  pulse  is 
recompressed  by  a  Treacy  compressor  which  comprises  two  gratings 
in  a  substractive  mounting.  The  gain  in  peak  power  should  be  of 
the  order  of  10^. 

6.2.  Characterization  of  the  source  oscillator . 

Three  diagnostic  tools  are  used  to  characterize  the  pulses. 

The  first  one  is  an  avalanche  potodiode  followed  by  a  sampling 
oscilloscope.  This  measure  can  only  ensure  that  the  laser  is 
pulsing  as  the  intrinseque  pulse  duration  is  much  smaller  than 
the  response  of  the  apparatus.  The  period  of  recurrence  can  be 
measured  and  is  equal  to  14ns  which  is  the  time  of  flight  of  the 
pulses  in  the  cavity. 

The  spectrum  of  the  pulse  is  qualitatively  checked  by  sending 
a  small  part  of  the  light  diffracted  by  a  grating  onto  a  screen. 

The  lengther  the  trace  on  the  screen  is,  the  shorter  is  the  pulse 
as  time  and  frequency  are  related  through  a  Fourrier  transformation. 

The  third  measure  allows  an  indirect  determination  of  the  pulse 
duration  by  means  of  an  autocorrelator.  This  autocorrelator  is 
made  up  of  a  folded  Michelson  interferometer  where  the  pulse  is 
divided  into  two  parts  on  a  semitransparent  layer.  One  of  the  two 
pulses  is  delayed  through  a  vibrating  prism.  The  two  pulses  are 
then  recombined  on  a  doubler  crystal  of  KDP.  At  the  output  of  the 
crystal  a  beam  of  ultra-violet  light  is  created  when  the  two  pulses 
are  overlapping.  This  light  is  detected  by  a  photomultiplier  and 
the  signal  is  sent  on  the  vertical  channel  of  an  oscilloscope.  A 
voltage  synchronized  with  the  vibration  movement  of  the  prism  is 
~ent  on  the  horizontal  channel  of  the  oscilloscope.  The  net  result 
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is  the  autocorrelation  trace  of  the  pulse  on  the  screen  of  the 
scope.  The  measure  of  very  short  time  is  thus  transformed  in  the 
measure  of  the  elongation  length  of  the  prism  movement.  The 
Figure  19  shows  an  autocorrelation  trace  of  the  oscillator.  The 
width  at  half  height  is  about  8,0fs  (femtosecond). 

The  mean  power  of  the  source  is  easily  measured  with  a  power- 
meter  of  the  calorimeter  type.  It  has  been  determined  to  be 
6mW  for  a  power  of  the  Argon  pump  laser  of  4W. 

Knowing  the  recurrence  period,  the  mean  power  and  the  pulse 
duration,  it  is  possible  to  determine  a  peak  power  of  lOOkW. 

7.  DISCUSSION  AND  CONCLUSIONS. 

Laser  annealing  has  now  a  long  history  and  an  abundant  litera¬ 
ture.  The  interesting  question,  which  still  remains  to  be  answered, 
is  what  are  the  elementary  processes  in  the  interaction  of  a  strong 
radiation  field  with  matter. 

The  process  of  laser  annealing  consists  of  two  sets  of  phenomena. 
The  first  concerns  the  effect  of  a  dense  radiation  field  on  matter 
creating  elementary  excitation  far  from  equilibrium.  The  initial 
hot  plasma  redistributes  through  interactions  between  carriers  and 
ultimately  thermalizes  giving  up  energy  to  heavier  particles  which 
are  the  lattice  constituents.  The  second  set  of  phenomena,  of  a 
completely  different  nature,  concerns  the  modifications  induced  to 
the  solid  as  a  result  of  the  creation  of  the  dense  hot  plasma.  If 
the  laser  pulse  is  short  enough,  the  e-h  density  reached  corresponds 
to  an  amount  of  broken  covalent  bonds  which  is  a  significant  portion 
of  the  total  number  of  bonds,  the  crystal  becomes  fluid  even  at 
T  =  0  K.  Ultimately,  the  solid  sets  in  a  new  phase  whose  structure 
depends  on  the  elementary  mechanisms  of  interaction  and  organization 
of  the  lattice  constituents.  We  shall  discuss  successively  these 
two  states  as  :  i)  direct  laser  effect  and  ii)  consequences  of  the 
laser  action. 


7.1.  Direct  laser  effect. 


The  incident  laser  energy  is  absorbed  by  electron-hole  pairs 
creation  and  by  free  carrier  excitation.  In  indirect  gap  semi¬ 
conductors,  electron-hole  pairs,  are  created  via  indirect  absor¬ 
ption  processes  involving  the  emission  and  absorption  of  phonons. 
Because  phonon  energies  are  much  smaller  than  photon  energy,  the 
amount  of  energy  transferred  to  the  lattice  during  absorption  is 
negligible  in  comparison  to  the  total  amount  absorbed.  The  rise 
of  carrier  density  leads,  in  turn,  to  increased  free  carrier 
absorption.  The  net  result  is  the  production  of  hot-electrons 
and  holes  far  from  equilibrium  which  subsequently  thermalize  with 

the  set  of  the  carriers  and  eventually  with  the  lattice.  The  obser- 
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vations  by  Shank  et  al'-  J  demonstrate  that  with  short,  90  femto¬ 
seconds  optical  pulse,  an  unstable  form  of  highly  excited  state  is 
created  near  the  surface  which  persists  for  a  fraction  of  picosecond. 

f  2  1  -  2  71 

For  very  short  impulsion  of  radiation,  from  the  experiments'-  J 

one  perceives  clear  indications  that  there  are  two  steps  of  the 
laser  action  on  a  crystal  well  separated  in  time.  In  the  first 
step,  the  electromagnetic  energy  transferred  from  the  laser  beam 
to  the  solid  is  retained  in  the  highly  excited  non-equilibrium 
electronic  state.  In  a  second  step,  this  energy  is  transferred  to 
the  lattice. 

Recently,  a  theoretical  model^®)  is  proposed  which  takes  into 

account  the  space  time  evolution  of  the  plasma  during  the  pulse 
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in  order, bo  explain  the  processes  in  densities  higher  than  10  cm 
reached  in  lOOfs  pulses.  The  novelty  in  this  model  is  the  argument 
that  free  carrier  absorption  is  dominated  by  e-h  collisions  with 
a  characteristic  relaxation  time  t  *  3.10~^s. 

A  different  situation  is  reached  with  much  longer  pulses  : 

31  -  3  - 1 

=  10ns  and  a  photon  absorption  rate  g  ^  10  cm  s  :  Auger 
recombination  becomes  the  dominant  recombination  mechanism  at 
these  densities. 

Most  of  the  laser  energy  is  absorbed  by  the  carriers  within  the 
absorption  depth,  liventually,  these  carriers  loose  their  energy 
to  the  lattice,  the  rise  of  the  lattice  temperature  then  depends 
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on  the  distance  they  have  diffused  before  substantial  phonon 
emission  occurs.  At  moderate  density,  the  phonon  scattering 
time  is  t  ^  10~^s.  Screening  does  not  affect  the  rate  of 
intervalley  phonon  emission  until  Ng  'v  10  cm  .  Because 
screening  increases  the  electron-phonon  scattering  time, 
it  not  only  decreases  the  rate  of  phonon  emission  but  also 
enhances  diffusion.  This  increases  the  volume  of  the  region 
in  which  the  energy  of  the  excited  carriers  is  transferred 
to  the  lattice.  Ojwing  to  the  extreme  non  linearity  of  the  hot 
carrier  effects,  it  is  impossible  to  make  an  accurate  estimate 
of  the  precise  temperature  to  which  the  lattice  is  heated  or 
to  determine  the  laser  power  threshold  above  which  melting  will 
occur. 


7.2.  Consequences  of  the  laser  action. 

The  equilibrium  observations  are  clear  :  an  amorphous  or 
glass  solid  is  transformed  into  crystal  under  laser  action  and 
a  crystal  submitted  to  very  short  laser  irradiation  is  transfor¬ 
med  into  amorphous  material.  An  implication  of  both  of  these 
transformations  is  that  melting  preceeds  the  transition.  Another 
alternative  is  that  the  phase  transition  is  directly  induced  in 
the  highly  excited  state. 

The  effect  of  a  dense  plasma  on  the  melting  temperature  is 

itself  an  interesting  problem  of  solid  state  theory.  This  question 

has  been  recently  addressed  by  Bok  and  Combescot  .  It  is  shown 

that  in  the  presence  of  a  dense  plasma,  the  melting  temperature 

of  a  solid  changes.  The  melting  temperature  decreases  with  increa- 

sing  plasma  density.  For  a  laser  pulse  of  IJ/cm  during  10ns,  it 

is  considered  that  e-h  plasira  reaches  a  steady  state  in  a  time 

shorter  than  the  laser  pulse.  Considering  the  plasma  expansion  due 

to  its  high  pressure,  its  collision  with  phonons  and  Auger  recom- 
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bination,  the  highest  plasma  density  is  of  the  order  10  cm 
This  density  is  nevertheless  considered  to  be  sufficient  to  consi¬ 
derably  reduce  the  melting  temperature  so  that  a  metallic  layer 
of  liquid  Silicon  is  formed  at  the  surface. 


The  role  of  a  high  plasma  density  in  laser  annealing  has  been 

discussed  by  Van  Vechten  et  al^  in  a  quite  different  way  but 

still  involving  electron-phonon  coupling  and  lattice  instability 

induced  by  this  interaction.  Above  a  critical  carrier  density 
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estimated  at  8 . 1 0  cm  ,  a  second  order  phase  transition  occurs. 

At  this  plasma  density,  the  bond  charges  will  be  so  depleted  that 
they  will  no  longer  be  able  to  stabilize  the  TA  phonon  modes ^  . 

The  crystal  will  no  longer  resist  shearing  stresses  and  will  beco¬ 
me  fluid.  This  fluid  is  distinct  from  the  normal  molten  phase  of 
Si  the  latter  being  the  result  of  a  strictly  first  order  phase 
transition  driven  by  the  atomic  motion  at  high  temperatures.  The 
assumption  of  Van  Vechten  is  that  the  plasma  is  supposed  to  direc¬ 
tly  induce  the  structural  transformation.  The  energy  is  retained 
in  the  electronic  system  instead  of  being  entirely  associate  with 
the  atomic  motion.  As  the  plasma  becomes  less  dense  due  to  expan¬ 
sion  and  to  transfer  of  energy  to  the  lattice,  the  material  will 
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pass  back  through  the  second  order  phase  transition  at  8.10  cm 
and  covalent  bonding  will  gradually  appear.  The  material  will  fina¬ 
lly  recrystallize  if  this  process  is  relatively  slow  or  will  soli¬ 
dify  in  the  amorphous  phase  if  the  process  is  very  fast. 

This  dense  plasma  phase  could  be  compared  to  the  highly  excited 
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Silicon  which  persists  for  a  fraction  of  a  picosecond''  *  ‘  .  The 

interpretation  of  the  laser  action  differs  nevertheless  with  regard 
to  the  following  step  ;  it  is  generally  considered  that  the  solid 
melts  after  the  initial  interaction  stage. 

In  conclusion,  few  points  appear  clear  today.  The  laser  inter¬ 
action  with  the  solids  results  first  in  the  creation  of  a  highly 
excited  non  equilibrium  phase  which  persists  for  a  fraction  of  a 
picosecond.  The  question  of  how  this  highly  excited  fluid  interac¬ 
tion  with  the  lattice  remains  still  open.  Further  investigation  in 
the  very  short  impulse  regime  are  certainly  desirable  to  clarify 
the  physical  processes  in  laser  annealing. 
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The  electrostatic  energy  per  cell  for  thin  slabs  of  ionic  crystals  is  calculated  as  a  function 
of  the  thickness  of  the  slab,  its  orientation,  and  the  position  of  the  cell  with  respect  to  the 
surfaces.  The  summation  method  used  is  simple  and  quite  general  to  allow  for  direct  appli¬ 
cation  to  slabs  of  any  structure.  The  depolarization  field  associated  with  slabs  of  certain 
orientations  is  also  calculated  and  discussed.  Numerical  results  are  given  for  slabs  of  the 
zinc-blende  and  rocksalt  structures. 


I.  INTRODUCTION 

The  properties  of  small  particles  of  solids 

(10 -  500  A)  and  thin  films  are  sometimes  very 

different  compared  to  those  of  the  bulk  material. 
Since  the  infinite  crystal  approximation  is  widely 
used,  calculation  of  these  properties  has  to  over¬ 
come,  in  most  cases,  the  problem  of  the  presence  of 
surfaces. 

In  particular,  lattice-dynamical  calculations  in 
finite  ionic  crystals  involve  the  evaluation  of  slowly 
and  conditionally  converging  sums,  depending  on 
the  shape  of  the  crystal.  Hence  quantities  such  as 
the  electrostatic  energy  per  cell,  the  internal  field, 
and  the  long-range  interaction  depend  on  both  the 
size  and  the  shape  of  the  finite  crystal. 

De  Wette  and  Schacher1  discussed  these  matters 
in  connection  with  a  planewise  summation  method 
employed  to  calculate  the  internal  field  in  general 
dipole  lattices.  Monkhorst  and  Schwalm2  treated 
the  electrostatic  energy  of  two-dimensional  periodic 
charge  distributions,  where  they  demonstrate  the 
cancellation  of  singularities  in  the  case  of  a  vanish¬ 
ing  dipole  moment  on  the  plane  of  the  film  and  cal¬ 
culated  the  Madelung  energy  of  some  ionic  layers. 

In  this  paper  (paper  I)  we  report  on  the  electro¬ 
static  energy  per  lattice  cell  in  a  thin  ionic  slab,  cal¬ 
culated  as  a  function  of  its  thickness  and  of  the  posi¬ 
tion  of  the  cell  with  respect  to  the  surfaces,  for  slabs 
of  different  orientations  in  the  zinc-blende  and  rock- 
salt  structures.  Surface  cells  show  up  to  10%  higher 
electrostatic  energies  in  certain  orientation,  assum¬ 
ing  they  maintain  the  same  structure  as  the  bulk. 
This  fact  could  be  of  importance  for  surface  recon¬ 
struction  and  other  surface  phenomena.  Moreover, 
the  depolarization  field  associated  with  certain 
orientations  may  also  influence  the  stability  of  the 
surface  layer,  through  the  rearrangement  of  the  elec¬ 
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tronic  charge  it  necessitates. 

In  Sec.  II  we  describe  the  planewise  summation 
method  employed,  a  variant  of  Ewald’s  method, 
which  assures  fast  convergence  and  which  is  very 
simple  and  quite  general  in  order  to  allow  for  direct 
application  in  any  case.  It  has  advantages  also  over 
the  traditional  Ewald’s  method  currently  used  in 
lattice-dynamical  calculations,  for  it  converges  fas¬ 
ter  and  it  provides  results  practically  independent  of 
the  dividing  point  R. 

In  Sec.  Ill  we  calculate  the  depolarization  field  as¬ 
sociated  with  certain  orientations  of  the  slab  as  a 
function  of  its  thickness.  Finally  in  Sec.  IV  we  give 
and  discuss  some  results  for  slabs  of  the  zinc-blende 
and  rocksalt  structures. 

II.  DEFINITIONS  AND  METHOD 

As  defined  in  the  case  of  an  infinite  structure,  the 
Madelung  constant  is  a  dimensionless  number  ex¬ 
pressing  the  electrostatic  energy  per  lattice  cell  with 
respect  to  some  unit  of  length  (usually  the  nearest- 
neighbor  separation  r0).  This  energy  is  of  course  the 
same  for  any  cell  in  an  infinite  structure,  but  in  a 
finite  crystal  it  depends  on  the  position  of  the  cell 
with  respect  to  the  surfaces. 

Hence  it  is  possible  to  define  two  “constants,”  one 
the  average  over  the  whole  finite  crystal  aav(N),  de¬ 
pending  only  on  the  total  number  N  of  cells  in  the 
crystal,  and  one  aU,N),  depending  both  on  the  posi¬ 
tion  of  the  cell  under  consideration  [l  =  Uihh'1] 
with  respect  to  some  reference  point  in  the  crystal 
[/  =  (000)]>  and  the  total  number  N  of  the  cells  in  it. 
In  the  case  of  a  slab,  N  is  the  number  of  the  cell 
layers  and  /  runs  from  ( —  oo ,  —  oc  ,0)  to  ( cc ,  oo  ,N ). 

Between  these  two  constants  the  following  rela¬ 
tion  holds: 
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Jrjad.M,  (1) 

N  l 

where  the  suni  extends  over  the  N  layers  of  the  crys¬ 
tal,  and  the  limit  for  N  -*  oo  should  give  the  value 
of  the  Madelung  constant  for  the  infinite  crystal. 

The  position-dependent  “constant”  can  be  defined 
as 


on  the  consequences  of  using  this  method. 

The  potential  at  a  point  x  due  to  all  ions  in  a 
crystal  slab  is  given  by* 


N  N 

v(x,n)=  2  2  2 


£ 


(4) 


,  N  N  N  N 

a(l,N)={  222  2  ■ 

/*  I /'  =  !  *=1 /r'=l  |*(/  )  X\l,K)  | 

i#r 


(2) 

where  £«  is  the  fraction  of  the  electronic  charge  at¬ 
tributed  to  ion  k. 

In  the  case  of  a  slab  parallel  to  the  plane  (hkl)  of 
the  crystal,  we  can  always  choose  such  a  unit  cell 
that  the  primitive  translation  vectors  ai  and  a2  lie 
on  the  plane  (hkl)  and  a3  lies  out  of  it.  The  thick¬ 
ness  of  the  slab  is  measured  by  the  number  N  of 
cells  superimposed  on  the  direction  perpendicular  to 
the  plane  (hkl).  Since  the  unit  cell  chosen  that  way 
may  be  a  multiple  one,  the  summation  over  k  has  to 
be  performed  over  all  atoms  in  that  cell. 

The  structure  is  considered  as  two  dimensional 
and  a  two-dimensional  reciprocal  lattice  is  associat¬ 
ed  with  it,  defined  by  the  vectors 


where  the  sum  over  l\  and  1'2  runs  to  infinity.  To 
find  a(l2,N )  one  takes  the  limit 


a(h,N)=~  2£  lim 


V(x) 


£ 

|  x(/3,#c)  —  X  i 


(5) 


where  x(/3,k-)  =  x(0,0,/3,k). 

Since  by  assumption  the  vectors  a]  and  a2  lie  on 
the  plane  of  the  slab,  we  can  write 

\W\,l2,ri,K')-x\2 

=  I  x(/i  ,/2  )-f  X||(/3,»f')—  X||  I  2 


+  I  x2(/'3 ,«•')— xx ! 2  ,  (6) 


a2X(a,xS2)  g  (a|Xa2)Xa! 
|a,Xa2|2  ’  2  |  a  i  X  a2 1 2 


We  note  also  that  since  the  structure  is  infinite  in 
two  dimensions,  a(l,N )  becomes  a  function  of  the 
index  /3,  which  gives  the  position  of  the  cell  with 
respect  to  the  surfaces  of  the  slab. 

To  evaluate  the  right-hand  side  of  Eq.  (2)  we  use  a 
variant  of  Ewald’s  method.  We  will  comment  later 
_ I 


where  ||  and  1  designate,  respectively,  the  com¬ 
ponents  of  the  vectors  parallel  and  perpendicular  to 
the  plane  (hkl). 

With  the  use  of  relation  (6)  and  the  integral  repre¬ 
sentation  of  1  /x, 

—  =  —7 —  f  exp(  —  x2p2)dp  ,  (7) 

x  V7T 

Eq.  (4)  becomes 


N  N 


v(x)=-y=r  f  2  2  £expl- l  xii/j.'f'i-xJ  V]  2  ■ 

Vw  ij-  i*'-1  i\.r2 


x(/’l  ,1'2  )+X||(/3,k')-X||  I  1p1]dp  . 


(8) 


The  sum  over  l\l\  in  the  above  expression  is  a  periodic  function  in  two  dimensions  and  can  be  expanded  in 
Fourier  series, 


2  exp[~  |x(/i,/i)  +  X||</3,*')-X|||  2p2J 


= - r  2  np\-ir2\V(hl,h2)\2/p2-2my(h],h1)-[xll(ri,K')-x]l]}  , 

saP  hvh2 


(9) 


where  sa=  |  atxa2|  and  y(hlth2)  is  a  vector  of  the  two-dimensional  reciprocal  lattice,  defined  by  relations 
(3). 

To  keep  in  close  analogy  with  Ewald’s  method,  one  should  split  the  integral  in  Eq.  (8),  into  two  integrals 
over  the  intervals  (0,B )  and  (/?,»)  using  both  sides  of  relation  (9).  This  leads  to  an  expression  containing  er¬ 
ror  functions  which  converges  rapidly,  but  its  numerical  value  depends  on  the  dividing  point  R  more  strongly 
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than  if  unnecessary  splitting  is  avoided. 

Since  we  are  going  to  take  the  limit  of  V(x)  for  x— *  x(/3,k)  [Eq.  (5)],  the  first  exponent  in  Eq.  (8)  will  vanish 
in  the  term  I'3,k'=13k.  It  will  vanish  also  for  1'3  =  l3  and  k'—k"=^k  for  all  ions  k"  lying  on  the  plane  (hkl). 
Hence  we  separate  out  these  terms  (/'3  =/ 3  and  k'=k")  and  we  use  only  the  right-hand  side  of  relation  (9)  in  all 
terms  of  Eq.  (8)  except  the  last  ones,  for  which  the  integral  has  to  be  split  and  use  of  both  sides  of  relation  (9) 
has  to  be  made.  We  get 


N  N  N  N 


V(x,N)=- 


—  2  2  2  2^2  f0  exp[ 

sa  /'  1  / J  =  l  «f'=l  *"=  1  A,. A, 


x1U'3,k')~x1  I  lp -ir 1  y(huh2)  |  /p  ]2 

P 


Xexpf  — 2jriy(A|tA2)*[X||(/j,#f')  — X||]J 


+  -7—  2  £*’■  2  exp[-  ;  xU\,r2)  +  X||(/3,k")  — X||  IV-  I  x1(/3,tc")-jf1|2p2]</p 

V  TT  .  ” 

"  K  I *  I* 


*  /  ,/ 


2£k"  2  L  exp[  —  ;  y(A,,A:)|:/p2-  |  x,(/3,*")-  xL\2p2]^f 

k"  ht.h-%  P 


Xexp)  —  2jny<Ji,,/jj )•[ xi;.(/3,#c")  —  X|i]|  , 


where  k"  runs  over  all  ions  in  the  unit  cell  lying  on  the  same  plane  ihkl)  with  ion  k  (k  included). 
Taking  the  limit  in  Eq.  (5)  using  the  function 


H(x)=  —  erfc(x)  =  —=■  —  f  e  ^dt  , 
r  Vff  x 


the  result3 


1  7  D 

lim  RHiR  |  xU3,k)-x  |  )  — - ! - =  -  — 

lC—lC\l},Kt  !  x(/,,K)-x  i  V  V 

■  > 

and  V(x,N)  from  Eq.  (10),  Eq.  (5)  becomes 
fn  I  Tx/Z  S  .V  V  V 


a(/3,iV)=  —  2$,  — 51  2  2  2  1  sV  2'  1  exp[- :  xI(/3,ic,/3,*c')iy-iT2ry(/i1,A2)i  J/p2] 


V)4 

p 


xexp[  —  2jrty(AlfA2)'X|j(/3,<c,/3,tf')l 


x((f,/',  i  ) 

«"  i\.r2 


+  -jr~2£«"  2  H<>Tr  y^i,h2)\  /R)exp[-2my[huh2)-X\l(K,K’')]-—r$K+T(l3.N) 

^S«  k"  h,.h,  *  17 


(13) 
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where  the  primed  summation  signs  mean  that  terms  with  l\  ,1'2  for  k"=k  or  (A  t A2 )  =  (00)  have  to  be  omitted 
and  T(l3,N)  is  just  the  sum  of  the  ( A  j  A 2 ) = (CX))  terms  given  by 

T(h,N)~^-  - 1 *x<  wu') i +^2&- /0* w 

a  '3  P  sa  P 

The  function  T(l3,N)  diverges  for  p  equal  to  zero.  Its  limit  for  p— ►  0  is 


TU3,N)=-^5-  Vff  2  2  2  2^l*i(WiX)|+|2^ 

s«  ij  «•  «"  «" 


2  2  2  2^+2^|'r  7  -o 


/',  '3  ..  «" 


7  \p-°[P 


due  to  the  neutrality  of  the  unit  cell. 

Evaluating  the  remaining  integral  in  Eq.  (13)  and  using  the  result  (ISa)  we  find 

rn  N  N  N  N  IJT 

aUi,N)m-±  2  2  2  2  2  2  - - — -exp[  —  2ir  |  y(A,,A2)  |  |  xlU3,k,1'3,k')  | 

*  /•  =I  i3=t  «■=!  «"=.i*,.a2  y(AiA2)  | 

'iw, 

— 2iriy(AjA2)*x(/j,#r,/j,if')] 


>  A  ‘  A- 11 


t  • 


I  -  • 


I.  ^JL 


r0vir 


+^~22 U* 2' I >l'i .*”> I >~ -7r  2 E- -  Jr- 2 2 


+  r~22£«£«"  2'  H(ir|y(hI,h2)|/#Oexp[-2ir.7(h„hJ)*Y(,r,ir")] 

«  «"  A,.*, 


2  2^1*1^)  I  - 


K 


where  we  dropped  the  subscript  j  |  as  of  no  importance  since  the  vector  \  ,k2 )  has  no  component  perpendic¬ 
ular  to  the  plane  of  the  slab,  and  where  we  also  reduced  the  first  sum  in  (15a)  since  for  the  pair  of  terms  with 
*=*,,  k'=k\  and  k'=k,  one  has 

I  xLU3,K\,i3,K\ )  |  +  |  Xx(/„»r; ,/'„*,)  |  -2  |  iTl/j./'j )  |  (17) 

for  all  V3  #/ 3,  and  hence  the  sum  over  k  and  k'  is  zero  for  all  i3  ^ l3  due  to  the  neutrality  of  the  unit  cell. 

Equation  (16)  is  our  final  result.  In  this  expression  the  sums  of  //(x)  functions  converge  rapidly  and  so  does 
the  first  sum  provided  that  |  \lU3,k,1'3,k')  |  is  not  too  small,  the  slowest  convergence  occurring  when  l\  =l3f 
For  N— ►  oo  and  l3  well  away  from  the  surface  a(/}, «  )  will  give  the  Madelung  constant  only  if  the  depolariza¬ 
tion  field  is  zero  as  we  will  show  in  the  next  paragraph. 

III.  THE  DEPOLARIZATION  FIELD 

The  function  V(x,N )  defined  by  Eq.  (4)  gives  the  potential  at  a  point  x  of  the  space  due  to  all  ions  in  the 
slab.  Its  values  at  a  point  P  ( lf )  outside  the  slab  is  given  by 


.t  t 


.•  t 
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V[xP,N)=  Vf(xP,N)  +  Ve{xPtN) 


where 


.  I  A  A  *  exp[-2ir|y(h„h2)j  |  xA(/j,*f')-?Px|  ] 

vf(xP,N)=—  I  2^2  - r: - 

5«  *i-*2  |y(«l»«2)| 

Xexpf  -2m7(h,,A2Hx||(/i,ic')-xJ,|,]} 


Vc(xP,N)=-~  2  2  £k-\WiS)-*n\ 
Sa 


as  obtained  by  using  only  the  form  of  the  first  sum 
in  Eq.  (10)  and  taking  the  limit  of  the  first  integral 
inEq.  (14). 

The  first  term,  Vf(xP,N),  in  Eq.  (17a)  represents 
the  local  fluctuations  of  the  potential  and  it  goes  to 
zero  for  a  point  P  very  distant  from  the  slab.  The 
second  term,  Vc(xP,N),  represents  a  constant  poten¬ 
tial  and  it  is  independent  of  the  position  of  the  point 
P  (considered  on  the  same  side  of  the  slab).  This  can 
be  easily  shown  by  taking  the  charge  fraction  §„  of 
a  particular  ion  k,  to  be 

k,— 2k-  (18> 


Then,  assuming  a  positive  xPl  and  a  unit  vector 
Vo=»iX?2/|  ?iXa2|,  Ve(xP,N)  from  Eq.  (17c) 
becomes 


Ve{xP,N)=  K(  +  oo,1V)=1VK(  +  ,I )  , 


where 


K(+fi)=— —  2  M*i(*i>-ii(*'>Ho 
s«  *•#*, 


A  Vc  -  V(  -  oe ,N ) -  V(  +  « ,N) 

=*—  N  2  eA1M)-*iMY*0  f  <22> 

stt  «•*«, 

*',*1 

which  is  locally  modified  by  the  difference  of  the 
corresponding  terms  (17b)  for  local  fluctuations. 
Defining  a  dipole  moment  for  the  unit  cell  by 

Pe  =  2  (23) 

the  potential  difference  AFC  of  Eq.  (22)  becomes 
AFC  =  — ATVvq  .  (24) 

*o 

The  distance  between  the  two  surfaces  of  the  slab  is 

L=(AT  — Dfli+d  ,  (25) 

where  flj.  =  |  a3’ v0 1  is  the  thickness  of  one  unit  cell 
and  d  is  the  distance  of  the  outermost  planes  of  ions 
belonging  to  the  same  cell  (Fig.  1). 

Assuming  that  the  potential  AFf  [Eq.  (24)]  drops 
uniformly  across  the  thickness  of  the  slab,  the  fol¬ 
lowing  field  is  created: 

__  A—  Nd  . 

where  va  =sa-ai  is  the  volume  of  the  unit  cell. 


is  the  potential  created  on  the  one  side  of  a  set  of 
parallel  uniform  charged  planes  through  the  posi¬ 
tions  of  the  ions  in  a  slab  of  one-unit-cell  thickness. 

Considering  now  the  potential  at  a  point  P'  on  the 
other  side  of  the  slab,  one  finds,  apart  from  the  term 
describing  the  local  fluctuations,  a  constant  term 

Ft(xr,iV)=K(-oo,iV)=-F(  +  oo,iV) .  (21) 

Hence  there  is  a  constant  potential  difference  AFC 
between  the  two  surfaces  of  the  slab,  given  by 


FIG.  1.  Schematic  representation  of  a  slab  consisting 
of  N  one-cell  layers. 


28 


LATTICE  DYNAMICS  OF  THIN  IONIC  SLABS.  I.  ... 


3395 


Increasing  the  thickness  of  the  slab  and  taking  the 
limit  of  Eq.  (26)  for  N-+  oo ,  we  find 

E^ooM— (IVv0)v0 .  (27) 

ya 

Hence  the  presence  of  surfaces  parallel  to  charged 
glanes  of  the  crystal  produces  a  macroscopic  field 
E|(  oo )  perpendicular  to  the  surfaces,  and  conse¬ 
quently  the  ions  are  no  longer  under  the  same  condi¬ 
tions  as  in  the  infinite  crystal,  where  no  surfaces  are 
present.3 

To  make  the  slab  equivalent  from  that  point  of 
view,  with  an  infinite  crystal  without  surfaces,  one 
has  to  apply  an  external  field  E„(JV)  equal  in 
strength  and  opposite  in  direction  for  the  total  mac¬ 
roscopic  field  in  the  slab  to  be  zero.  By  doing  that, 
the  energy  per  cell  changes  by  the  amount 

^(M=t2^W),  (28) 

K 

where  Ftx(«,JV)  is  the  potential  at  ion  site  k  due  to 
the  externally  applied  field. 

Using  again  Eqs.  (18)  and  (23)  we  find  for  the 
change  of  the  energy  per  cell 

We{N)= 


lit  NaL 
va  (N  —  DaL+d 


(Pc'V. 


(29) 


which  is  a  negative  quantity  or  zero  if  Pc  lies  on  the 
plane  of  the  slab.  The  factor  depending  on  N  goes 
to  unity  as  N— *  oo ,  and  in  the  limit  we  get 


fVe(oo)=-—(Pc-v0)2. 

va 


(30) 


Hence  the  value  of  the  electrostatic  energy  per  unit 
cell,  obtained  by  Eq.  (16)  for  N— *  oo  and  l}  well  in¬ 
side  the  slab,  is  higher  than  that  given  by  the 
Madelung  constant,  by  We(  oo ), 


aM=  lim  a(N/2,N)  +  r0lVc(ac)  ,  (31) 

<V—*  oo 


where  we  multiply  Wc(  oo )  by  rQ,  the  nearest- 
neighbor  distance  in  order  to  be  consistant  with  the 
other  terms. 


IV.  RESULTS  AND  DISCUSSION 


Using  a  variant  of  the  Ewald’s  summation 
method  we  calculated  the  Coulomb  energy  per  unit 


Rocksalt 


Thickness 
(Number  of  cells) 

FIG.  2.  Electrostatic  energy  per  cell  for  the  central 
layer,  as  a  function  of  the  thickness  of  the  slab,  for  slabs 
of  the  zinc-blende  structure  oriented  parallel  to  (100), 
( llT),  and  (110)  planes. 


Thickness 
(Number  of  cells) 

FIG.  3.  Electrostatic  energy  per  cell  for  the  central 
layer,  as  a  function  of  the  thickness  of  the  slab,  for  slabs 
of  the  rock-salt  structure  oriented  parallel  to  (100),  ( 111"), 
and  (110)  planes. 
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TABLE  I.  Electrostatic  energy  per  cell  for  central  cells  of  thin  slabs. 


Number 

of 

layers  (AT) 

Central 

cell 

(JV/2) 

a(N/2,N)  for  zinc  blende* 
(100)  (110)  (lit) 

a(N/2,N)  for  rocksalt* 

(100)  (110)  (llT) 

1 

1 

0.96795 

1.34448 

1.358  32 

1.615  54 

1.33129 

0.78646 

2 

1 

0.96333 

1.50820 

1.38484 

1.68232 

1.57097 

0.74311 

3 

2 

0.95871 

1.67192 

1.41136 

1.749  10 

1.81064 

0.69975 

4 

2 

0.95829 

1.65347 

1.41135 

1.74832 

1.77588 

0.70006 

5 

3 

0.957  88 

1.63501 

1.41134 

1.747  54 

1.741 12 

0.70037 

6 

3 

0.95788 

1.63670 

1.41134 

1.74755 

1.74470 

0.70037 

7 

4 

0.957  89 

1.638  38 

1.74756 

1.74828 

8 

4 

0.95789 

1.63821 

1.747  56 

1.74788 

9 

5 

1.63803 

1.74749 

10 

5 

1.63805 

1.747  53 

*The  minus  sign  is  omitted. 


cell  of  thin  slabs  of  the  zinc-blende  and  rocksalt 
structures  for  cells  of  the  central  layer  (or  layers 
when  N  is  even).  In  Fig  2  we  show  the  dependence 
of  a(N/2,N)  on  N  for  thin  slabs  of  the  zinc-blende 
structure  oriented  parallel  to  (100),  (111),  and  (110) 
planes.  For  each  orientation,  a{N/2,N)  converges 
very  rapidly  to  a  different  limit.  The  only  limit 
which  coincides  with  the  Madelung  constant 
(flM  =  1.63806)  of  an  infinite  crystal  is  that  for  the 
slab  par;  lei  to  (110)  planes  since  these  planes  are 
through  both  positive  and  negative  ions  and  hence 
neutral.  The  limit  for  a  slab  parallel  to  ( 111)  planes 
is  higher  by  ir/iV^,  which  is  exactly  the  energy  of 
the  unit  cell  given  by  Eq.  (30)  (multiplied  by  r0).  In 
the  case  of  the  slab  parallel  to  (100)  planes  the  limit 
is  higher  by  rrV'i/i  in  accordance  with  Eq.  (30).  In 
both  of  the  latter  cases  the  infinite  planes  forming 
the  slab  are  through  only  one  kind  of  ion  (positive  or 
negative)  and  hence  there  is  a  net  dipole  moment 
perpendicular  to  the  plane  of  the  slab. 

A  common  feature  to  all  three  cases  is  that 
a(N /2,N )  tends  to  its  limiting  value  through  an  ex¬ 
tremum,  a  maximum  in  the  case_of  the  (100)  slab 
and  a  minimum  in  the  cases  of  ( 1  ll)  and  ( 1 10)  slabs, 
in  an  oscillatory  way.  For  the  (110)  slab  the 
minimum  is  most  pronounced,  lying  more  than  2% 
lower  than  the  limiting  value  a*. 

In  Fig.  3  the  dependence  of  a(N/2,N)  on  N  is 
shown  for  thin  slabs  of  Jhe  rocksalt  structure  orient¬ 
ed  parallel  to  (100),  (111),  and  (110)  planes.  In  both 
(100)  and  (110)  cases,  a(N /2,N)  tends  to  the  value 
of  the  Madelung  constant  (1.74756)  since  both 
these  planes  are  neutral,  through  a  minimum  which 
is  almost  4%  lower  than  the  limiting  value  for  the 
(110)  slab.  For  slabs  parallel  to  ( 1  ll)  planes  the  lim¬ 
it  is  higher  by  ir/3  in  accordance  with  Eq.  (30). 

In  Table  I  we  summarize  the  numerical  results  for 
the  six  slabs  of  different  orientations  and  structures 


previously  described.  The  most  rapid  convergence 
of  a(N/2,N)  to  its  limiting  value  is  observed  for 
slabs  of  ( 1  ll)  orientation  in  both  structures,  and  the 
slowest  is  for  slabs  oriented  parallel  to  (1 10)  planes. 

In  Table  II  we  give  the  values  of  the  electrostatic 
energy  per  cell  a(l3,N)  for  cells  at  different  dis¬ 
tances  from  the  surface  for  a  thin  slab  of  the  zinc- 
blende  structure,  oriented  parallel  to  (100),  (110), 
and  (111)  planes.  It  can  be  seen  that  the  limiting 
value  is  reached  within  even  the  second  layer  for  the 
(111)  slab,  while  for  the  (110)  slab  it  is  reached  at  a 
deeper  layer. 

In  Fig.  4  we  give  the  depedence  of  the  electrostat¬ 
ic  energy  per  cell  [r0W''c(JV),  Eq.  (29)]  due  to  a  mac¬ 
roscopic  field  equal  in  strength  and  opposite  in 
direction  to  the  field  created  by  the  slab  of  N  layers. 
Among  the  cases  considered  here,  such  a  field  is 
created  by  the  (111)  and  (100)  _slabs  of  the  zinc- 
blende  structure  and  by  the  (111)  slab  of  rocksalt 
structure.  The  corresponding  limiting  values  are 
reached  very  slowly  but  the  convergence  is  the  same 


TABLE  II.  Electrostatic  energy  per  cell,  in  a  slab  of  20 
layers,  for  layers  in  different  distances  from  the  surface. 


Layer 

from 

surface 

U,) 

a(/3, 20)  for  zinc  blende' 
(100)  (110) 

1 

(llT) 

0 

0.96292 

1.49127 

1.38482 

l 

0.95830 

1.65499 

1.41134 

2 

0.95789 

1.63654 

1.41134 

3 

0.95789 

1.63822 

1.41134 

4 

1.63805 

5 

1.63807 

6 

1.63806 

The  minus  sign  is  omitted. 
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Thick  n«st  (numbtr  of  colls) 

FIG.  4.  Electrostatic  energy  of  a  cell  in  the  macroscop- 
ic  field  of  the  slab  as  a  function  of  its  thickness,  for  ( 1  lT) 
and  (100)  slabs  of  the  zinc-blende  structure  (left-hand 
scale)  and  (111)  slab  of  the  rocksalt  structure  (right-hand 
scale). 


for  the  (100)  zinc-blende  and  the  (111)  rocksalt 
slabs.  For  slabs  with  50—70  layers  the  calculated 
values  are  within  1%  of  their  limits. 

It  is  worthwhile  to  note  that  although  for  N— *  oc 
Eq.  (2)  is  the  definition  of  the  Madetung  constant, 
its  limit  for  W— >■  oo  depends  on  the  orientation  of 
the  slab.  This  is  due  to  the  fact  that  the  series  in  the 
right-hand  side  of  Eq.  (2)  is  not  absolutely  conver¬ 
gent  and  hence  its  limit  depends  on  the  way  of  sum¬ 
mation.  Any  method  summing  first  on  infinite 
planes  cannot  avoid  creating  a  constant  potential 
outside  the  planes  and  hence  a  depolarization  field  if 
the  dipole  moment  of  the  cell  has  a  nonzero  com¬ 
ponent  perpendicular  to  the  plane  of  summation. 

Monkhorst  and  Schwalnr  calculated  the  electro¬ 
static  energy  of  a  layer  of  the  NaCl  structure  paral¬ 
lel  to  the  (100)  plane  in  different  distorted  configu¬ 
rations.  Our  result  for  one  layer  of  the  undistorted 


lattice  is  in  exact  agreement  with  their  correspond¬ 
ing  result  (namely  the  value  —1.615  54),  The  in¬ 
crease  of  the  energy,  as  their  layer  is  separated  into 
two  oppositely  charged  sheets,  is  due  to  the  same 
reason  we  already  explained  and  the  term  E}  in  their 
energy  of  the  layer  is  equivalent  to  our  result  of  Eq. 
(28). 

A  final  remark  concerns  their  suggestion  to  treat 
the  electrostatic  energy  of  systems  of  point  charges 
in  terms  of  coefficients  <(tKK>  of  the  products. 
For  diatomic  ionic  crystals  there  is  only  one  such 
coefficient,  the  Madelung  constant  itself  (since 
£2=—  £i)-  For  crystals  with  more  atoms  per  unit 
cell,  with  at  least  three  different  such  an  expres¬ 
sion  will  save  considerable  computation  time  and 
will  point  out  relations  between  structure  factors  in 
structures  and  superstructures.  We  have  found  that 
the  Madelung  constant  of  the  chalcopyrite  structure 
(which  is  a  superstructure  of  the  zinc  blende)  can  be 
expressed  as 

a ch  =  —  2. 1 39  25( £1 + £j) — 2. 273  78£j£j  ,  (32) 

where  £1  and  are  the  two  cation  charge  fractions. 
For  £|=£2=  1,  ach  equals  4  times  the  Madelung 
constant  of  the  zinc-blende  structure,  since  the  unit 
cell  for  the  ideal  chalcopyrite  is  4  times  larger. 

After  this  paper  was  submitted  for  publication,  V. 
V.  Avilov4  proposed  formulas  to  calculate  the  elec¬ 
trostatic  energy  of  planar  lattices.  Apart  from  an 
obvious  typing  error,  the  potential  of  a  planar  lattice 
of  charges  at  a  distance  |  z  |  from  its  plane, 
represented  as  a  Fourier  series  in  the  plane,  is  given 
in  this  paper  as  exactly  the  same  as  in  our  calcula¬ 
tions.  In  our  opinion,  the  uniform  background 
charge  distribution  of  opposite  sign  assumed  to  re¬ 
move  the  divergence  of  the  sum  results  in  an  unreal¬ 
istic  picture  for  a  slab  or  an  infinite  crystal  with  two 
parallel  surfaces.  It  does  in  fact  implicitly  compen¬ 
sate  for  the  depolarization  field  through  a  multilayer 
(or  a  two-layer)  sandwich  of  uniform  charge  distri¬ 
butions  of  opposite  polarity,  but  without  its  source 
and  importance  being  demonstrated. 
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Please  make  the  following  changes  to  our  paper. 

(1)  Page  3391,  in  Eq.  (2),  instead  of 

N  N  N  N 

7  2  2  2  2  . 

/'-i  — >  «'-i 


t  2  2  2  ■ 

*  /'  «' 

(2)  Page  3391,  in  Eqs.  (4)  and  (8),  the  sum  over  *'  runs  from  1  to  the  number  of  ions  in  the  unit  cell  (not  to  N ). 

(3)  Page  3392,  in  the  first  line  of  Eqs.  (10)  and  (13),  instead  of 

-  r~  N  N  N  N 

—  2  2  2  2  ■ 

Sa  /'-i  'j-1  ,'.i  «",i 

-'3 


2Vrr  "  -  .  , 

—  2  XfK  ...  . 


(4)  Page  3393,  in  the  first  line  of  the  text,  instead  of . . .  terms  with  l{,li  for  . . . ,  write  . . .  terms  with  )  -  (0, 0)  for 

(5)  Page  3393,  in  Eqs.  (14),  (15a),  and  (15b),  instead  of 


2  2  2  X(* 


2  2  «*'•••  • 

(6)  Page  3393,  in  the  first  line  of  Eq.  (16),  instead  of 

N  N  S  N 

2  2  2  2  2  2' 

*  (j-i  'j-1 1  «"-i  *j*2 

/'#i3  *!-*" 


2  2  2  2 

*  - 1  «'  »,.Hj 


29  6372 


PHYSICAL  REVIEW  B 


VOLUME  28,  NUMBER  6 


15  SEPTEMBER  1983 


Lattice  dynamics  of  thin  ionic  slabs.  II.  The  long-range  forces 

G.  Kanellis*  J.  F.  Morhange,  and  M.  Bafkanski 
Laboratoire  de  Physique  des  Solides  de  I'Universite  Pierre  et  Marie  Curie  associe  au  Central  National 
de  la  Recherche  Scientifique,  4  Place  Jussieu,  F-75230  Paris  Cedex  05,  France 
(Received  15  March  1982;  revised  manuscript  received  1 1  March  1983) 

The  Coulomb  interaction  is  calculated  for  thin  ionic  slabs  of  any  structure  and  orienta¬ 
tion.  For  ions  lying  on  different  planes,  parallel  to  the  plane  of  the  slab,  the  direct-space 
double  sum  is  transformed  into  a  rapidly  converging  double  sum  in  reciprocal  space  only, 
while  for  ions  lying  on  the  same  piane,  a  two-dimensional  version  of  Ewald’s  method  is 
used.  Finally  we  comment  on  the  macroscopic  field  associated  with  long-wavelength  vibra¬ 
tions  of  an  infinite  lattice  on  the  basis  of  a  two-dimensional  summation  method 


I.  INTRODUCTION 

i  iieu  is  great  interest  in  the  physics  of  surfaces 
and  thin- film  phenomena  and  much  attention  has 
bee; i  d/av.ii  to  the  study  of  surface-localized  vibra¬ 
tional  modes  of  crystals.  To  investigate  these  modes 
in  !m>  .-.labs  of  ionic  crystals  by  direct  lattice- 
dv  au  ical  calculations,  one  needs  to  calculate  the 
h  -ug-winge  interaction  between  plane  lattices.  This 
task  involves  the  evaluation  of  dipole  lattice  sums 
wilt  an  infinite  number  of  terms,  in  two  dimensions. 
Since  h'ttice  sums  in  general  are  conditionally  con¬ 
verged,  two-dimensional  summation  methods  in 
tl  ems.  Ives  imply  conditions  leading  to  different  re- 
which  have  to  be  interpreted. 

Niiboer  and  de  Wette1  and  de  Wette  and  Schach- 
e.  discussed  these  matters  when  calculating  the 
ii  i-M\d  field  in  dipole  lattices  and  they  pointed  out 
hev.  the  sums  are  influenced  by  the  orientation  of 
the  .’ianes  on  which  the  summation  is  performed 
wuh  respect  to  the  dipole  direction.  They  propose  a 
pianewi.se  summation  method  which  can  be  used  to 
calculate  the  long-range  interaction  in  thin  slabs,  but 
thev  do  not  derive  expressions  for  such  a  case. 

Detailed  calculations  on  finite  ionic  slabs  have 
been  arried  out  by  Tong  and  Maradudin.1  These 
authors  treated  in  detail  the  long-range  interaction 
■I.  the  case  of  a  slab  of  the  rocksalt  structure  parallel 
t-  fie  (100*  plane,  using  a  method  proposed  by 
Mackenzie.  The  same  case  has  also  been  treated  by 
Biyksiii  a.-.  I  Firsov4  and  by  Jones  and  Fuchs5  for  re¬ 
tarded  and  unretarded  interaction,  but  on  different 
hjsev 

In  this  i  a  per  we  derive  explicit  expressions  for  the 
coupling  coefficients  due  to  unretarded  long-range 
foi"i>  bt  i ween  the  plane  lattice4  ir  an  ionic  slab, 
whatever  its  structure  and  orientation.  This  is  done 
folk- a  mg  a  procedure  similar  to  the  one  presented  in 
Me  |  my  thing  paper,  referred  to  as  paper  !,  where  we 


calculated  the  electrostatic  energy  per  cell  in  thin 
ionic  slabs.  We  believe  the  method  we  propose  here, 
a  variant  of  Ewald’s  method,  is  quite  general, 
straightforward,  and  simple,  and  hence  more  con¬ 
venient  to  use. 

Finally  we  comment  on  the  macroscopic  field  as¬ 
sociated  with  long- wavelength  vibrations  of  an  in¬ 
finite  lattice.  In  fact  a  planewise  (two-dimensional) 
summation  method  for  the  Coulomb  interactions 
can  never  give  the  formal  results  of  Ewald’s  method 
for  the  macroscopic  field  without  additional  as¬ 
sumptions  concerning  the  effect  due  to  the  presence 
of  surfaces. 


II.  COULOMB  INTERACTIONS 
BETWEEN  IONS 


We  consider,  as  in  paper  I,  a  crystal  slab  consist¬ 
ing  of  N  one-cell  layers,  parallel  to  the  ihkl)  plane  of 
the  crystal  structure.  By  appropriate  choice  of  the 
unit  cell,  the  primitive  translation  vectors  at  and  a2 
can  be  taken  to  lie  on  the  (hkl)  plane,  while  a,  lies 
out  of  it.  The  slab  is  considered  as  a  two- 
dimensional  periodic  structure  and  a  two- 
dimensional  reciprocal  lattice  is^  associated  with  it. 
defined  by  the  vectors  b(  and  b2  [Eq.  (3)  of  paper 

I] 

The  field  at  a  point  x  due  to  all  point  dipoles 
P  (/',*')  at  lattice  sites  x(/',*'),  /'  —  (/',  is 

given  by 


,•  K.  e  oxa6xp 


v  exp[2iriy •  x(/',k')} 
g  ,•  |  xU\k')-  x  | 


where  a,P—  (1,2,3)  number  the  axis  of  some  Carte¬ 
sian  coordinate  system  Ox^xiXit  x'  runs  over  all 
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ions  in  the  unit  cell,  /'  runs  over  all  unit  cells  from 
(-oo,-oo,0)  to  ( oo , oo , JV ),  and  y  is  a  two- 
dimensional  wave  vector. 

We  note  that  the  index  /3  labels  the  one-unit-cell 
layers  of  the  slab.  Hence  the  index  (13,k)  labels  the 
atoms  in  the  long  unit  cell  of  the  two-dimensional 


periodic  structure  defined  by  the  vectors 
(a|,a2,JVaj).  To  evaluate  the  sum  over  l\  ,l\  on  the 
right-hand  side  of  Eq.  (1)  we  proceed  the  same  way 
as  in  evaluating  the  double  sums  in  paper  I. 

Using  the  integral  representation  of  1/x  [Eq.  (7) 
of  paper  I]  and  again  setting 


I  x  1 2=  I  x(/', )  +  X||(/3 X,,  |2+  I  Xil  \l\K')- X,  1 2  ,  (2) 

where  ||  and  1  designate  the  components  of  the  vectors  parallel  and  perpendicular  to  the  plane  (hkl),  respec¬ 
tively,  relation  (1)  becomes 

ejx)=  2  /<>*-  V+2ir,y'?] 


X  — 7=-  2  «Pi  -  I  x(i|,/'2)-»-X||(/3,#r')-X|||  V  (3) 

+  2ir/y-[x(/i  ,/2)+  Xj|(/3,K')-3f||]j</p  , 

since  the  vector  y  lies  on  the  plane  of  the  slab. 

The  sum  over  l\  ,/'2  on  the  right-hand  side  of  Eq.  (3)  is  a  periodic  function  in  two  dimensions  and  can  be  ex¬ 
panded  in  Fourier  series 

-7=-  2  exPl  -  I  x(/’i,/2)  +  X||^3»'f'>-*||l  V+2iri7-[x(/',,/i)-l-x,|(/3,ic')-X||]| 

n‘\A 

2  exPl  I  yMi,A2)+y  1 2/p2— 2^/y(Ai,A2)-[xM(/3,#c')— x,|]}  ,  (4) 

saP  Aj 


where  sa  =  |a,Xa2|  and  y(huh2)  is  a  vector  of  the 
two-dimensional  reciprocal  lattice. 

Here  again,  as  explained  in  paper  I,  one  can  intro¬ 
duce  a  dividing  point  R  on  the  p  axis  and  evaluate 
the  integral  in  Eq.  (3)  separately  over  the  interval 
(0,R)  using  the  right-hand  side  of  relation  (4)  and 
over  the  interval  (R,oo)  keeping  the  integrand  as  it 
is. 

Splitting  the  integral  in  Eq.  (3)  is  necessary  only  if 
;  Xi(/3,k')  - Xi  |  =0  for  some  i.e.,  if  the  vec¬ 

tor  x(1'3,k')~x  lies  entirely  on  the  plane  ( hkl) 
through  the  ion  (/3,»c'). 

In  calculating  the  Coulomb  interaction  between 
ions,  one  is  interested  in  the  field  at  ion  sites  ?(/ 3.x). 
It  will  happen  then  to  be  |  x2(13,k')  —  Xi(litK)  |  =0 
for  =  litK",  where  k"  labels  all  ions  of  the  unit 
cel!  on  the  same  (hkl)  plane  through  ion  k  (k  includ- 


Hence  we  can  consider  the  field  Ea(  x )  as  being 
the  sum  of  two  components, 

Ea(x)=E{"(x  )+£^2’(x)  ,  (5) 

where  E{J](x)  is  the  field  at  point  x  due  to  all  ions 
outside  the  plane  ( hkl)  through  that  point  and 
£q2i(x  )  is  the  field  at  the  same  point  due  to  all  ions 
lying  on  the  plane  (hkl).  We  distinguish  the  follow¬ 
ing  two  case':: 

Case  A:  |  xl(!',,k')  —  x2  |  =£0.  Since  we  are  later 
going  to  replace  x  by  x(13,k),  this  will  be  the  case 
for  U'3,k')^(Ij,k")  as  already  indicated.  We 
proceed  by  evaluating  the  integral  in  Eq.  (3),  using 
only  the  right-hand  side  of  relation  (4).  Then 
for  £„"(?  )  we  have 


eJ *(x  )=  2  2P0(/ - ~exp(2rriy  x  ) 

e  dx°dxe 


x  2  “ — ,  ,  exPl  -2tt|  y(huh2)+y  |  |  x[(l'i,K')-xi 
A, ,*2  !  y(A,,h2)+y  | 

-2jriy(h,,h2)-[X||(/3,»c')-X||]|  . 


-•  .  • 


•  .  • 


t-  « 


(6) 
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Before  carrying  out  the  differentiation  in  Eq.  (6),  we  note  that  the  unit  vector  v0  perpendicular  to  the  plane 
(hkD  is,  by  the  assumption  for  the  choice  of  the  unit  cell, 

aiXa2 

Vq  — - =cosa  x0i+cosj8x02+cosyx03  ,  (7) 

|a,Xa2| 

with  respect  to  the  coordinate  system  Ox\X2x 3. 

Hence 

0  0 

- —  |  xlU'},k')-x i  I  =sg- — [x(/'3)»c')-vn-x-v0]=  -s?cosa  ,  (8a) 


where 


|  +  1  if  xU'3,k')-v0>  x-v0 
1  —  1  if  x(/3,/f•,)•  v0  <  x-v0 


i=°. 

Carrying  out  the  differentiation  and  evaluating  the  derivatives  for  x  —  x(/3,» c),  we  get 
E'aU j*)=  2  'Sj,p^3,K')Qap{li,K;l’J,K,\y)exp[2iriS  x(I3,K)]  , 

l\ 


n  n  . i'  <i->  ^  [ya(hl,h2)+ya][yg(hl,h2)+yfl] 

Qap{h>K’h’K  !  y)  =  — -  2,  — — - - — — - - - -+  |  y(A,,A2)  +  y  |cosacos# 

s«  *,.*2  |y(A,FA2)+y| 

+  i  Ly  a  <  h  1  ,A  2 )  +ya  }cos/3 +  [y  fi(  h  i ,  h  2 )  +  y  p  ]cosa  { 


Xexp{  —  27r  |  y(huh2)+y  |  |  xx(/3,# c,/3,*')  | — 2tt* y ( A , , A 2 )- x(/3,#c,/^ ,«-' )]  . 


In  the  last  equation  we  dropped  the  index  ||  as  of  no  importance,  since  the  vector  y(A,,A2)  lies  on  the  plane  of 
the  slab. 

it  is  worthwhile  noticing  that  the  above  expression  for  Qag(h,>c;r},K'  |  y)  is  a  regular  function  of  |  y  |  hav¬ 
ing  a  well-defined  limit  for  |  y  |  — *0.  This  can  be  easily  verified  by  separating  the  ( A  t  A 2 )  ==  (00)  terms  on  the 
right-hand  side  of  Eq.  (10).  Their  sum  tends  always  to  zero  for  y  — >0  regardless  of  the  direction  of  the  vector 

y- 

Case  B.  |  Xi(/'3,K')—  xx  |  =0.  This  will  be  the  case  for  x  =  xU3,k)  and  (/3,/r’)  =  (/3,*")  where  k"  runs  over 
all  ions  of  the  unit  cell  lying  on  the  same  (hkl)  plane  with  ion  k.  In  this  case  we  split  the  integral  in  Eq.  (3) 
into  two  integrals  over  the  intervals  (0,B)  and  (R,oo)  using  both  sides  of  relation  (4).  We  get  for  E„'Cx  ), 

*2  //(/?  |  x(/j,/i,/3,#c'  )_x  |  )exp[2«'ry-x(/'1,/i,r3,#c">J 

k"  0  OXaOXp  J,  J, 

+  2  Hirr\yifi„h2)  +  y\/R) 

"  *0*2 

Xexpj27ri[y(A„A2)+y]-x- 2iriy(A,,A2)-x(/3,K")|  ,  (11) 


where  we  used  the  function 
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H(x)=~~  f"e-,2dt 
Vn  x 


(12) 


to  express  both  integrals. 

The  function  H(x)  diverges  for  x=0.  Hence  the  term  U\l'2)=(00)  on  the  right-hand  side  of  Eq.  (11)  will 
diverge  for  k"=k.  Performing  the  differentiation  and  replacing  x  by  x(/3,x)  we  get  for  Of ), 


fL2,(s>=2 

p 


2  Pp(h,K")Qafi(h,K\h>*"  I  y)+Qafi(l}>K'>l)yK  I  7 ) 

y,K 


exp[2ir/y-x(/3,<f)]  , 


(13) 


where 

Qap(h>K\li>K"  \  y)=^32  Ht tfilR  I  x(k,1\,12,k")\  )exp[2ir/y-x(x,/',,/i,K")] 
l\A 

-1T~  2  \ya^uh2)+ya\\y»^\M)-\-yp\H(iT\ y(hx,h2)+y  I  /R) 
Ks“  Mi 


Xexp[2iriy(h1,A2),x(«f",x)]  , 


(14) 


QapUiyK'JiyK  I  7>  is  the  diverged  term,  to  be  evaluated,  and 


Haf>(  I  *  |  )  = 


dxadxp 


H  (  |  x  |  )  . 


(15) 


To  find  Qa0(h>K>h>K  I  7)  one  has  to  subtract  from  the  right-hand  side  of  Eq.  (11)  the  contribution  in  the  po¬ 
tential  at  ion  site  (/3,x)  due  to  the  ion  (/3,x)  itself. 

This  is  done  by  taking  the  limit 


lim 


dxadxfi 


MR  |  x(0,0,/3,*)-x  | )_ - - - 

|  x(/3,x)— X  I 


=  lim 

V  -►Ti/j.k) 


dxadxp 


R|x(/3,k)-x|  V9 


2  |T(/j,«)-"7| 

—  f  exp(  —t2)dt 

IT 


3  Vff 


bag  ,  (16) 


where  b„p  is  Kronecker’s  6.  By  naming  the  function 
in  the  square  brackets  of  Eq.  (16)  as  H^p(x)  and 
considering  the  above-found  result  as  its  limit  for 
x— *0,  we  can  retain  Eq.  (14)  as  the  general  formal 
expression  for  Qap(l3,K\l3,K"  \y)  bearing  in  mind 
that  in  the  case  of  k=k",  one  has  to  replace  in  the 
term  (/i/2)=(00)  the  function  H^x)  by  H^p(x). 

We  note  again  that  the  right-hand  side  of  Eq.  (14) 
is  a  regular  function  of  1 7 1 .  The  sum  of  the 
A  i  Ar2  =  (00)  terms  goes  to  zero  for  |  y  |  — »0,  what¬ 
ever  the  direction  of  the  vector  7  is. 

For  the  coefficients  Qap(l3’K>l3>K"  I  71  given  by 
Eqs.  (10)  and  (14),  the  following  symmetry  relations 
hold: 

Qath,K\l'iS  I  7>=G*(W3.*'  I  7  >  >  d7a) 

QapU3,xA,K'  I  y)=Q,ap(l3,K’;l3tK  |  7> ,  U7b) 
Qae(h,KU3,K'\?)=QUh,*-,li>K'\ -7)  •  (170 


Relation  (17a)  can  be  verified  by  inspection  from 
Eq.  (10)  and  by  calculating  H„p(x)  for  Eq.  (14), 
Hapix)  being  also  symmetric  in  xa,xp.  Relation 
(17b)  is  also  true  for  Eq.  (10)  because  of  the  condi¬ 
tion  (8b)  while  for  Eq.  (14)  the  summation  indices 
l\  ,l2  have  to  be  replaced  by  — 1\ ,  — 12.  Finally,  re¬ 
lation  (17c)  is  verified  for  both  Eqs.  (10)  and  (14)  by 
replacing  the  indices  huh2  by  —h\,—h2. 

III.  ON  THE  MACROSCOPIC  HELD 

In  the  three-dimensional  case  the  macroscopic 
field  is  produced  by  the  infinite  charged  planes, 
through  the  positions  of  the  vibrating  ions,  perpen¬ 
dicular  to  the  direction  of  the  wave  vector  q,  for 
long-wavelength  vibrations  (q— »0).  Its  value  turns 
out  to  be  equal  to  a  nonregular,  for  the  q— >0  term, 
of  the  sum  of  the  (hih2h3)=(000)  terms  in  the  ex¬ 
pression  for  the  exciting  field  at  a  dipole  site.6  It  is 
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given  by 


gmic=_i£_J - (18) 

|ql  |ql  * 

where  va  is  the  volume  of  the  imit  cell  and  p(*)  is 
the  dipole  moment  associated  with  the  ion  k. 

In  the  preceding  paragraph  we  have  shown  that 
such  a  term  exists  neither  in  Eq.  (8)  nor  in  Eq.  (11) 
for  the  exciting  field.  Since  the  macroscopic  field, 
according  to  Eq.  (18),  is  on  the  same  direction  with 
vector  q,  we  have  to  conclude  that  there  is  not  a 
macroscopic  field  lying  on  the  plane  of  the  slab. 
This  is  of  course  what  one  should  expect,  since  in 
the  present  case  there  exist  no  infinite  oscillating 
charged  planes  perpendicular  to  the  plane  of  the 
slab,  but  rather,  plane-charged  “ strips ”  infinite  in 
one  dimension.  Consequently  in  slab-shaped  crys¬ 
tals,  longitudinal  and  transverse  phonons  propagated 
on  the  plane  of  the  slab  must  be  of  the  same  fre¬ 
quency,  as  far  as  no  displacements  perpendicular  to 
the  plane  of  the  slab  are  involved. 

On  the  other  hand,  there  are  infinite  oscillating 
charged  planes  parallel  to  the  plane  of  the  slab. 
Hence  a  macroscopic  field  must  appear  perpendicu¬ 
lar  to  that  plane  due  to  vibrations  along  this  direc¬ 
tion.  To  calculate  that  field  we  consider  the  poten¬ 
tial  distribution  inside  the  slab  (Fig.  1).  We  have 
shown  in  paper  I  [Eq.  (22)]  that,  depending  on  the 
orientation,  a  constant  potential  difference  AFc(A0 


HG.  1.  Potential  distribution  in  a  thin  slab.  (1)  Static 
potential  created  by  the  charged  planes,  (2)  external  poten¬ 
tial  compensating  the  former  one,  (3)  final  distribution  re¬ 
sulting  from  the  superposition  of  (1)  and  (2),  (4)  external 
potential  when  the  slab  is  considered  as  part  of  an  infinite 
structure,  (5)  potential  distribution  in  an  infinite  crystal, 
i.e.,  superposition  of  (1)  and  (4). 


is  created  between  its  two  free  surfaces. 

We  would  like  to  emphasize  at  this  point  that  the 
interaction  due  to  electrostatic  potentials  calculated 
in  a  preceding  paragraph  is  correct  and  accurate  as 
far  as  the  slab  is  considered  to  be  taken  out  of  an  in¬ 
finite  structure.  Moreover,  no  question  arises  for 
the  zero-wave-vector  limit  perpendicular  to  the 
plane  of  the  slab,  since  no  periodic  boundary  condi¬ 
tion  has  been  imposed  along  that  direction.  What 
we  are  going  to  show  is  that  the  well-known  macro¬ 
scopic  field  in  the  infinite-crystal  case  is  due  to  in¬ 
finite  distant  planes  (missing  in  the  case  of  a  slab) 
and  moreover  how  this  field  could  be  taken  into  ac¬ 
count  especially  in  the  case  where  the  surfaces  of  the 
slab  are  not  free. 

To  derive  an  appropriate  expression  for  the  con¬ 
stant  potential  difference  AVC  between  the  two  sur¬ 
faces  of  the  slab,  we  recall  Eq.  (17c)  of  paper  1  for 
the  constant  term  Vc  of  the  potential  at  a  point 
P{  Up )  outside  the  slab, 

Vc(xp,N)=~—  X  X  &-l*i(/3.*'>-*,il  , 

Sfl  = 


(19a) 

where  is  the  fraction  of  the  electronic  charge  at¬ 
tributed  to  ion  k'.  At  a  point  P'i\p ■)  on  the  other 
side  of  the  slab,  the  same  term  of  the  potential  will 
be 

vc{xp.,n)=~—  x  x^rw^'>-vti  • 

/'  =1*^=1 


(19b) 

Assuming  j  XpL  |  >  |  x^/j,*')  ]  and  xp-  —  —  xp, 
we  find  for  the  constant  potential  difference 

AVclN)~Vc(Xp-,N)—  Vc(xp,N) 

N  S 

=  -r2  (20) 

=  1«=< 

We  note  that  the  above  expression  is  completely 
equivalent  to  the  one  given  by  Eq.  (22)  of  paper  I. 

If  we  consider  the  slab  as  a  part  in  the  infinite 
crystal,  the  above  potential  difference  is  almost  com¬ 
pensated  by  the  result  of  the  crystal,  resulting  in  the 
regular  potential  distribution  of  the  infinite  struc¬ 
ture.  Instead,  we  will  consider  the  slab  free  and  an 
external  potential  to  be  applied,  equal  in  magnitude 
and  opposite  in  sign,  to  the  one  given  by  Eq.  (20), 
which  we  assume  to  drop  uniformly  across  the  slab. 
The  resulting  potential  distribution  is  illustrated  in 
Fig.  1  for_a  slab  of  the  zinc-blende  structure,  parallel 
to  the  (111)  direction,  consisting  of  four  layers. 

In  the  general  case  where  the  unit  cell  may  con- 
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tain  many  atoms  we  have  to  assume  that  these 
atoms  are  distributed  over  several  parallel  (hkl) 
planes.  Let  us  denote  the  position  vector  of  an  atom 
lying  on  the  first  surface  plane  by  x(  lx,)  and  the  po¬ 
sition  vector  of  an  atom  lying  on  the  other  surface 
plane  of  the  slab  by  x(Nk„).  Then  its  thickness  will 
be  given  by 


L  =[x(N,k„)  —  x(  1,Ki)]‘v0 


Having  applied  the  above-mentioned  external  po¬ 
tential,  the  total  potential  at  ion  site  U3k),  apart 
from  the  fluctuating  term  already  encountered  in 
Eq.  (6)  for  y  =  0,  will  be 


Sir  !l  i 1  An  [jt(/j,if')—  x(l,Ki)]-v0  * 

VV3,k)=-^-J,  2  £AZ(  KS)-*(  !,*,)]•  v0+  7-  „  t  ■■  ■  *  I  2i^U's,K)-v0,  (22) 


l\~  1*'“' 


sa  [x(N,k„)—  x(l,/c,)]-v0/>  i ' 


where  the  first  sum  extends  over  all  ions  of  the  unit  cell  from  the  first  surface  plane  up  to  those  U3,k")  lying  on 
the  same  plane  with  ion  (/3,/c).  In  Eq.  (22)  we  assume,  for  reference,  that  the  potential  on  the  first  surface 
plane  is  equal  to  zero. 

The  total  electrostatic  energy  of  the  ions  in  a  long  unit  cell  of  the  slab,  due  to  the  above-described  macro¬ 
scopic  potentials,  will  be 

ifc,Kt/,.«).  (23) 

/j  =  1  l 

and  the  interaction  force  constant  between  the  planes  (/3,k)  and  ( / 3  ,*'),  due  to  their  vibrations  as  solid  planes, 
is  given  by 

~  ,  (24) 

dxa(li,K)dxp(l},K  ) 

or,  by  using  Eqs.  (22)  and  (23)  and  taking  the  second  derivative, 


-y-^-cosacos/S+^^/j.^/j.K')  , 

SaLt 


where 


AFe(JV) 


AVC{N) 


^  (8/jlS|t«r,  ^liN^KKH  )(8/'  j .V ^  cosa  COS/8 


and  where  5’s  are  Kronecker’s  8’s. 

From  Eqs.  (25a)  and  (25b)  it  is  evident  that  for  a 
given  slab  there  exists  a  constant  long-range  interac¬ 
tion  between  any  pair  of  ion  planes,  depending  on 
the  thickness  of  the  slab  [first  term  in  Eq.  (25a)], 
plus  an  interaction  between  each  surface  plane  and 
the  rest  of  the  slab.  This  latter  interaction  given  by 
the  term  **>a'Jjf( ^3 ,#f')  in  Eq.  (25a)  is  a  conse¬ 
quence  of  the  change  in  the  thickness  of  the  slab  due 
to  vibrations  of  the  surface  planes. 

We  should  note  that  the  interaction  between  any 
pair  of  ion  planes  exists  even  if  we  consider  the  slab 
as  a  part  of  an  infinite  crystal  and  as  we  shall  see 
below  it  gives  an  alternative  for  the  calculation  of 
the  macroscopic  field  following  a  two-dimensional 
summation  method.  For  a  slab  which  consists,  in 
part,  of  a  different  structure  (e.g.,  heterostructure), 


the  interaction  of  Eq.  (25a)  has  to  be  calculated  on 
the  basis  of  a  more  realistic  assumption. 

The  total  force  on  the  ion  plane  U3k)  due  to  the 
displacements  of  all  other  ion  planes  is  (a  com¬ 
ponent) 

FaU3,K)=-  2  (26a) 


Fadi,K)=- 


^-^fCosa  2  ^cos/5up(/j,»r') 
tL 


•w’.V'.’V'- 


V'-.; v  -j 


-  2  <t>tfUJ,K;r3,K')uf}{i'3,K') .  <26b) 

r}y,0  v*-;!  v*v- 

To  make  a  comparison  with  the  three-dimensional 
case,  we  take  the  limit  of  Fa(/3,*c)  for  very  thick  -  -  '' 

slabs.  We  recall  [Eq.  (25)  of  paper  I]  that  *  :/~r 
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L=(N-\)aL+d,  (27) 

where  aL  is  the  thickness  of  a  one-unit-cell  layer  and 
d  is  the  distance  between  the  two  outermost  ion 
planes  belonging  to  the  same  unit  cell.  For  N—+  eo 
we  can  introduce  the  cyclic  boundary  condition, 
neglecting  the  deviations  near  the  surfaces,  and  we 
can  put 

ua(/),K)  =  ua(Kkxp[  —  2jr/y1*x(/3)]  ,  (28a) 

where  y3  is  the  wave  vector  perpendicular  to  the 
plane  of  the  slab.  Then  for  long-wavelength  vibra¬ 
tions 

yj->0  and  «a(/3, *)=«<,(*)  .  (28b) 

Taking  the  limit  of  the  first  term  in  Eq.  (26b)  for 
IV  —  oo ,  we  get 

lim  Fa(/3,K)=Fa(ic) 

■V  — *  00 

Air  N  cosa  . 

=  lim  ! - — — - 

.v  -x  j  i„  (N  —  ])aL+d 

(tc')*v0  (29) 

k' 

since  the  second  term  in  Eq.  (26b)  tends  obviously  to 
zero,  and  finally 

F„(k)~  -  —  cosa^2p(K')-v0  ,  (30a) 

va  K' 

where  p  (#c')  is  the  point  dipole  at  ion  site  *c'. 

Hence  the  total  force  on  ion  k  is 

F(»c)  =  fKEmac,  (30b) 

where 

En’iic=  —  —  v0v0'2p(ic)  .  (31) 

va  k 

The  field  of  Eq.  (31)  is  exactly  the  macroscopic  field 
of  Eq.  ( 1 8)  in  the  limit  of  zero  wave  vector  from  the 
direction  perpendicular  to  the  plane  of  the  slab. 

It  is  also  possible  to  obtain  a  more  general  expres¬ 
sion  for  the  interaction  constant  between  ion  planes 
of  the  slab  for  wave  vector  y^tO,  in  the  same  sense 
as  in  Eq.  (25a),  if  necessary.  This  can  be  done  by  us¬ 
ing  the  ih  i  /i  2 )  =  (00)  term  of  Eq.  (6)  where  the  dou¬ 
ble  differential  and  the  sum  over  (3  is  neglected, 
Ppil'ssK')  is  replaced  by  £*•,  and  a  summation  is  per¬ 
formed  over  all  U'3,k')  in  order  to  take  an  appropri¬ 
ate  expression  for  the  potential. 

IV.  DISCUSSION 

We  have  calculated  the  Coulomb  interaction  in 
the  case  of  ionic  slabs  of  arbitrary  orientation  in  any 


structure.  As  far  as  interactions  between  plane  lat¬ 
tices  on  different  planes  are  concerned  the  infinite 
two-dimensional  sums  in  direct  space  are  converted 
into  sums  in  reciprocal  space  only,  which  converge 
rapidly,  depending  on  the  distance  of  the  plane  lat¬ 
tices.  In  the  case  of  a  slab  of  the  zinc-blende  struc¬ 
ture  parallel  to  the  (111)  plane,  ten  to  twelve  terms 
are  needed  for  an  accuracy  of  the  order  of  10~6  (in 
units  of  e1/va)  to  calculate  the  interaction  of  the 
closest  plane  lattices,  while  three  of  four  terms  are 
enough  to  calculate  the  interaction  of  more  distant 
planes.  On  the  other  hand,  this  interaction  drops  off 
very  rapidly  with  the  distance  between  the  interact¬ 
ing  planes,  becoming  about  8  orders  of  magnitude 
weaker  for  the  fourth  neighborhood,  for  zero  wave 
vector,  in  the  above  example. 

For  interactions  between  plane  lattices  lying  on 
the  same  plane,  a  two-dimensional  version  of  the 
Ewald’s  method  of  summation  is  used.  Comparing 
our  results  with  the  formulas  given  by  Tong  and 
Maradudin  for  NaCl  we  find  that  in  case  A,  they 
are  practically  the  same,  while  in  case  B  they  are 
different  since  Mackenzie’s  method  sums  both,  in 
direct  (over  infinite  ion  lines)  and  reciprocal  space 
(one-dimensional  transform  for  each  line).  We  con¬ 
sider  the  formulas  given  here  as  accurate  and  of  easy 
and  direct  use  in  any  case.  We  also  note  that  the 
dependence  of  the  interaction  coefficients  of  Eq.  (14) 
on  the  dividing  point  R  is  unimportant.  A  value  of 
R  equal  to  ( tr/sa )' n  is  quite  adequate  without  any 
further  test. 

Finally  we  should  add  that  the  given  formulas  for 
the  interaction  coefficients  contain  the  macroscopic 
field.  If,  for  instance,  one  uses  the  present  formulas 
to  calculate  the  Coulomb  interaction  in  the  usual 
single-cell  approach  of  the  infinite-crystal  case,  sum¬ 
ming  in  direct  space  as  implied,  what  will  be  ob¬ 
tained  are  the  longitudinal  solutions  for  the  corre¬ 
sponding  direction.  A  two-dimensional  summation 
method,  such  as  the  present  one,  implies  the  pres¬ 
ence  of  surfaces  even  at  infinity,  and  does  not 
separate  out  explicitly  the  macroscopic  field.  A  slab 
of  infinite  thickness  is  not  equivalent  to  an  infinite 
crystal,  since  it  presents  an  infinite  potential  differ¬ 
ence  between  its  two  surfaces.  Compensation  of  this 
potential  results  in  a  definite  field.  The  variation  of 
that  field,  due  to  vibrations  of  the  parallel  to  the 
surface’s  ion  planes,  is  the  well-known  macroscopic 
field  responsible  for  the  higher  frequency  of  longitu¬ 
dinal  modes  in  polar  crystals.  The  calculation 
presented  above  shows  how  this  field  could  be  calcu¬ 
lated  in  more  complicated  cases,  such  as  the  recently 
developed  heterostructures. 
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Erratum:  Lattice  dynamics  of  thin  ionic  slabs.  II.  The  long-range  forces 
[Phys.  Rev.  B  2|,  3398  (1983)} 

G.  Kaneiiis,  J.  F.  Morhange,  am)  M.  Balkansld 


Please  make  the  following  changes  to  our  paper. 

(1)  Page  3399,  in  the  first  line  of  Eq.  (6),  instead  of 

2  2-  ■ 

iyK  *iy*  * 

<j.«  .<j.« 


2  2- 

i-.'  f 


(2)  Page  3400,  in  Eq.  (9),  do  the  same  as  above. 

(3)  Page  3400,  in  the  first  line  of  Eq.  (10),  instead  of 


4ir*  y 

i.  M, 

(4)  Page  3401,  in  Eq.  (13),  the  quantity  in  the  large  parentheses  must  read 

2Pa(/J> K\h, k"I7)  +p0{I}. K )Qm0USl «;/j,  s|y)J  . 

(5)  Page  3401,  second  line  of  Eq.  (14),  instead  of 

ill 

Rs,  ‘ 

write 

4 wi 

Rs, 

(6)  Page  3402,  in  Eqs.  (19a),  (19b),  and  (20),  the  sum  over  k'  runs  from  1  to  the  number  of  ions  in  the  unit  cell  (not  to 

AO. 

(7)  Page  3402,  10th  line  from  the  bottom  of  the  right  column,  instead  of  ...  by  the  result  of  the  crystal  ....  write  ...  by 
the  rest  of  the  crystal .... 

(8)  Page  3403,  in  Eq.  (22),  delete  N  on  the  rightmost  summation  symbol  (over  *'). 

(9)  Page  3403,  in  Eq.  (23),  delete  N  on  the  summation  symbol  over  k. 

(10)  Page  3403,  last  two  lines  of  the  left  column,  instead  of  . . .  which  consists,  in  part,  of  a  ... ,  write  . . .  which  consists 
part  of  a  ...  . 
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Lattice  dynamics  of  thin  ionic  slabs.  III.  Application  to  GaAs  slabs 
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The  lattice  vibrations  of  thin  slabs  (up  to  23  layers)  of  GaAs  are  calculated  for  zero  wave 
vector,  on  the  basis  of  a  rigid-ion  model  fitted  on  the  phonon  dispersion  curves  of  the  infi¬ 
nite  crystal.  It  is  shown  that  all  the  modes,  except  the  surface  ones,  fall  on  the  branches  of 
the  dispersion  curves  of  the  infinite  crystal,  while  the  surface  modes  seem  to  be  combina¬ 
tions  of  those  missing  near  the  zone-edge  modes.  Comparison  is  made  with  former  calcula¬ 
tions  and  the  influence  of  the  short-range  interaction  near  the  surfaces  and  that  of  possible 
macroscopic  fields  is  discussed. 


I.  INTRODUCTION 

Normal  modes  of  vibrations  of  thin  ionic  slabs 
have  been  investigated  theoretically  by  many  au¬ 
thors,1-9  mainly  on  the  basis  of  a  rigid-ion  model, 
and  the  existence  of  surface  modes  have  been  well 
established.  Jones  and  Fuchs7  calculated  the  unre¬ 
tarded  modes  of  a  thin  NaCl  slab;  they  developed  a 
theory  for  the  infrared  optical  properties  of  ionic 
slabs  and  discussed  previously  published  results  on 
surface  modes.  Since  most  of  the  work  mentioned 
above  has  been  done  on  NaCl  slabs  a  more  substan¬ 
tial  comparison  is  possible  between  their  results. 
Benedek9  calculated  surface  dispersion  curves  and 
phonon  densities  for  thin  ionic  slabs  on  the  basis  of 
a  breathing  shell  model  using  Green’s-function  for¬ 
malism.  This  approach  suitably  extended  proves  to 
give  very  good  results  compared  to  those  from  direct 
calculations  on  thin  slabs. 

A  very  general  approach  to  the  effect  of  surfaces 
on  the  vibrational  modes  of  crystalline  solids  is 
given  by  Feuchtwang10  based  on  the  assumption  of 
finite-range  interaction.  Although  in  ionic  crystals 
long-range  electrostatic  interaction  is  important, 
infinite-range  forces  arise  only  in  the  presence  of  a 
macroscopic  field.  Hence  his  results  must  be,  in 
general,  valid.  Different  features  may  appear  to  ap¬ 
ply  only  to  modes  depending  strongly  on  forces  due 
to  such  fields.  The  influence  of  a  macroscopic  field 
on  the  vibrational  modes  of  a  slab  will  be  discussed 
briefly  in  the  last  paragraph  of  Sec.  IV. 

Dispersion  relations  for  surface  modes  are  given 
by  Fuchs  and  Kliewer1  who  found  all  long- 
wavelength  optical  modes  of  an  ionic  slab  to  have  ei¬ 
ther  TO  or  LO  frequencies.  Tong  and  Maradudin3 
treated  in  detail  the  case  of  a  NaCl  slab  parallel  to 
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the  (100)  plane  and  pointed  out  the  importance  of 
some  approximations  made  in  the  former  work. 
Further  comments  on  this  point  are  given  by  Jones 
and  Fuchs.7  We  also  found  the  proper  description 
of  the  short-range  interaction  between  atoms  near 
the  surface  to  be  of  fundamental  importance  in  cal¬ 
culating  the  correct  frequencies  of  the  surface 
modes. 

Among  the  considerable  experimental  works  on 
infrared  absorption  or  Raman  scattering  on  thin 
films  and  powders,  which  show  vibrational  states  ei¬ 
ther  between  the  TO  and  LO  frequencies  of  the  in¬ 
finite  crystal  or  below  the  TO  frequency,  we  men¬ 
tion  the  infrared  transmission  and  reflection  mea¬ 
surements  on  thin  (up  to  68  pm)  films  on  GaAS  by 
Cochran  et  o/.“  and  Fray  et  a/.*2  They  observed  in 
a  25  pm  thin  film,  lattice  absorption  on  the  TO  and 
LO  frequencies  and  also  two  other  strong  peaks  on 
frequencies  lying  between  these  two.  Each  of  those 
peaks  seems  to  consist  of  three  or  four  fine-structure 
peaks.  These  features  are  attributed  to  the  size, 
shape,  and  orientation  of  the  specimens  used. 

Raman  spectra  on  laser-annealed  GaAs  (Ref.  13) 
show  a  gradual  transition  from  the  amorphous  spec¬ 
trum  to  the  crystalline  one  involving  a  continuous 
variation  of  degree  of  order  versus  the  annealing 
density  energy.  The  above  experimental  results  and 
also  analogous  results  on  other  materials  reveal  the 
need  for  a  complete  calculation  of  infrared  and  Ra¬ 
man  spectra  of  small  crystallites  in  order  to  justify 
the  hypothesis  of  size  effects. 

In  the  present  paper  we  report  on  calculations  of 
vibrational  modes  of  thin  GaAs  slabs  parallel  to 
(111)  planes,  on  the  basis  of  a  rigid-ion  model.  This 
model  has  been  used  by  Kune14  to  fit  measured  pho¬ 
non  dispersion  curves,  hence  our  results  are  directly 
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comparable  to  the  ones  of  an  infinite  crystal.  The 
long-range  interaction  has  been  calculated  in  a  form¬ 
er  paper'5  referred  to  as  paper  II. 

In  the  next  section  we  briefly  review  the  equations 
of  motion  and  the  resulting  dynamical  matrix  to  be 
solved.  In  the  third  section  we  describe  in  short  the 
rigid-ion  model  used  for  the  calculation  and  the  pro¬ 
cedure  applied  to  calculate  the  short-range  interac¬ 
tion  of  atoms  lying  near  the  surfaces.  Finally,  in  the 
last  section  we  describe  the  results  obtained  and  we 
discuss  them  in  comparison  with  other  results. 


II.  EQUATIONS  OF  MOTION 
AND  DYNAMICAL  MATRIX 


We  consider  a  crystal  slab  of  zinc-blende  structure 
(particularly  GaAs)  parallel  to  the  (111)  plane. 
Choosing  a  Cartesian  coordinate  system  Ox  tx2x3, 
whose  axes  are  parallel  to  the  edges  of  the  cubic  fee 
unit  cell  and  its  origin  on  an  ion  site  (for  instance,  a 
Ga  site),  the  primitive  translation  vectors  are 


*1 

«2 
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~  2 

0  1  1 

1  0  1 

*01 

*02 

83. 

1  1  0 

*03 

where  a  is  the  lattice  constant.  Vectors  a[  and  a2  lie 
on  the  plane  (111)  while  a3  lies  out  of  it. 

The  equations  of  motion  for  a  lattice  are16 

m„ua(l,K)=-'2t'2l<l>afif,l,KU  ,K')UpU',K')  , 

/',«■  p 

(2) 

where  m,  is  the  mass  of  the  »rth  kind  of  ion,  ua{l,K ) 
is  the  ath  Cartesian  component  of  the  displacement 
from  the  equilibrium  position,  ^a/j( /,*;/',*')  are  the 
atomic  force  constants,  and  /  =(/ i,/2,/3)  labels  the 
unit  cells. 

Applying  cyclic  boundary  conditions  along  the 
directions  3]  and  a2  [on  the  infinite  (111)  plane]  we 
put 


UaU,K)  = 


Va(litK) 


m 


1/2 


exp[  —  ia)t  +  27riy-j?(/,*)]  , 


(3) 


where  y  is  a  two-dimensional  wave  vector.  Equa¬ 
tion  (2)  becomes 


W2Ua(/3,K)=  2  ^LDap(li>K\l\K'  I  y)Vp(r3K')  , 

I\y  p 


(4) 


2VWi,*'l7)=  2; 


Xexp[2my-T(/,<c;/V)] 

(5) 

are  the  elements  of  the  6NX6N  dynamical  matrix 
(N  is  the  number  of  layers  in  the  slab).  Since  no 
periodic  boundary  condition  is  used  along  the  finite 
dimension,  any  further  reduction  of  the  dynamical 
matrix  will  be  likely  possible  only  from  symmetry 
considerations. 

It  has  been  shown17  that  for  the  above  chosen 
orientation  of  the  slab  one  can  use  a  new  coordinate 
system  0x',x2x3,  related  to  the  old  one  associated 
with  the  crystallographic  unit  cell,  by  the  transfor¬ 
mation 
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whose  Ox  |  and  Ox  2  axes  are  coplanar  with  the  slab, 
while  axis  Ox\  is  perpendicular  to  it.  Since  it  is  al¬ 
ways  possible  to  find  such  a  coordinate  system, 
whatever  the  orientation  of  the  slab  might  be,  let  us 
denote  the  corresponding  transformation  matrix  by 
H  [in  the  present  case  H  is  the  3x3  unitary  matrix 
used  in  Eq.  (6)]. 

With  the  use  of  transformation  H  and  the  repre¬ 
sentation  of  the  space  group  G  of  the  three- 
dimensional  structure  with  respect  to  the  old  coordi¬ 
nate  system,  one  can  construct  a  new  space  group  G' 
appropriate  for  the  two-dimensional  structure  of  the 
slab,  by  transforming  the  representation  of  G  to  the 
new  coordinate  system  and  by  picking  up  those  ele¬ 
ments  which  act  only  parallel  to  the  plane  of  the 
slab. 

Group  G'  can  then  be  used  to  provide  the  form  of 
the  force-constant  matrices  and  the  relations  be¬ 
tween  the  elements  of  the  dynamical  matrix.  In  the 
present  case  and  for  wave  vector  y  =  0,  we  find  that 
the  interaction  between  the  plane  lattices  assumes 
the  general  form 


fi(/3,<c;/3,x')= 


A  B  —B 
B  A  -B 
-B  -B  A 


(7) 


in  the  old  coordinate  system  0X|X2x3,  while  in  the 
system  0xix'2x3  it  takes  the  form 


fi’(/3,*;/3,*')  = 


A  —B  0  0 

0  A  -B  0 
0  0  A+2B 


(8) 


where 


>  f 
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through  the  transformation 

The  form  of  Eq.  (8)  implies  that  for  the  so-elmscn 
orientation  of  the  slab,  the  solutions  of  l:.q.  (4)  fall 
into  two  groups,  one  doubly  tlegcneralr  comi.sling  of 
the  in-plane  solutions  (x-y  nuxlcs)  ami  one  nondo- 
generate  consisting  of  the  out-of-plane  solutions  (2 
modes). 

The  matrix  whose  elements  arc  defined  by  l:.q.  (5) 
expresses  the  interaction  between  plane  lattices  and 
has  to  be  calculated  on  the  basis  of  some  model.  In 
the  next  section  we  give  a  brief  description  of  the 
applied  rigid-ion  model.  Since  both  short-  and 
long-range  forces  are  taken  info  account,  it  is  cus¬ 
tomary  to  consider  the  dynamical  matrix  as  consist¬ 
ing  of  two  parts,  the  short-range  (sr)  and  the 
Coulomb  parts; 

QUuX'Jf  —  ,k) 

+QcU},kJ'},k)  .  (10) 

The  form  of  the  interaction  matrices  between  the 
plane  lattices  of  the  slab  for  zero  wave  vector,  is 
given  in  the  Appendix. 

The  Coulomb  part  can  be  expressed  in  terms  of 
the  Q  coefficients  calculated  in  paper  II  [Eqs.  (10) 
and  (14)],  and  for  the  general  wave  vector  it  assumes 
the  form 

x  2  S*'Qaf>lli,Kin,K"  I  0) 

(mKmK-V'1 


(11) 

where  £*  is  the  charge  fraction  attributed  to  ion  k. 
Values  of  the  Q  coefficients  for  zero  wave  vector  are 
also  given  ir.  the  Appendix  (Table  I). 


TAIll.K  I.  Values  of  die  Coulomb  coefficient  (i  [in 
iiuils  of  C/.r)' 


r 


a 


1 

1 

1 

2 
2 
1 

1 

2 
2 


2 

1 

2 

1 

2 

1 

2 

1 

2 


0 

1 

1 

1 

1 

2 
2 


5.591  35 
0.14341 
0.03349 
1.43762 
0.14341 
0.00040 
0.00018 


2  -0.00175 

2  -0.00040 


III.  THE  MODEL 

In  the  present  case  we  use  the  rigid-ion  model 
(RIM)  developed  and  applied  to  several  binary  com¬ 
pounds  of  the  zinc-blende  structure  by  Kune.14 
Apart  from  the  effective  charge  q*  the  model 
parameters  are  ten  tensorial  force  constants,  A  and  B 
for  first-neighbor  central  and  noncentral  interaction, 
Ci,  D 1,  Ei,  Fu  C2,  D2,  E2,  and  F2,  for  second- 
neighbor  central  and  noncentral  interaction,  for  the 
two  different  kinds  of  ions. 

All  of  the  above  parameters  have  been  fitted  to 
experimentally  known  phonon  dispersion  curves,  the 
elastic  constants,  and  the  piezoelectric  constant.  In 
order  to  use  the  above  model  in  the  case  of  a  slab  the 
following  adaptations  have  to  be  made: 

(a)  Long-range  forces  have  to  be  recalculated  on 
the  basis  of  new  suitable  formulas  as  already  men¬ 
tioned. 

(b)  Short-range  interaction  of  the  near-the-surface 
atoms  has  to  be  modified,  so  as  to  take  the  missing 
ions  into  account.  For  the  short-range  part  of  the 
dynamical  matrix  the  following  procedure  has  been 
followed.  From  the  ten  fitted  tensorial  force  con¬ 
stants,  the  values  of  a  set  of  ten  valence  force  field 
(VFF)  force  constants  are  deduced,  namely  the  £,  A., 
p,  p,  a,  v,  kg,  kg,  kn,  and  k'„  according  to  the  fol¬ 
lowing  model  for  the  potential  energy: 


<*>— ro%  2  Aty  +  —  2  (^rij)Z  +  r\P  2  Ar<*  + 

Ga-As  Ga-As  Ga-Ga 

2  2  <Afy*>2 

L  AvGa-Aa  L  Ga-As-Ga 

where  r0  and  r,  are  the  first-  and  second-nearest- 
neighbor  distances,  respectively. 

The  interaction  of  the  near-the-surface  ions  is  cal¬ 
culated  on  the  basis  of  the  above  VFF  model.  It 


2  (Ar,*)2+r,o-  2  &rjl  +  T  2  <Ar,,)2 

Ga-Ga  As-As  2  As-As 

~~  2  A/>Ar,7  +  —Z~  2  ArijbrJk  ,  (12) 

^  As-Ga-As  Ga-As-Ga 

r— - — - 

should  be  noted  that  the  ten  VFF  parameters  are  not 
independent,  but  have  to  fulfill  the  equilibrium  con¬ 
dition.  For  the  given  values  of  the  tensorial  force 
constants  this  condition  does  not  hold.  Hence  we 
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FIG.  1.  Variation  of  the  frequencies  of  six  modes  vs 
the  thickness  of  the  slab. 


decompose  the  tensorial  force  constants  into  the 
VFF  parameters  disregarding  this  condition. 

By  so  doing,  only  the  interaction  of  atoms  lying 
near  the  surfaces  may  be  influenced.  Since  the 
correct  description  of  the  above  interaction  is  of 
critical  importance  in  calculating  the  surface  mode 
frequencies,  one  needs  a  faithful  VFF  model  in  or¬ 
der  to  obtain  unambiguous  results  for  these  modes. 
In  the  next  section  we  comment  on  the  results  based 
on  the  above  consideration  and  on  those  obtained 
when  the  short-range  interaction  of  surface  atoms  is 
considered  to  be  the  same  as  for  atoms  in  the  interi¬ 
or  of  the  slab. 

IV.  RESULTS  AND  DISCUSSION 

We  have  calculated  the  frequencies  and  the  eigen¬ 
vectors  for  thin  slabs  of  GaAs  oriented  parallel  to 
the  (111)  plane  with  thicknesses  from  2  up  to  25 
cells  (10—130  A  thick)  for  zero  (two-dimensional) 
wave  vector. 

The  main  feature  displayed  by  the  solutions  is 
that,  by  increasing  the  number  of  layers  of  the  slab 
(i.e.,  its  thickness)  all  of  the  solutions  tend  rapidly  to 
certain  limits,  while  the  new  solutions  appearing  in 


each  thicker  slab  are  interpolated  between  the  solu¬ 
tions  of  the  thinner  one.  For  a  slab  20  layers  thick 
all  the  frequencies  are  within  \%  of  their  limits. 

In  Fig.  1  we  show  the  variation  of  the  frequencies 
of  six  modes  versus  the  thickness  of  the  slabs:  the 
highest-frequency  z  mode  (first  LO),  the  highest- 
frequency  x-y  mode  (first  TO),  and  four  surface 
modes,  i.e.,  two  surface  x-y  modes  (TS1  and  TS2) 
and  two  surface  z  modes  (LSI  and  LS2).  We  see 
that  the  surface  modes  TS 1  and  TS2  have  practical¬ 
ly  constant  frequency  for  any  slab  as  thin  as  three 
layers,  while  the  frequencies  of  the  rest  of  the  modes 
tend  very  rapidly  to  their  limiting  values. 

The  two-dimensional  Brillouin  zone  correspond¬ 
ing  to  the  slab  structure  is  a  section  through  the 
center  of  the  three-dimensional  zone  of  the  fee  lat¬ 
tice,  perpendicular  to  the  A  direction  (—  £,  —  £,  £). 
Comparing  the  frequencies  of  the  above  modes  with 
the  frequencies  of  the  modes  belonging  to  the  disper¬ 
sion  branches  of  the  A  direction  of  the  infinite  crys¬ 
tal,  we  find  that  the  first  LO  and  first  TO  modes 
tend  to  have  frequencies  equal  to  those  of  the  LO 
and  TO  modes  of  the  T  point  of  the  infinite  struc¬ 
ture,  respectively.  The  surface  mode  TS  1  has  a  fre¬ 
quency  almost  equal  to  the  TO  mode  at  the  L  point, 
while  modes  TS2,  LSI,  and  LS2  have  frequencies 
which  fall  into  the  gap  at  the  L  point  of  the  infinite 
crystal. 

Before  making  any  further  comments  on  the  sur¬ 
face  modes  and  on  the  rest  of  the  modes,  we  will 
turn  our  attention  to  the  eigenvectors.  For  a  GaAS 
slab  N  layers  thick,  there  are  2 N  x-y  modes,  which 
are  doubly  degenerate  (the  transverse  modes),  and 
IN  z  modes  (the  longitudinal  ones).  From  the  ?\ 
modes  in  each  configuration,  N  are  optical  ana  the 
remaining  N  are  acoustic.  Among  the  .V  optical 
modes  in  each  case  there  are  two  whose  amplitudes 
decay  exponentially  along  the  finite  dimension  of 
the  slab,  from  the  one  surface  to  the  other,  and  are 
therefore  called  surface  mo  '  *s.  For  wave  vectors 
different  from  zero  the  so-called  Rayleigh  waves16 
appear,  among  the  acoustic  modes.  For  each  optic 
mode  m  we  construct  the  difference  of  the  reduced 
displacements, 


u(l  3,1  ( m) 

Vm, 


u(l},2  |  m ) 

VWi 


(13) 


where  1  denotes  the  Ga  atom,  2  the  As  atom,  /, 
numbers  the  layers  of  the  slab,  and  M  is  the  mass. 
Accordingly,  for  the  acoustic  modes,  we  construct 
the  sum  of  the  same  quantities. 


„  u(/,,  1  |  m)  «(/ 3,2  |  m ) 

'S'-'— ✓sr+-?sr 


(14) 


In  Fig.  2  we  plot  /£P(/.J  and  f£U ,)  vs  z 
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FIG.  2.  For  a  slab  of  ten  layers:  (a)  Relative  reduced 
ionic  displacements  for  some  x-y  optic  modes.  Number¬ 
ing  starts  with  the  mode  of  highest  frequency,  (b)  Sum  of 
reduced  ionic  displacements  for  some  x-y  acoustic  modes. 
Numbering  starts  with  the  zero-frequency  mode. 


layor 

FIG.  3.  For  a  slab  of  ten  layers:  (a)  Relative  reduced 
ionic  displacements  for  some  z  optic  modes.  Numbering 
starts  with  the  mode  of  highest  frequency,  (b)  Some  of  re¬ 
duced  ionic  displacements  for  some  z  acoustic  modes. 
Numbering  starts  with  the  zero-frequency  mode. 


the  coordinate  along  the  finite  dimension  of  the  slab, 
for  some  x-y  modes  (transverse  modes),  and  in  Fig.  3 
we  plot  the  same  functions  for  some  z  modes  (longi¬ 
tudinal).  It  is  evident  from  these  figures  that  all  the 
optic  modes,  except  the  surface  ones,  can  almost  be 
described  (neglecting  the  sign)  as 

/°»(z)^^ sin-^p-  ,  (15) 

where  L  is  the  thickness  of  the  slab,  while  the  acous¬ 
tic  modes  can  be  described  as 


,,c,  i  2  (m  —  l)irz 

f”iz)~VW C  - L -  6 

(m  =  1,2, . . . ,  N,  numbers  the  optic  modes  starting 
from  the  highest-frequency  one  and  the  acoustic 
modes  starting  from  the  zero-frequency  one). 

We  note  that  the  above  simple  trigonometric  ex¬ 
pressions  do  not  describe  exactly  the  functions  of  the 
displacement  defined  by  Eqs.  (13)  and  (12),  but  rath- 
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FIG.  4.  Phonon  dispersion  curves  of  GaAs  along  the  A 
direction  (solid  lines),  x-y  modes  (solid  circles)  and  z 
modes  (open  circles)  for  a  slab  of  20  layers  thick.  Broken 
lines  show  the  position  of  surface  modes.  The  positions  of 
TO  and  LO  missing  modes  near  and  on  the  zone  boun¬ 
dary  are  shown  by  arrows. 


er,  give  a  good  picture  of  what  these  functions  look 
like.  Moreover,  although  the  functions  sketched  in 
Figs.  2  and  3  refer  to  a  slab  ten  layers  thick,  the  pic¬ 
ture  is  the  same  for  slabs  of  any  thickness. 

From  Eqs.  (15)  and  (16)  we  see  that  there  is  an 
implicit  dependence  of  the  modes  on  a  wave  vector 
along  the  finite  dimension  of  the  slab,  although  such 
an  assumption  has  not  been  made.  All  the  modes, 
except  the  surface  ones,  on  which  we  will  comment 
later  on,  seem  to  be  characterized  by  the  values  of  a 
wave  vector 


There  is  no  optic  mode  with  n  =0,  i.e.,  there  is  no 
mode  of  infinite  wavelength.  Since  the  surface 


modes  cannot  be  described  by  the  same  function  /£p, 
only  N  —  2  optic  modes  could  be  attributed  to  the 
wave  vectors  of  Eq.  (17).  On  the  other  hand,  all  the 
N  acoustic  modes  can  be  regarded  as  corresponding 
to  some  wave  vector  yz.  We  should  note  that  the 
solution  of  zero  frequency  for  y2  = 0  has  been  im¬ 
posed  by  using  relation  (A5)  when  calculating  the 
self-terms  of  the  dynamical  matrix. 

According  to  the  above  observations,  we  plot  in 
Fig.  4  the  dispersion  curves  of  GaAs  along  the  A 
direction,  as  they  have  been  calculated  by  Kune14  on 
the  basis  of  the  same  model  and  we  put  on  the  same 
graph  the  solutions  for  a  slab  20  layers  thick.  We 
have  chosen  this  thickness  where  each  solution  has 
reached  its  limit  within  less  than  1%.  Acoustic 
modes  start  from  zero  wave  vector,  while  optic 
modes  start  from  wave  vector  equal  to  0.5/20.  We 
see  that  all  the  modes,  except  the  surface  ones,  fall 
on  the  corresponding  branches:  x-y  modes  on  the 
transverse  branches  and  z  modes  on  the  longitudinal 
ones.  There  are  four  optic  modes  missing.  Instead 
there  are  four  surface  modes  shown  by  dotted  lines 
on  the  graph  along  the  zonewidth,  although  they 
seem  to  belong  to  a  value  of  the  wave  vector  very 
close  to  the  zone  boundary.  Hence  we  conclude  that 
the  surface  modes  could  be  regarded  as  superposi¬ 
tions  of  the  missing  optical  modes,  and  perhaps  the 
corresponding  acoustic  modes  of  the  zone  edge,  in 
each  configuration  (transverse  or  longitudinal).  The 
longitudinal  surface  modes  (z  modes),  considered  as 
being  superpositions  of  modes  of  higher  frequency, 
appear  to  have  lower  frequencies  than  the  transverse 
ones.  They  also  show  considerably  less  decay  along 
the  finite  dimension  of  the  slab.  Both  of  these 
features  could  be  explained,  qualitatively  at  the  mo¬ 
ment,  by  the  weaker  short-range  forces  near  the  sur¬ 
faces  and  by  the  assumption  that  these  modes  are 
strongly  damped  by  the  long-range  forces.  We 
should  note  at  this  point  that  the  slowest  conver¬ 
gence  of  the  frequencies  to  their  limiting  values, 
with  increasing  thickness  of  the  slab,  is  observed  for 
the  longitudinal  modes  of  shorter  wavelength  (wave 
vector  near  the  zone  edge).  Indeed,  the  frequencies 
of  the  surface  modes  depend  very  strongly  on  the 
short-range  forces  assumed  for  the  surface  layers.  If 
we  restore  on  these  layers  the  same  short-range 
forces  as  for  the  rest  of  the  slab,  then  all  the  modes, 
except  the  surface  ones,  tend  faster  to  the  same  lim¬ 
iting  frequencies.  The  surface  modes  have  higher' 
limits  as  follows:  The  TS1  mode,  253.3  cm-1  (in¬ 
stead  of  252.5  cm-1),  the  TS2  mode,  247.3  cm-1 
(instead  of  235  cm-1),  the  LSI  mode,  250.7  cm-1 
(instead  of  230.7  cm-1),  and  the  LS2  mode,  249.6 
cm-1  (instead  of  210.5  cm-1).  We  see  that  the  long¬ 
itudinal  surface  modes  have  a  much  stronger  depen¬ 
dence  on  the  short-range  forces  assumed  for  the  sur- 
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face  layers,  but  still  at  least  one  of  them  should  have  cussed  by  Tong  and  Maradudin,3  we  would  like  to 

higher  frequency  than  the  transverse  modes  if  no  note  that  their  results  for  “bulk”  modes  concern 

damping  due  to  long-range  forces  was  existing.  those  modes  which  lie  on  the  corresponding 

Hence  the  correct  description  of  the  short-range  branches  but  very  near  to  the  T  point.  Hence  for  a 

forces  near  the  surfaces  is  of  critical  importance  in  thick  slab  there  are  many  modes  which  have  the 

calculating  the  surface  mode  frequencies.  Of  course,  zone-center  TO  and  LO  frequencies  with  small  wave 

the  eigenvectors  are  also  influenced  by  the  change  of  vectors  mv/L.  For  optic  modes  we  have  found 

the  forces  near  the  surface,  in  particular  the  eigen-  only  the  sine  dependence  of  the  displacements  with 

vectors  of  the  longitudinal  modes,  but  they  retain  odd  and  even  values  of  m.  Their  point,  that  this 

their  main  features.  Stronger  forces  between  the  dependence  is  such  that  there  are  exactly  m  half 

surface  layers  result  in  deeper  penetration  of  those  waves  across  the  thickness  of  the  slab,  is  correct 

modes  in  the  slab.  even  for  very  thin  slabs  and  for  the  acoustic  modes 

On  the  basis  of  the  above  observations  it  is  clear  also.  The  quantization  of  the  wave  vector  along  the 

that  the  continuum  of  frequencies  for  “bulk”  modes  finite  dimension  is  of  course  the  same  for  all  disper- 

obtained  by  several  authors  when  calculating  solu-  sion  branches. 

tions  for  different  wave  vectors  on  the  plane  of  a  The  surface  modes  found  by  these  authors  reduce, 

slab,  is  the  projection  of  the  dispersion  surfaces  of  for  zero  wave  vector  on  the  plane  of  the  slab,  in  un- 

the  infinite  crystal  on  the  two-dimensional  Brillouin  localized  surface  modes,  corresponding  to  zone- 

zone  appropriate  for  the  slab.  Comparing  our  re-  center  modes  of  the  infinite  crystal.  It  has  been 

suits  with  those  of  Tong  and  Maradudin,3  who  cal-  shown  by  Feuchtwang10  that  these  modes  are  a  spe¬ 
culated  the  solutions  for  a  slab  of  NaCl  parallel  to  cial  case  of  bulk  modes.  Since  the  above  authors 

(001)  plane,  we  could  note  the  following:  The  fre-  treated  the  case  of  long  waves,  these  modes  may  be 

quency  of  the  transverse  surface  mode  in  NaCl,  ly-  the  only  ones  which,  for  finite  wave  vectors  (#0)  on 

ing  just  below  the  lower  limit  of  the  “bulk”  optical  the  plane  of  the  slab,  became  surface  modes.  Lucas4 

modes,  is  in  agreement  with  our  results.  The  near  and  Jones  and  Fuchs7  attribute  the  above  behavior 

degeneracy  of  both  pairs  of  surface  modes  they  of  these  modes  to  the  neglect  of  the  changes  of  the 

found  can  be  explained  by  the  fact  that  the  two  sur-  forces  acting  on  atoms  near  the  surfaces.  The  re¬ 
faces  of  the  NaCl  slab  in  the  above  orientation  are  suits  of  the  latter  author  at  y  =  0  are  in  agreement 

completely  equivalent.  Moreover,  the  thickness  of  to  those  of  Lucas,  and  Tong  and  Maradudin,  as  to 

this  slab  is  an  integer  multiple  of  the  lattice  con-  the  number  and  type  of  surface  modes, 

stant,  while  in  our  case  it  is  not.  This  feature  may  A  final  remark  concerns  the  possible  influence  of 

be  of  importance  when  superimposing  plane  waves,  macroscopic  fields.  We  have  shown  in  paper  II  that 

due  to  the  phase  difference  it  may  introduce.  As  far  vibrations  involving  ionic  displacements  perpendicu- 

as  it  concerns  the  higher-frequency  surface  modes  lar  to  the  plane  of  the  slab  give  rise  to  a  potential 

for  zero  wave  vector  in  the  case  of  the  NaCl  slab,  difference  between  its  two  surfaces.  This  potential 

they  can  be  considered  as  unlocalized  surface  modes  does  not  influence  the  vibrations  of  a  free  slab,  ex- 

whose  wave  vector  perpendicular  to  the  plane  of  the  cept  if  it  results  in  additional  surface  charges,  or  if 

slab  can  assume  values  either  at  the  center  or  at  the  the  slab  is  considered  in  some  polarizable  environ- 

zone  boundary10  in  the  present  case.  The  corre-  ment.  In  such  cases,  assuming  complete  compensa- 

sponding  modes  in  our  case  are  of  the  same  type  as  tion  of  the  above  potential,  all  the  frequencies  of  the 

those  found  by  Wallis.18  This  difference  may  be  due  z  modes  tend  slower  to  the  same  limits,  while  an  ad- 

to  the  fact  that  the  NaCl  slab  consists  of  identical  ditional  phase  difference  is  introduced  in  the  dis- 

“neutral”  planes.  placements  of  atoms  in  neighboring  cells.  This 

Fuchs  and  Kliewer1  treated  the  case  of  long-wave  phase  difference  moves  all  the  optic  modes  one  step 

optical  vibrations  in  a  slab  in  the  electrostatic  ap-  toward  the  zone  boundary.  Hence  for  a  slab  thicker 

proximation.  Apart  from  the  influence  on  their  re-  than  ten  layers,  the  relative  displacements  show  the 

suits  due  to  approximations  concerning  the  short-  same  pattern  as  described  by  Eq.  (15),  but  now 

range  forces  near  the  surfaces  and  the  replacement  modes  with  m  =0,  m  =  1,  and  m  =jV  —  1  are  miss- 

of  infinite  sums  by  integrals,  which  have  been  dis-  ing. 


APPENDIX 

The  form  of  force-constant  matrices  (interaction  between  individual  ions)  for  the  zinc-blende  structure  is  the 
following: 
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-B  A 

c, 

o, 

c2 

d2  e2 

4>(/,l  ;/',!)  = 

*>> 

C, 

Ex 

,  $(/,2;/',2)=  D2 

c2  e2 

F, 

-e2 

E2  f2 

where  /  =  ( /j ,/2,/3 )  and  l'  =  (lx,l2,l3  +  1). 

For  the  surface  layers  Iq  =  U\,12,  1 )  and  lN  =(l\,l2,N)  the  force-constant  matrices  in  the  present  case  take  the 
form 


4  |  £ [  fi|  A  3  flj  £4 

d>(/o, l;/o.2)=  fl|  4|  B|  ,  4>(//v,  1  ;//v *2 )  =  £3  ^3  £4 

£2  £2  A2  B 3  £3  45 


The  decomposition  of  the  tensorial  force  constants  into  parameters  of  the  VFF  model  used  can  be  found  in 
Ref.  14  for  all  interactions,  except  the  ones  between  ions  in  the  surface  layers,  which  are  listed  below: 

A\--A—krr/  6,  A2  =A  +  4kg/3  +  /c„.-/6  ,  A2=A — krr' /  6,  A^~Akg/^-\-kn' /(>,  (A3) 

£,=£-*„,/ 6,  B2  =  fi-2k9/3  +  kn./6,  B^B-k^/6,  BX  =  B -2k'e/3,  +  k^ /6  . 

The  interaction  matrices  for  plane  lattices  for  zero  wave  vector  (submatrices  of  the  dynamical  matrix)  as¬ 
sume  the  following  forms  (with  /'3  =  /3  + 1): 

3/1  -B  B  A  B  -£ 

£(/3.1;/j,2)~  — £  3/1  £  ,  D(/3,2;/'3l)~  £  4  -£  ,  (A4a) 


B  B  3A 


\-B  -B  A 


2Ci+£i  Dx  —D\ 

Dx  2Cj-)-£|  —D\ 

,  £(/3,2;/j,2)~ 

2C2  +  F2 

D2 

2  C2+F2 

-d2 

—  D2 

—  D\  ~D\  2Ci+Fi 

-d2 

-d2 

2  C2  +  F2 

Self  terms,  for  zero  wave  vector,  are  calculated  on  the  basis  of  the  equation16 

2^i«-,iy)=o . 

The  corresponding  interaction  matrices  for  planes  of  the  surface  layers  1/3  =  1  and  l3=N)  take  the  form 


24 1  +  42 

-b2 

243+44  —  £4  £4 

fi(l.l;1,2)~ 

-B2 

24 1  +  4  2 

£2  ,  Q(N,l;N,2)~ 

—£4  243+44  £4 

*2 

*2 

24 1  +4  2 

£4  £4  243+44 

All  the  Coulomb  interaction  matrices  assume  for  zero  wave  vector  the  form 
0  0-0 

g(0, *;/’,*')=  p  0  -p  .  (A7) 

-P  -p  0 

The  values  of  0  are  given  in  Table  I  for  /'= 0,  1,  and  2.  All  more  distant  interactions  are  less  than  10' 5  [in 
units  of  ( Ze)2/va ]. 
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Please  make  the  following  changes  to  our  paper. 
(1)  Page  3413.  in  Eq.  (A2),  instead  of 
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Systematic  measurements  by  light  scattering  of  the  linewidth  and  frequency  shift  of  the  q  =0  op¬ 
tical  phonon  in  silicon  over  the  temperature  range  of  5-1400  K  are  presented.  Both  the  linewidth 
and  frequency  shift  exhibit  a  quadratic  dependence  on  temperature  at  high  temperatures.  This  indi¬ 
cates  the  necessity  of  including  terms  in  the  phonon  proper  self-energy  corresponding  to  four- 
phonon  anharmonic  processes. 


I.  INTRODUCTION 

Experimental  studies  of  the  inelastic  scattering  of  light 
by  crystals  have  provided  a  great  deal  of  information  con¬ 
cerning  the  optical  modes  of  vibration  at  the  center  of  the 
Brillouin  zone.  In  pure  materials  one  finds  typically  that 
both  the  line  center  and  the  linewidth  vary  with  tempera¬ 
ture.  This  temperature  dependence  can  be  attributed  to 
the  anharmonic  terms  in  the  vibrational  potential  energy.' 

If  one  restricts  oneself  to  cubic  anharmonicity  in  second 
order,  the  damping  constant  which  characterizes  the 
linewidth  is  proportional  to  the  absolute  temperature  T  in 
the  high-temperature  limit,  but  when  one  includes  quartic 
anharmonicity  to  second  order  and/or  cubic  anharmonici¬ 
ty  to  fourth  order,  the  damping  constant  involves  terms 
proportional  to  T 2  in  the  high-temperature  limit.2  For  the 
case  of  silicon.  Hart,  Aggarwal,  and  Lax3  have  measured 
the  frequency  shift  of  the  line  center  and  the  damping 
constant  over  a  range  of  temperatures  from  20  to  770  K. 
They  found  that  their  data  for  the  frequency  shift  agree 
rather  well  with  the  theoretical  calculations  of  Cowley4 
based  on  cubic  anharmonicity  to  second  order,  but  their 
data  for  the  damping  constant  show  significant  deviations 
from  Cowley's  results.  Hart  et  al.  were  able  to  show, 
however,  that  their  data  for  the  damping  constant  can  be 
fitted  satisfactorily  by  the  cubic  anharmonic  model  of  Kle- 
mens5  if  the  zero-temperature  value  of  the  damping  con¬ 
stant  is  properly  chosen. 


Recently,  Tsu  and  Hernandez*1  have  reported  measure¬ 
ments  of  the  frequency  shifts  of  both  one-phonon  and 
two-phonon  Raman  lines  for  silicon  over  the  temperature 
range  20— 900‘C.  Where  their  results  overlap  with  those 
of  Hart  e‘  al.,  the  agreement  is  good.  No  data  on  the 
linewidth  is  presented  by  Tsu  and  Hernandez. 

In  the  present  paper,  measurements  of  the  light  scatter¬ 
ing  spectrum  of  silicon  are  reported  for  the  temperature 
range  between  5  and  1400  K.  The  temperature  depen¬ 
dences  of  the  frequency  shift  and  damping  constant  of  the 
Raman  active  LO  phonon  are  analyzed  in  terms  of  cubic 
and  quartic  anharmonic  contributions.  It  is  found  that  at 
the  higher  temperatures,  cubic  anharmonic  terms  to 
second  order  are  not  sufficient  to  fit  the  data,  but  the  in¬ 
clusion  of  higher-order  terms  involving  cubic  and/or 
quartic  anharmonicity  makes  possible  a  satisfactory  fit. 

II.  ANHARMONICITY  IN  LIGHT  SCATTERING 
BY  OPTICAL  PHONONS 

For  a  system  whose  equilibrium  atomic  positions  are 
specified  by 

R'°'(/,ic)  =  R(/)  +  R(ic)  ,  (2.1) 

where  R(/)  =  /,f,-f  /2  r j,  the  r,  are  primitive 

translation  vectors,  the  /,  are  integers,  and  R(x)  is  a  vector 
of  the  basis;  the  vibrational  Hamiltonian  can  be  written  as 


H=  7  X  TT'/’al/.^+T  2  2 

le.a  l.e.a  l.e.B 

+  12  2  2 

t.*,o  r.e’.fi  r,*",r 


l,K.a  l\i r.ir *\y 


(2.2) 


where  u  ( /.  )  is  the  displacement  of  atom  I,k  from  its 
equilibrium  position  and  <t>ag,  and  <1 are  the  har¬ 
monic,  cubic  anharmonic,  and  quartic  anharmonic  force 
constants,  respectively.  The  first  two  terms  are  the  har¬ 
monic  Hamiltonian  H0.  The  remaining  terms  are  the 
anharmonic  Hamiltonian  HA,  We  diagonalize  the  har- 


I - - 

monic  Hamiltonian  by  means  of  the  normal-coordinate 
transformation 


T1U,k)= 


31 

2  M.N 


1/2 


2 


WjjUUl 

«v//1 


(2.3) 
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P  (/,«)= -i 


2  N 


1/2 


(2.4) 


and  N  is  the  number  of  unit  cells  in  the  crystal.  The  field 
operators  A~  ^  and  B~  j  are  specified  in  terms  of  the  pho¬ 
non  creation  and  annihilation  operators  b~  .  and  b-,  .  by 
the  relations 

A*.>=b~S J+b*-*.j  ’ 


Here  is  the  normal-mode  frequency  for  wave  vector 


B-  .  =  b _  -bf  .  . 
q q -  q ./ 


(2.5) 

(2.6) 


q  and  branch  index  j,  WU  |  q.j)  is  the  polarization  vector  After  making  the  normal-coordinate  transformation,  the 

for  the  normal  mode,  MK  is  the  mass  of  an  atom  of  type  k,  contributions  to  the  Hamiltonian  take  the  forms 


J 


tfo=7  2*"-  .M-  A-  +Bl  B „  .)=  X  Abl  b~  ,+ j) , 

u  4  q  j  q  ,/  q  ,y  q  ,j  qj  ^  q  ,y  q  ,j  q  tj  * 

q  J  ? 

?./ . 

The  anharmonic  coefficients  F  are  given  by 

13/2 


^(q.y';q'./;q"/')=  ^ 


2# 


(tu_  .0)-,.,  .„)  1/2#A(q+q '  +  q  ”) 
q./  q  ./  q  J 


x2  2  2  4W0,*;/', 

*.a  T  .t’.p  T"X\r 


"D 


1/2 


^q.;;q',y';q"J";q"7"')=^; 


JL 

2N 


(a •  •  •  oi  -  ...  y... )  ~ 1 /2N  A(  q  4-  q  '  +  q  "  -rq  "' ) 


(2.7) 


(2.8) 


(2.9) 


X2  2  2  <t>a^O,K;P,K);P-,K";P",K'") 

*.*  T'J.fi  T"W,t 


Wa[K\q,j)Wg(K>  \q',y)Wr(K"  \q",j")Wb(K"  1  q  "’J"') 
'  (Af,  •  •  •  Mk-)w1 

Xgifq  -R  tn+ 7  "•  rh")+ ?  lin) 


where 


A(q)  = 


1  if  q=G 
0  otherwise 


(2.10) 


and  G  is  a  vector  of  the  reciprocal  lattice. 

In  the  case  of  light  scattering,  the  efficiency  for  Stokes 
scattering  by  zone-center  LO  phonons  in  a  homopolar 
crystal  is  given7  by 


e*LV 


d2S  _ 
dcodSl  2  rrf^m  *a  lMNc  *a> 


o.j 


ft>s 

<u, 


X  (no  1 )  I B  ( j'.I.S )  | 


r(0J;o)) 


[«— n(0,;;o»)]2+rl(0  ,j;a>) 


(2.11) 


where  R  ( j  ;I,S)  is  the  Raman  tensor,  the  branch  index  j 
refers  to  the  longitudinal  optical  branch,  if/  ( rfs > is  the 
wave  vector  of  the  incident  (scattered)  radiation,  a  is  the 
lattice  constant,  L  is  the  crystal  thickness,  id-~  ,  is  the 
zone-center  LO-phonon  frequency,  and  »o  is  the  mean 
number  of  LO  phonons. 

The  resonant  frequency  fl(0, j;co)  in  Eq.  (2.11)  deter¬ 
mines  the  scattering  line  position  and  is  given  to  first  ap¬ 
proximation  by 


fl( 0, ;';&>)  =  cu -j  +  A(0,j;o))  . 


(2.12) 


The  quantities  A(0 ,j\(o)  and  HO ,;';<u)  specify  the  real 
and  imaginary  parts  of  the  proper  self-energy,  BtO, /;&>), 
according  to  the  relation1 


lim  P(0,j ',(!>+ ie)=  —  fHS[M0,j;(o)— iT(0,j;a>)] 
»— o* 


(2.13) 
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where  &  denotes  the  principal  value.  In  Eqs.  (2.14)  and 
(2.15)  we  have  written 


j  ,  i  =1,2,3 
1 


n,  =  n  -  .  = 

1  n  l 


sW. 

e  ''-1 


,  x  =1,2,3 


(2.16a) 

(2.16b) 


where  f}—l/kBT.  The  various  contributions  to  the  fre¬ 
quency  shift  and  damping  constant  are  shown  diagram- 
matically  in  Fig.  1 .  In  addition,  there  are  other  diagrams 
not  shown  in  Fig.  1  which  can  give  nonzero  contributions 
due  to  the  fact  that  the  atoms  in  silicon  do  not  lie  at 
centers  of  inversion  symmetry. 

A  specific  remark  should  be  made  about  the  tempera¬ 
ture  behavior  of  A  and  F.  At  high  temperatures,  i.e.,  tem¬ 
peratures  larger  than  the  Debye  temperature,  the  cubic 
anharmonic  terms  in  A  and  T  given  by  Eqs.  (2.14a)  and 
(2.15a),  respectively,  vary  linearly  with  T.  The  quartic 
anharmonic  term  in  A  corresponding  to  Eq.  (2.14b)  also 
varies  linearly  with  T,  but  the  quadratic  terms  correspond¬ 
ing  to  Eqs.  (2.14c)  and  (2.14d)  vary  quddratically  with  T. 
The  quartic  anharmonic  term  in  T  corresponding  to  Eq. 
(2.15b)  also  varies  quadratically  with  T.  Additional  T2 
contributions  to  both  A  and  T  arise  from  terms  corre¬ 
sponding  to  the  diagrams  in  Fig.  2. 

The  light  scattering  process  can  be  viewed  as  involving 
the  absorption  of  a  photon  fko,,  the  emission  of  a  photon 
fkos,  and  the  creation  of  an  optical  phonon  0 j  which  then 
decays  via  anharmonicity  into  two  phonons,  three  pho¬ 
nons,  etc.  The  production  of  two  and  three  phonons  is 
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(o) 


(b) 


o 
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QUARTIC 


CUBIC 


(Q) 
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FIG.  2.  Diagrams  representing  higher-order  contributions  to 
the  proper  self-energy  of  the  Raman-active  LO  mode  in  silicon. 


shown  diagrammatically  in  Figs.  3(a)  and  3(c).  At 
nonzero  temperatures,  processes  can  also  occur  in  which 
the  decay  of  the  optical  phonon  0 j  is  accompanied  by  the 
absorption  of  another  phonon,  and  the  emission  of  one  or 
more  phonons,  as  shown  in  Figs.  3(b)  and  3(d)  for  the 
cases  of  one  and  two  emitted  phonons,  respectively. 


III.  EXPERIMENTAL  RESULTS 


Light  scattering  measurements  have  been  performed 
with  a  CODERG  PHO  spectrometer  and  an  excitation 
laser  on  single  crystal  nondoped  silicon  with  a  resistivity 
of  100  flcm  and  oriented  with  a  (111)  face  perpendicular 
to  the  incident  beam.  In  view  of  the  large  temperature 
range  explored  the  temperature  was  regulated  in  a  liquid- 
He  cryostat  for  low  temperatures,  an  electrically  heated 
furnace  for  the  intermediate  temperatures,  and  by  laser 
heating  at  high  temperatures. 

The  sample  temperature  was  measured  by  a  platinum 
resistor  for  low  temperatures,  by  a  thermocouple  in  the  in¬ 
termediate  range,  and  by  an  optical  pyrometer  at  high 
temperatures.  Verification  of  the  measured  temperature 
was  made  by  two  additional  methods.  The  first  used  the 
integrated  ratio  of  the  Stokes  to  anti-Stokes  Raman  peaks. 
The  intensity  of  the  Stokes  and  anti-Stokes  peaks  being 
proportional,  respectively,  to  n0+l  and  "o>  the  intensity 
ratio  is 


ftOjQ 

k,T 


(3.1) 


(f)  - 1 - \  QUARTIC  where  a)0  is  the  Raman  frequency.  (We  omit  the  subscript 

^ — o'  j  from  here  on.) 

FIG.  1.  Diagrams  representing  contributions  to  the  frequency  The  second  method  was  based  on  the  black-body  radia- 
shift  A  and  damping  constant  T  for  the  Raman-active  LO  mode  tion  of  the  sample.  If  we  admit  that  the  sample  is  a  black 
in  silicon.  body  we  can  apply  Planck’s  law  for  the  power  emitted  per 
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FIG.  3.  Diagrams  representing  three-  and  four-phonon 
anharmonic  processes  contributing  to  the  decay  of  the  Raman- 
active  LO  mode  in  silicon 
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and  with  aioj)—  1  get  the  sample  temperature. 

The  shape  and  the  position  of  the  peak  due  to  scattering 
by  the  Raman-active  LO  mode  in  silicon  vary  for  different 
temperatures.  In  Fig.  4  we  show  for  comparison  two  spec¬ 
tra  taken  at  295  and  1140  K.  From  these  spectra  we  can 
deduce  the  values  of  T  and  fl  for  these  two  temperatures: 
for  T  =295  K,  r(0,LO)  =  4  cm'1,  and  ft(0,LO)  =  520 
cm~!  and  for  T  =  1 140  K,  f(0,LO)=14  cm"'  and 


FIG.  4.  First-order  Raman  spectra  for  silicon  at  295  and 
1 140  K.  The  hue  pos.'ion  of  the  zone-center  LO  mode  shifts 
from  320  cm~  295  K  to  498  cm-'  at  1140  K  and  the 
linewidth  frvn  4cm'1  at  295  K  to  14cm'1  at  1140  K. 
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FIG.  5.  Anti-Stokes  to  Stokes  intensity  ratio  vs  temperature, 
considering  the  correction  as  discussed  in  the  text.  The  closed 
circles  represent  points  for  which  temperature  was  measured  by 
a  Pt  resistor,  a  thermocouple  and  for  temperatures  above  600’  C 
by  an  optical  pyrometer.  The  open  circles  are  points  obtained  by 
heating  the  sample  with  the  laser  and  the  oven;  their  tempera¬ 
ture  is  not  precisely  determined.  The  theoretical  curve 
expl  ~fki)0/kBT )  is  represented  by  a  solid  lint 


fl(0,LO)=498  cm-1.  The  Stokes  to  anu-»toKes  ratio  of 
the  intensity  of  these  peaks  gives  the  temperature  by  using 
the  solid  line  representation  given  in  Fig  5.  In  this  lie-m 
the  black  circles  represent  the  Pt  resistor,  inermocouole 
and  pyrometric  measurements  of  the  temperature.  !  •’ 

der  to  reach  the  melting  point,  we  heated  th  samp'"  with 
the  oven  and  the  laser,  by  increasing  the  po"  e-  beam 
This  is  represented  by  the  circles.  We  should  non-  th-  for 
these  points  we  were  unable  to  measure  the  tempetature 
accurately.  The  calibration  of  the  solid  lint  applies  alter 
the  measured  Raman  intensities  have  been  correct  e.3 
the  actual  absorption  coefficient  and  the  frequence  depen¬ 
dence  of  the  Raman  efficiency.  In  applying  these  cor-?;- 
tions  the  expression  for  the  intensity  ratio  hecom 
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where  ay,aAs,gj  are  the  absorption  constants  at  the  ire- 
quencies  <y/,euAs,tus  (incident  beam,  anti-Stokes,  and 
Stokes)  and  S(tu/, oi5)  and  S(to;,o>AS)  are  the  Raman  cross 
sections  at  the  involved  frequencies.  Practically  all  the 
points  obtained  by  pyrometric  measurements  are  above  the 
curve  given  by  the  Raman  intensity  ratio.  This  indicates 
that  the  temperature  determined  by  this  method  is  sys¬ 
tematically  higher  than  that  obtained  by  other  measure¬ 
ments.  A  better  knowledge  of  the  correction  factors  is 
therefore  necessary  in  order  for  this  method  to  be  used  tor 
temperature  measurements 

The  damping  constant  and  the  frequency  shift  have 
been  investigated  systematically  as  a  function  of  tempe:  a- 
ture.  Figure  6  gives  the  temperature  variation  of  the 
damping  constant  FtT)  between  5  and  1400  K.  The 
dashed  curve  represents  F(  T I  calculated  from  the  relation5 


nn=no) 


1  + 


(3.4) 


FIG.  6.  Temperature  dependence  of  the  damping  constant  T 
for  the  Raman-active  LO  mode  in  silicon.  The  solid  curve  gives 
the  theoretical  fit  using  both  three-  and  four-phonon  processes. 
The  dashed  curve  gives  the  theoretical  fit  using  only  three- 
phonon  processes.  The  open  and  closed  circles  have  the  same 
significance  as  in  Fig.  5. 


where  x  =tiat0/2kBT  and  r(0)=  1.40  cm-1.  Equation 
(3.4)  is  an  approximate  expression  for  the  temperature 
dependence  of  the  damping  constant  based  on  three- 
phonon  processes  (cubic  anharmonicity  in  second  order) 
and  the  simple  Klemens5  model.  It  seriously  underesti¬ 
mates  the  damping  constant  at  high  temperatures.  We  at¬ 
tribute  this  discrepancy  at  least  in  part  to  the  neglect  of 
four-phonon  processes  associated  with  the  diagrams  in 
Figs.  1(f)  and  2(a)— 2(c). 

It  is  of  interest  to  investigate  whether  this  discrepancy 
can  be  eliminated  by  generalizing  Eq.  (3.4)  to  include  the 
contribution  of  four-phonon  processes.  Following  the  ap¬ 
proach  of  Klemens5  we  write  the  kinetic  equation  for  the 
net  rate  of  decay  of  an  incident  phonon  into  three  thermal 
phonons  in  the  form 

— -(8n0)  =  — &[  (8n0+n0)(« i  -1- 1  )(n 2  -+- 1  )(zi 3 -♦- 1 ) 

at 

—  (8n0+«o+  l)"i”2«3]  .  (3-5) 

where  8n0  is  the  deviation  of  the  incident  phonon  occupa¬ 
tion  number  from  its  thermal  equilibrium  value  n0  and  B 
is  a  constant.  Using  the  equilibrium  condition 

flo(n  1  +  l)(«i  +  l)(»j  +  l) — (no+l)«in2n3=0  > 

(3.6) 

we  can  rewrite  Eq.  (3.5)  as 

— (8n0)=  — Bi  fl|n2+ninj-Fn2fl3  +  ni-t-n2 

at 

+  l)8n0  •  (3.7) 


Energy  conservation  can  be  satisfied  in  the  simple  Kle¬ 
mens  fashion  by  setting  a>x=a)1=<oi=<oo/l.  Consequent¬ 
ly,  U|=n2=n3.  The  generalization  of  Eq.  (3.4)  to  four- 
phonon  processes  then  takes  the  form 


T(T)=A 


3 

(e'-l)2 


(3.8) 


where  y  =fkoQ/ikBT  and  A  and  B  are  constants.  In  the 
high-temperature  limit,  the  factors  multiplying  A  and  B  in 
Eq.  (3.8)  vary  as  T  and  T1,  respectively. 

Equation  (3.8)  has  been  used  to  fit  the  experimental  data 
presented  in  Fig.  6  by  suitably  choosing  the  constants  A 
and  B.  The  best  values  of  A  and  B  are  found  to  be  1.295 
and  0.105  cm-1,  respectively,  and  the  resulting  plot  of 
HD  vs  T  is  given  by  the  solid  curve  in  Fig.  6.  We  see 
that  the  agreement  between  the  calculated  curve  and  the 
experimental  points  is  now  quite  good. 

The  experimental  results  for  the  line  position  ft(  T)  as  a 
function  of  T  are  shown  in  Fig.  7.  Also  shown  is  the  fit  to 
the  data  (solid  curve)  specified  by  the  expressions 

in  n =*><,+ Am  0.9) 


where  ct>o>  C,  and  D  are  constants  with  the  values  528, 
—2.96,  and  —0.174  cm-1,  respectively.  Equation  (3.10)  is 
the  analog  of  Eq.  (3.8)  and  specifies  the  contributions  of 
three-phonon  and  four-phonon  processes  to  the  frequency 
shift.  The  agreement  between  the  experimental  points  and 
the  solid  curve  is  seen  to  be  good. 

If  we  try  to  fit  the  experimental  data  with  three-phonon 
processes  only  by  omitting  the  term  in  Eq.  (3.10)  with  the 
factor  D,  we  obtain  the  dashed  curve  in  Fig.  7  with 
a»o=529  cm-'  and  C  =  —4.24  cm-1.  Although  this 
curve  fits  the  data  well  at  temperatures  up  to  600  K,  it  is 
clearly  inadequate  at  higher  temperatures.  This  demon¬ 
strates  the  necessity  of  including  terms  corresponding  to 
four-phonon  processes  in  the  expression  for  A(  T). 

In  principle,  the  four-phonon  contributions  in  Eqs.  (3.8) 
and  (3.10)  should  include  terms  arising  from  difference 
processes  of  the  type  represented  by  Fig.  3(d).  We  have 
omitted  such  terms  on  the  grounds  that  their  inclusion 


FIG.  7.  Temperature  dependence  of  the  line  position  fl  for 
the  Raman-active  LO  mode  in  silicon.  The  solid  curve  gives  the 
theoretical  fit  using  both  three-  and  four-phonon  processes.  The 
dashed  curve  gives  the  theoretical  fit  using  only  three-phonon 
processes.  The  open  and  closed  circles  have  the  same  signifi¬ 
cance  as  in  Fig.  5. 
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would  Simply  introduce  additional  terms  varying  as  T  and 
as  T 2  in  the  high-temperature  limit  and  would  not  add 
any  new  qualitative  features. 

IV.  DISCUSSION 

We  have  seen  in  the  previous  section  that  the  extension 
of  the  Klemens-Hart-Aggarwal-Lax  model3,5  to  include 
four-phonon  processes  provides  a  good  fit  to  the  experi¬ 
mental  values  for  the  frequency  shift  and  damping  con¬ 
stant  of  the  Raman-active  mode  in  silicon  up  to  1400  K. 
This  fit  is  achieved  by  suitably  choosing  two  adjustable 
parameters  for  each  of  the  two  quantities.  Since  one 
would  expect  the  contribution  of  four-phonon  processes  to 
be  small  compared  to  that  of  three-phonon  processes,  the 
ratios  B/A  and  D/C  should  be  small.  The  actual  values 
of  these  ratios  are  0.08  and  0.06,  respectively  ,  so  this  ex¬ 
pectation  is  fulfilled. 

It  would,  of  course,  be  desirable  to  carry  out  a  first- 
principles  calculation  of  the  frequency  shift  and  damping 
constant  arising  from  both  cubic  and  quartic  anharmonici- 
ty.  However,  such  a  calculation  is  by  no  means  trivial. 
The  principal  difficulty  is  that  a  simple  model  such  as  a 
nearest-neighbor  model  is  inadequate  to  describe  either  the 
harmonic  or  the  anharmonic  properties  of  silicon.  This 
was  shown  a  quarter  century  ago  by  Herman8  and  by  Lax9 
for  the  harmonic  properties  and  very  recently  by  Wanser 


and  Wallis10  for  anharmonic  properties.  Long-range 
forces  are  necessary  for  a  proper  description  and  can  be  in¬ 
troduced  via  a  shell  model,"  a  bond  charge  model,12  or  a 
model  containing  dipole-dipole  and/or  quadrupole- 
quadrupole  interactions.9, 10, 13 

Cowley4  has  carried  out  a  calculation  of  the  frequency 
shift  and  damping  constant  for  silicon  at  temperatures  up 
to  500  K  using  a  shell  model  for  the  harmonic  forces  and 
a  nearest-neighbor  model  for  the  anharmonic  forces. 
Reasonable  agreement  with  the  experimental  data  was  ob¬ 
tained  by  Cowley  for  the  frequency  shifts,  but  not  the 
damping  constant.3  We  are  currently  engaged  in  making  a 
calculation  of  these  quantities  using  the  long-range  force 
model  of  Wanser  and  Wallis  generalized  to  quartic  anhar- 
monicity. 
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1 .  INTRODUCTION 

One  of  the  most  fascinating  debates  in  the  field  of  Semiconductors  Physics 
in  recent  years  .centers  on  the  fundamental  interpretation  of  laser  anneal¬ 
ing.  The  question  of  the  energy  transfer  from  an  intense  beam  to  a  disor¬ 
dered  . material  ,  such  as  amorphous  silicon,  resulting  in  the  crystalliza¬ 
tion  of  the  amorphous  substance  has  been  approached  from  two  different 
points  cf  view,  both  referring  to  a  set  of  fairly  clear  experimental  re¬ 
sults.  On  one  hand,  a  claim  has  been  made  that  the  laser  beam  simply  heats 
the  sample  up  to  melting  the  amorphous  material  which  on  cooling  crystal li 
zes  from  melt  (j-4j  .  On  the  other  hand,  one  has  considered  that  amorphous 
to  crystalline  phase  transition  can  occur  at  low  temperatures  without 
passing  through  the  molten  state  t5J  .  Many  arguments  have  been  dcvelopped 
in  support  of  these  interpretations.  A  l^rge  number  of  publications  have 
appeared  in  the  literature. 

In  addition  to  the  fundamental  interest,  laser  annealing  has  been 
considered  to  have  a  strong  potential  in  the  technology  of  semiconductor 
doping  by  ion  implantation.  For  all  these  reasons,  the  attention  of  a  lar¬ 
ge  audit  nee  has  been  focussed  on  this  problem. 

2.  STRUCTURAL  ANALYSIS  AFTER  LASER  ANNEALING 


2.!.  Amorphous  -  Crystalline  transformation 


The  aim  of  laser  annealing  is  to  achieve  the  transformation  of  the  amor¬ 
phous  material  obtained  by  ion  implantation  into  perfecl  crystal.  There¬ 
fore,  we  shall  first  discuss  the  results  of  the  "'-uctural  analysis  after 
annea ling. 

One  of  the  methods  which  seems  to  be  most  appropriate  to  define  the 
dogi ee  of  crystallinity  of  the  material  is  light  scattering.  The  light  ‘ 
scattering  spectrum  of  amorphous  silicon  for  example,  is  very  different 
firm  that  of  the  perfect  crystal. 

Laser  annealing  performed  with  Ruby  Q-switched  laser  delivering  pulses 
of  a’^out  one  hundred  nanoseconds  with  energy  density  of  the  order  of  one 
J/crr  is  described  by  Morhange  At  this  energy  density,  the  diameter 

of  the  power  laser  beam  was  of  the  order  of  one  mm.  The  annealed  region 
is  explored  with  an  Argon  laser  probe  beam  focussed  to  a  diameter  of  80  pin 
Tin  ’pr  :num  is  used  lor  Raman  spectroscopy.  With  this  arrangement,  it 
was  poss  ’.  U  to  explore  the  annealed  region  in  the  xy  plane  by  moving  the 
sample  will,  regard  to  the  probe  beam  and  in  the  z  direction  by  varying  the 


2 


frequency  of  the  probe  beam  inducing  changes  in  the  penetration  depth  in 
the  a-S'i  sample. 

Two  parameters  characterize  the  annealed  Si  :  the  frequency  position  of 
the  Raman  line  and  its  half  width.  For  monocrystalline  Si,  the  Raman  acti¬ 
ve  normal  mode  is  observed  at  520,5  cm" '  with  a  half  width  of  3  cm-*.  When 

.  o 

ion  implanted  samples  having  a  thickness  of  the  amorphous  layer  of  5000  A 
are  investigated,  the  annealed  region  does  not  show  a  uniform  single 
crystal  structure.  At  the  center  of  the  irradiated  region,  the  Raman  peak 
frequency  is  519  cm-*,  shifted  by  1,5  cm”  towards  lower  frequencies.  The 
width  of  the  peak  is  6  cm-*  significantly  broader  then  the  single  crystal 
peak.  When  the  probe  beam  is  moved  towards  the  periphery  of  the  power 
laser  irradiated  region,  the  frequency  shift  increase  and  reaches  5,5  cm 
at  the  edge  of  the  annealed  area.  The  band  is  also  broader  and  has  a 
width  of  10  cm-*  near  the  edge.  These  results  are  displayed  in  Figure  1, 
where  the  circles  represent  the  extent  of  the  power  laser  beam  and  the 
dots  indicate  the  position  of  the  probe  beam  for  the  spectrum  represented 
on  the  left  side  of  the  figure. 


Figure  1.  Raman  spectra  coming  from 
various  points  along  a 
diameter  of  a  laser  annealed 
region.  From  J.F.  MORHANGE 
et  al  .  Proceedings  of  the 
-  Material  Research  Society 

Annual  Meeting,  Boston, 

1978  (American  Institute 
of  Physics,  New-York  1979) 
p.  429. 
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These  results  could  be  viewed  as  follows.  Even  at  the  center  of  the 
annealed  region,  laser  annealing  does  not  lead  to  a  large  single  crystal 
in  the  way  it  would  be  obtained  in  equilibrium  epitaxial  regrowth. 
Instead , , the  annealed  material  consists  of  large  pol ycristal 1 i tes  whose 
dimensions  decrease  as  one  approaches  the  interface  between  annealed  and 
amorphous  material.  The  recrystallization  occurs  as  a  result  of  random 
nucleation  in  the  amorphous  layer.  The  dimensions  of  the  crystallites 
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depend  on  energy  distribution  in  the  incident  power  laser  beam.  The  crys¬ 
tallite  dimensions  determine  the  frequency  shift  of  the  Raman  active  nor¬ 
mal  mode  band.  An  account  of  the  normal  mode  frequency  shift  as  a  function 
of  the  dimension  of  the  crystallite  is  given  by  a  simple  lattice  dynamics 
calculatibn  developped  by  Kanellis  W  .  The  model  used  in  this  calcula¬ 
tion  is  that  of  a  thin  slab  limited  iii  one  direction  and  infinite  in  the 
plane  perpendicular  to  that  direction.  In  the  case  of  Si,  the  results 
obtained  are  shown  in  Figure  2. 


Figure  2.  Frequency 
variation  of  the  higher- 
frequency  optical  modes- 
as  a  function  of  the 
number  of  cells. 

From  G.  KANELLIS  et  al . , 
Phys.  Rev.  B2J,  1543  (1980) 


In  this  figure  the  frequency  shift  of  the  two  high  frequency  modes  are 
represented.  These  are  the  surface  modes  which  in  the  limit  of  infinite 
crystal  tends  toward  the  Raman  active  mode  at  the  center  of  the  Brillouin 
Zone.  These  calculations  show  that  for  crystallite  having  dimensions 
smaller  than  80  A,  a  noticeable  frequency  shift  should  occur.  Indeed,  the 
experimental  observations  show  that  in  the  region  near  the  crystalline- 
amorphous  interface  a  significant  frequency  Bhift  is  measured.  From  this 
theoretical  model  and  the  experimental  results,  we  draw  the  conclusion 
that  in  laser  annealing,  the  crystallization  occurs  randomly. 

An  amorphous  network  ic  not  a  perfectly  regular  structure.  The  energy 
supplied  by  the  laser  beam  at  threshold  is-  sufficient  to  soften  bonds  in 
the  less  favourable  topological  situation  and  allow  nucieation.  Crystal- 


lographic  ordering  develops  around  the  randomly  distributed  nucleus  embed¬ 
ded  in  the  amorphous  background  creating  microstalline  clusters  whose  size 
depends  on  the  energy- density  distribution. 

First  order  phase  transitions  such  as  melting  or  crystallization  take 
place,  in  general,  via  riucleation  and  growth  of  microphases.  Unstable 
phases  of  Si  have  also  been  prepared  in  high-pressure  experiments  LO. 
Phillips  has  given  evidence  of  metastable  phases  in  laser-induced  and 
thermally . reversible  microcrystallization  in  the  chalcogenide  glass  formers 
GeSe2*  Recent  light  scattering  experiments  ^.1 0^  shows  the  existence  of 
microcrystal 1 ites  embedded  in  the  glass  (GeSe2)  during  laser  annealing.  The 
micro  crystal  line  clusters  initially  formed  are  free  to  rotate  and  form 
larger  clusters,  a  fact  which'is  expressed  by  Raman  line  narrowing  and 
demonstrates  a  precursor  effect  in  the  laser  induced  glass  to  crystal 
transition.  These  experiments  although  achieved  in  different  materials 
and  in  a  very  different  time  scale  might  be  suggestive  of  the. processes 
occuring  in  silicon. 

2.2.  Picosecond  laser  induced  patterns  on  silicon  single  crystal  surface 


It  is  generally  believed  that  monocrystalline  silicon  being  a  well  defined 
target'would  yield  unambiguous  indication  about  the  mechanism  of  the 
transformations  produced  by  a  laser  beam  on  a  semiconductor. 


The  surface  of  a  silicon  single  crystal  may  be  transformed  into  an 
amorphous  layer  by  a  single  picosecond  pulse  Jj  )J  .  Optical  observation 
show  the  formation  of  an  amorphous  ring  pattern  on  a  (111)  crystal 
surface  after  exposure  to  a  laser  pulse  at  532  nm  with  an  average  duration 
of  30  ps  and  with  a  spot  having  a  size  of  the  order  of  5 . 1 0“ ^  cm2.  At 
lower  intensities,  the  ring  diameter  becomes  smaller,  the  amorphous  region 
coalesces  to  a  spot  at  the  center  and  disappears  below  a  critical  intensity 
threshold.  At2532  nm,  the  threshold  is  reported  to  be  0,18  J/cm  for  (100) 
and  0,08  J/cni  for  (111)  surfaces.  Amorphous  rings  rather  than  central 
spots  are  formed  for  intensity  levels  exceeding  0,24  J/cm2  for  (100)  and 
0,12  J/cm2  for  (1)1)  surfaces.  The  center  in  this  case  is  a  single  crystal 
with  the  same  orientation  as  the  substrate.  The  amorphous  nature  of  the 
rings  was  deduced  from  an  inr ease  of  the  reflectivity  compared  to  crystal¬ 
line  silicon  and  was  confirmed  by  electron  diffraction  with  a  transmission 
electron  microscope. 


The  annular  pattern  of  the  amorphous  phase  is  interpreted  in  terms  of 
cooling  rate  and  crystal-growth  speed.  It  is  first  assumed  that  above 
threshold,  the  laser  pulse  induces  melting  of  the  silicon  ta.rget. 


Following  laser  pulse  melting,  two  limiting  cases  are  then  proposed  : 
just  above  the  threshold  intensity,  a  layer  of  depth  somewhat  smaller  than 
5.10  cm  is  cooled  in  t  =  3.1 0—  s  and  the  condition,  d/U^  >t  ,  for 
nucleatior.  of  an  amorphoBs  solid  phase,  is  fulfilled.  Well  above  the  thres¬ 
hold  intensity,  the  temperature  of  the  molten  layer  rises  far  above  T 
and  the  cooling  rate  in  the  critical  region  below  0.8  T^  becomes  longer. 
This  gives  a  description  of  the  annular  pattern  of  the  amorphous  phase. 


2.3.  Picosecond  laser  annealing  of  implanted  silicon 
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implanted  amorphous  silicon. 

With  a  laser  beam  of  I .06  ym,  beside  the  annular  amorphous  region  and 
recrystall izcd  center,  one  also  observes  a  re'crystallized  ring.  Both  the 
ring  and  the  center  are  polycrystals.  The  energy  fluence  level  to  form^the 
ring  is  0.22  J/cm  .  It  is  0.35  J/cnt  for  the  a-Si  region  and  0.85  J/cm 
or  the  recrystallized  region. 


Rozgonyi  et  al  £l3j  have  re-examined  the  structural  modification  of 
amorphized  silicon  surfaces  following  picosecond  laser  irradiation  using 
cross-section  transmission  electron  microscopy  and  showed  that  the  centcr 
is  always  a  dislocation  free  single  crystal  encircled  by  poly-silicon; 
ring.  This  observation  is  consistent  with  the  results  of  Morhange  et  al 

4 J  obtained  by  Raman  spectroscopy. 

More  recently,  Nissim  et  al  jj5j  have  performed  light  scattering  measu¬ 
rements  analogous  to  that  reported  by  Morhange  et  al.  £l4j  scanning  the 
irradiated  area  with  a  probe  laser  beam  focussed  on  a  I  ym  diameter  spot. 
When  a  single  30ps  pulse  at  1.06  and  0.532  ym  wavelengths  from  a  mode- 
locked  neodymium  :  Yttrium  Aluminium  Garnet  laser  is  used,  a  multi  annular 
recrystallization  pattern  is  observed  on  implanted  silicon.  At  high 
incident  energies,  single  crystal  silicon  is  observed  in  the  central  spot 
and  in  the  first  recrystallized  ring  of  the  annealed  area.  With  an  irra¬ 
diation  at  1^06  ym,  the  threshold  of  laser  induced  damage  was  found  to  be 
above  2  J/cm  ,  The  multi  annular  pattern  has  been  ascribed  to  multiple 
melting  resolidification  process  during  the  pulse  duration  leading  to  the 
superposition  of  basic  structures  of  different  sizes. 


For  the  undet standing  of  this  complex  structures,  it  should  also  be 
remembered  that  periodic  surface  structures  on  solids  may  result  from 
inhomogeneous  energy  deposition  associated  with  the  interference  of  the 


incident  beam  with  a  surface  scattered  field  as  discussed  by  Van  Driel 
et  al  {j6  ]  . 


3.  TIME  RESOLVED  ANALYSIS  OF  PULSED  LASER  IRRADIATION  OF  SEMICONDUCTORS 


For  the  discussion  of  the  data  obtained  by  short  laser  irradiation  pulses 
one  could  distinguish  three  different  time  ranges  :  i)  very  short  time 
scale,  t  <<  10-^  s,  where  the  system  is  in  a  far  from  equilibrium  state,, 
the  seat  of  highly  non  linear  processes,  ii)  intermediate  time  scale, 
t  "v  10-1  s  when  different  elementary  relaxation  processes  take  place  and 
iii)  long  time  range  t  >  10”  s  when  thermalization  occure  and  the 
system  tends  toward  equilibrium. 

V 

3.1.  Investigations  in  the  very  short  time  range 

2 

A  radiation  beam,  of  energy  density  of  i  J/cm  ,  is  bound  to  produce  a 
strong  perturbation  on  the  material  at  the  instant  of  interaction.  The 
photon  density  is  such  that  a  hot  electron  plasma  might  be  produced. 

The  experiments  in  the  shortest  time  scale  are  these  of  optical  pulse 
induced  phase  transitions  in  silicon  described  by  Shank,  Yen  and  Hirlimann 
Ll73  .  They  reported  the  first  observations  of  optically  induced  reflec¬ 
tivity  changes  in  silicon  with  90  fs  optical  pulses. 

The  results  of  reflectivity  measurements  are  shown  in  Figure  3.  The 
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energy  E  is  defined  as  the  excitation  energy  density  where  visual  evi¬ 
dence  of  amorphous  layer  formation  is  observed.  This  energy  corresponds 
to  0. 1  J/cm^. 


Figure  3.  Transient  reflectivi¬ 
ty  data  in  Silicon  at  probe 
wavelength  I  pm  and  pump 
wavelength  620  nm  at  various 
incident  energies. 

From  C.V.  SHANK  et  al . ,  Phys. 
Rev.  Letters  50,  454  (1983). 
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Shank  et  al  no  discuss  the  results  shown  in  Figure  3  in  the  follo¬ 
wing  way.  Optical  excitation  of  Si  with  a  2.0  eV  optical  pulse  results 
in  the  generation  of_a  dense  electron-hole  plasma  within  the  optical 
absorption  depth  :  a  1  =  3  pm.  They  suppose  that  during  the  earlier  time 
following  excitation,  the  reflectivity  is  dominated  by  the  electron-hole 
plasma.  As  the  pulse  intensity  is  increased,  the  energy  is  transfered  to 
the  crystal  lattice  and  the  crystal  melts.  It  is  also  supposed  that 
melting  begins  at  the  surface  and  moves  inward  into  the  bulk. 

Probably,  the  most  important  contribution  of  this  investigation  is  the 
demonstration  that  when  the  crystal  is  excited  with  a  short  optical  pulse 
of  90  fs  a  form  of  unstable  highly  excited  silicon  is  created  which  per¬ 
sists  for  a  fraction  of  picosecond.  A  challenging  programme  now  is  to 
determine  the  properties  of.  this  material  with  at  least  10  %  of  the  avai¬ 
lable  electrons  excited. 

3.2.  Intermediate  time  range.  Picosecond  irradiation 

Most  of  the  wort  on  time  resolved  spectroscopy  is  in  the  range  of  the 
picosecond  for  the  simple  reason  that  tins  is  the  time  scale  for  which 
short  pulse  laser  sources  are  readily  available. 

Experiments  in  the  picosecond  range  on  transmission  and  reflectivity 
were  performed  by  l,iu  et  al.  8j  at  the  fundamental  and  doubled  frequency 
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of  a  mode-locked  Nd  YAG  laser  producing  a  30  ps  single  pulse.  In  this 
vork  are  presented  measurements  on  self  reflectivity  and  self-transmission 
at  X  =  532  nm  for  increasing  energy  fluence  with  20.ps  pulses.  The  self¬ 
reflectivity  of  bulk  silicon  with  (111)  surface  starts  to  increase  when  the 
energy  fluence  of  the  pump  beam  reaches  0.2  J/cm  .  The  initial  rise  of  the 
reflectivity  from  the  crystalline  value  of  0.37  to  a  maximum  value  of  0.48 
is  in  agreement  with  changes  in  the  index  of  refraction  due  to  melting  of  a 
thin  -surface  layer.  In  this  case,  it  is  supposed  that  the  hot  e-h  plasma 
transfers  sufficient  energy  for  melting  within  the  duration  of  the  pulse 
iltself.  Below  the  critical  fluence  level  of  0.2  J/cm  ,  the  photoexcited 
e-h  plasma  causes  a  decrease  in  the  real  part  n  of  the  complex  refractive 
index. 

A  pump  and  probe  technique  was  also  used  where  a  picosecond  excitation 
pulse  is  followed  by  a  slightly  focussed  weak  probe  with  a  variable  time 
delay. 

At  100  ps  delay,  the  probe  pulse  is  temporarily  completely  separated 
from  the  pump  pulse2and  an  abrupt  rise  in  reflectivity  at  the  threshold 
fluence  of  0.2  J/cm  is  observed  DO  .  This  discontinuity  in  the  car¬ 
rier  density  should  be  associated  with  local  structural  changes.  The 
reflectivity  rises  to  0.75  +  0.03  which  is  characteristic  of  molten  silicon 
at  the  probe  wavelength  of  1.064  um.  This  behaviour  is  observed  even  at 
zero  delay  indicating  that  melting  occurs  within  the  pulse  duration  of 
20  ps. 

3.3.  Long  time  phenomena 

At  long  time,  the  energy  is  transferred  to  the  lattice  which  eventually 
melts.  For  the  long  time  phenomena,  we  are  faced  with  three  major  problems 

i)  what  is  the  mechanism  of  melting, 

ii)  do  intermediate,  metastable  phase  exist  in  the  further  thermalization 

iii)  what  is  the  microscopic  mechanism  of  resolidification. 

There  are  two  experimental  approaches  to  these  questions  : 

a)  determination  of  the  temperature  of  the  system, 

b)  direct  observation  of  the  structural  phase  transf ormation  during  the 
evolution  of  the  system  towards  thermalization. 

We  have  discussed  in  paragraph  2  the  structural  observations  and  the 
insight  they  could  be  for  the  dynamical  evolution  of  the  system.  We  shall 
now  focuss  on  considerations  on  the  determination  of  the  temperature  of  the 
system  :  the  electron  gas  temperature  and  the  lattice  temperature. 

3.4.  Energy  transfert  and  carrier  density 

Of  particular  interest  for  the  understanding  of  the  laser  action  on  a  solid 
is  the  analysis  of  the  initial  events  following  the  irradiation  by  an  inten 
so  optical  pulse.  Femtosecond  spectroscopy  has  already  shown  to  be  a 
valuable  method  for  such  an  analysis  £.17,  19,  20j  .  Measurements  on  the 
time  dependent  reflectivity  leads  to  estimation  of  the  electron  hole 
density,  initially  created  as  well  as  its  evolution  with  time.  Such 
measurements  are  also  suggestive  for  models  for  the  energy  transfer  during 
and  after  the  irradiation. 


The  questions  to  be  considered  are  the  following  : 


i)  what  is  the  electron-hole  (e-h)  density  resulting  from  the  excitation", 

ii)  what  are  the  interaction  processes  behind  the  excitation, 

iii)  how  is  the  energy  stored  in  the  excited  carriers  and  transferred  to 
the  rest;  of  the  system. 

i)  electron-hole  density  produced  f emtosecondjgulse 

The  high  density  electron-hole  plasma  in  silicon  created  by  a  90  femtosecond 
pulse  has  been  investigated  £ 1 7  )  by  measuring  the  time  dependent  reflectivi¬ 
ty  over  a  20  picoseconds  time  scale,  at  various  laser  energies.  Such  short 
pul.ses  lead  to  the  possibility  to  break  so  many  covalent  bonds  that  the 
crystal  becomes  fluid  even  at  T  =  0  K,.  The  carrier  ^|nsity  created  at  this 
time  scale  was  estimated  to  be  of  the  order  of  5.10z  cm~’  under  an  inci¬ 
dent  energy  density  of  0.063  J/cm^.  The  phenomena  produced  in  such  very 
short  pulses  can  be  viewed  in  the  following  way  [20  3  .The  part  of  the  ini¬ 
tial  laser  beam  absorbed  over  a  penetration  depth  of  d  =  3  pm  creates  an 
electron  hole  plasma  with  a  decreasing  density  profile.  After  the  surface 
density  reaches  a  fraction  of  N  ,  the  plasma  density,  the  refli  tivity  and 
the  penetration  depth  fall,  becEuse  of  the  decrease  of  the  real  part  of 
the  dielectric  constant  as  well  as  of  the  increase  of  tha  induced  free 
carrier  absorption.  For  high  enough  power,  the  surface  density  goes  beyond 
N  within  the  pulse,  yielding  an  instantaneous  reflectivity  larger  than  RQ 
tRe  initially  reflected  beam.  At  that  high  surface  density  the  laser  beam 
becomes  a  vanishing  wave  which  can  create  e-h  pairs  only  very  close  to  the 
surface. 

The  reflectivity  during  the  pulse  R(t)  is  obtained  through  a  general 
resolution  of  the  Maxwell  equations,  taking  into  account  the  variation  of 
the  dielectric  constant  via  the  modification  of  the  e-h  density  profile. 

The  experimental  data  [20}  are  fitted  with  a  relaxation  time  x  =  3. 

‘  '  high  as  2.10 

total  number  of 


ii)  electron-gastemgerature 

An  interesting  investigation  on  the  photoexcited  electron  distribution  is 
performed  by  non-linear  photoemission  from  silicon  by  Bensoussan  and 
Moison  Ql]  .  They  show  that  in  Si,  at  moderate  fluences  an  equilibrium 
distribution  coexists  with  the  electrons  at  very  high  energies  in  the 
conduction  band,  generated  by  two  and  three  quantum  processes.  The 
equilibrium  distribution  is  described  by  a  well-defined  temperature  which 
differs  significantly  from  the  lattice  temperature.  Because  of  the  fast 
photogeneration  of  the  carriers,  and  of  the  two  photon  absorption  and 
biparticle  Auger  recombination  processes  in  which  carriers  are  sent 
continuously  high  in  the  bands,  the  electron  gas  can  reach  an  internal 
equilibrium  characterized  by  a  temperature  Tg  higher  than  the  lattice 
temperature  T  .  The  experimentally  measured  thermal  emission  is  interpreted 
by  the  Richardson  equation.  The  temperatures  deduced  for  different  fluences 
are  shown  in  Figure  A. 


s.  For  A 
N  <  8,1 
valence 


A  22°,J  1)1,1  the  Piasma  frequency  density  can  be  as 
,  1CT  cm~3  which  is  indeed  very  high  in  view  of  the 
e  band  electrons  (2.1 02^  cm~^) . 
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Figure  4.  Electron  and  lat¬ 
tice  temperatures  vs.  photon 
flux . 

From  M.  BENSOUSSAN.  Procee¬ 
dings  of  the  16th  Internation 
al  Conference  on  the  Physics 
of  Semiconductors.  North- 
Holland ,  Editor  M.  Averous 
(1982)  p.'  405. 
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iii)  latticc_temgerature 

Lattice  temperature  at  equilibrium  can  be  inferred  from  light  scattering 
measurements  of  the  frequency  shift  and  line  width  of  the  L0  mode  22 
and  Stokes  to  anti-Stokes  ratio  [23]  .  The  Stokes  to  anti-Stokes  ratio 
indicates  a  lattice  temperature  of  1400  K  for  a  delay  of  150  ns  after 
excitation  with  a  power  beam  of  0.8  J/crrr  at  532  nm.  Raman  measurements 
are  evidently  taken  after  the  high  reflectivity  falls  off. 

4.  DISCUSSION  AND  CONCLUSIONS 

Laser  annealing  has  now  a  long  history  and  an  abundant  literature.  The 
interesting  question,  which  still  remains  to  be  answered,  is  what  are  the 
elementary  processes  in  the  interaction  of  a  strong  radiation  field  with 
matter. 

The  process  of  laser  annealing  consists  of  two  sets  of  phenomena.  The 
first  concerns  the  effect  of  a  dense  radiation  field  on  matter  creating 
elementary  excitation  far  from  equilibrium.  The  initial  hot  plasma  redis¬ 
tributes  through  interactions  between  carriers  and  ultimately  thermalizes 
giving  up  energy  to  heavier  particles  which  are  the  lattice  constituents. 
The  second  set  of  phenomena,  of  a  completely  different  nature,  concerns 
the  modifications  induced  to  the  solid  as  a  result  of  the  creation  of  the 
dense  hot  plasma.  If  the  laser  pulse  is  short  enough,  the  e-h  density 
reached  corresponds  to  an  amount  of  broken  covalent  bonds  which  is  a 
significant  portion  of  the  total  number  of  bonds,  the  crystal  becomes 
fluid  even  at  T  =  0  K.  Ultimately,  the  solid  sets  ;i n  a  new  phase  whose 
structure  depends  on  the  elementary  mechanisms  of  interaction  and  organi- 


/  AD-A146  979 
UNCLASSIFIED 


LASER  ANNEALING  OF  ION  IMPLANTED  SEMICONDUCTORS(U) 
PARIS-6  UNI V  (FRANCE)  J  F  MORHANGE  ET  AL.  JUL  84 
DAJA37-81-C-0587 

F/G  20/12 


2/2 


NL 


I 


zation  of  the  lattice  constituents.  We  shall  examine  successively  these 
two  states  as  :  i)  direct  laser  effect  and  ii)  consequences  of  the  laser 


•action. 


4.1.  Direct  laser  effect 

The  incident  laser  energy  is  absorbed  by  electron-hole  pairs  creation  and 
by  free  carrier  excitation.  In  indirect  gap  semiconductors,  electron-hole 
pairs  are  created  via  indirect  absorption  processes  involving  the  emission 
and  absorption  of  phonons.  Because  phonon  energies  are  much  smaller  than 
photon  energy,  the  amount  of  energy  transferred  to  the  lattice  during  j 
absorption  is  negligible  in  comparison  to  the  total  amount  absorbed.  The 
rise  of  carrier  density  leads,  in  turn,  to  increased  free  carrier  absorp¬ 
tion.  The  net  result  is  the  production  of  hot-electrons  and  holes  far  from 
equilibrium  which  subsequently  thermalize  with  the  set  of  the  carriers  and 
eventually  with  the  lattice.  The  observations  by  Shank  et  al  Ql  demons¬ 
trate  that  with  short,  90  femtoseconds  optical  pulse,  an  unstable  form  of 
highly  excited  state  is  created  near  the  surface  which  persists  for  a 
fraction  of  picosecond. 

For  very  short  impulsion  of  radiation,  from  the  experiments  C”.  .  2qJ 
one  perceives  clear  indications  that  there  are  two  steps  of  the  laser 
action  on  a  crystal  well  separated  in  time.  In  the  first  step,  the  electro¬ 
magnetic  energy  transferred  from  the  laser  beam  to  the  solid  is  retained 
in  the  highly  excited  non-equilibrium  electronic  state.  In  a  second  step, 
this  energy  is  transferred  to  the  lattice. 

Recently,  a  theoretical  model  D°3  is  proposed  which  takes  into  account 
the  space  time  evolution  of  the  plasma  during  the  pulse  in  order  to  explain 
the  processes  in  densities  higher  than  I02  cm”  reached  in  100  fs  pulses. 
The  novelty  in  this  model  is  the  argument  that  free  carrier  absorption  is 
dominated  by  e-h  collisions  with  a  characteristic  relaxation  time  t  ** 
3.10"16  s. 


A  different  situation  is  reached  with  much  longer  pulses  :  t,  1 
and  a  photon  absorption  rate  g  “vlOJ  cm-^  s”  :  Auger  recombination 
becomes  the  dominant  recombination  mechanism  at  these  densities. 
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Most  of  the  laser  energy  is  absorbed  by  the  carriers  within  the  absorp¬ 
tion  depth.  Eventually,  these  carriers  loose  their  energy  to  the  lattice, 
the  rise  of  the  lattice  temperature  then  depends  on  the  distance  they  have 
diffused  before  substantial  phonon  emission  occurs.  At  moderate  density, 
the  phonon  scattering  time  is  t  ^  10-*^  s.  Screening  does  not  affect  the 
rate  of  intervalley  phonon  emission  until  Ne  'v  10^'  cm-^.  Because  screening 
increases  the  electron-phonon  scattering  time,  it  not  only  decreases  the 
rate  of  phonon  emission  but  also  enhances  diffusion.  This  increases  the 
volume  of. the  region  in  which  the  energy  of  the  excited  carriers  is  trans¬ 
ferred  to  the  lattice.  Owing  to  the  extreme  non  linearity  of  the  hot 
carrier  effects,  it  is  impossible  to  make  an  accurate  estimate  of  the 
precise  temperature  to  which  the  lattice  is  heated  or  to  determine  the 
.laser  power  threshold  above  which  melting  will  occur. 

4.2.  Con'sequenees  of  the  laser  action 

The  equilibrium  observations  are  clear  :  an  amorphous  or  glass  solid  is 
transformed  into  crystal  under  laser  action  and  a  crystal  submitted  to 
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yery  short  laser  irradiation  is  transformed  into  amorphous  material.  An 
implication  of  both  of  these  transformations  is  that  melting  preceeds  the 
transition.  Another  alternative  is  that  the  phase  transition  is  directly 
induced ’in  the  highly  excited  state. 

The  effect  of  a  dense  plasma  on  the  melting  temperature-  is  itself  an 
interesting  problem  of  solid  state  theory.  This  question  has  been  recently 
addressed  by  Bok  and  Combescot  [24 J .  It  is  shown  that  in  the  presence  of 
a  dense  plasma,  the  melting  temperature  of  a  -solid  changes.  The  melting 
temperature  decreases  with  increasing  plasma  density.  For  a  laser  pulse  of 
J  J/cm  during  10  ns,  it  is  considered  that  e-h  plasma  reaches  a  steady 
state  in  a  time  shorter  than  the  laser  pulse.  Considering  the  plasma 
expansion  due  to  its  high  pressure,  i.ts  collision  with  phono^  and  Auger 
recombination,  the  highest  plasma  density  is  of  the  order  10^  cm  .  This 
density  is  nevertheless  considered  to  be  sufficient  to  considerably  reduce 
the  melting  temperature  so  that  a  metallic  layer  of  liquid  silicon  is  formed 
at  the  surface. 

The  role  of  a  high  plasma  density  in  laser  annealing  has  been  discussed 
by  Van  Vechjten  et  al .  in  a  quite  different  way  but  still  involving 
electron-phonon  coupling  and  lattice  instability  induced  by  this  interac¬ 
tion.  Above  a  critical  carrier  density  estimated  at  8.10^  cm”  ,  a  second 
order  phase  transition  occurs.  At  this  plasma  density,  the  bond  charges 
will  be  so  depleted  that  they  will  no  longer  be  able  to  stabilize  the  TA 
phonon  modes  .  The  crystal  will  no  longer  resist  shearing  stresses 

and  will  become  fluid.  This  fluid  is  distinct  from  the  normal  molten  phase 
of  Si  the  latter  beeing  the  result  of  a  strictly  first  order  phase 
transition  driven  by  the  atomic  motion  at  high  temperatures.  The  assumption 
of  Van  Vechten  is  that  the  plasma  is  supposed  to  directly  induce  the 
structural  transformation.  The  energy  is  retained  in  the  electronic  system 
instead  of  beeing  entirely  associate  with  the  atomic  motion.  As  the  plasma 
becomes  less  dense  due  to  expansion  and  to  transfer  of  energy  to  the 
lattice,  the  material  will  pass  back  through  the  second  order  phase 
transition  at  8.1(r  cm--*  and  covalent  bonding  will  gradually  appear.  The 

material  will  finally  recrystallize  if  this  process  is  relatively  slow 
or  will  solidify  in  the  amorphous  phase  if  the  process  is  very  fast. 

This  dense  plasma  phase  could  be  compared  to  the  highly  excited  silicon 
which  persists  for  a  fraction  of  a  picosecond  0  7,  20]  .  The  interpretation 
of  the  laser  action  differs  nevertheless  with  regard  to  the  following  step  ; 
it  is  generally  considered  that  the  solid  melts  after  the  initial  inter¬ 
action  stage. 

In  conclusion,  few  points  appear  clear  today.  The  laser  interaction 
with  the  solids  results  first  in  the  creation  of  a  dense  plasma  which 
persists  for  a  fraction  of  a  picosecond.  Melting  seems  to  occur  after  the 
excitation  pulse.  The  mechanism  of  melting  is  not  clear  and  consequently 
the  mechanism  of  solidification  is  not  clear  either.  Further  investigation 
in  the  very  short  impulse  regime  ard  certainly  desirable  to  clarify  the 
physical  processes  in  laser  annealing. 
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